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The phenomenon of Bose-Einstein condensation of dilute gases in traps is reviewed from a the- 
oretical perspective. Mean-field theory provides a framework to understand the main features of 
the condensation and the role of interactions between particles. Various properties of these systems 
are discussed, including the density profiles and the energy of the ground state configurations, the 
collective oscillations and the dynamics of the expansion, the condensate fraction and the thermo- 
00 ' dynamic functions. The thermodynamic limit exhibits a scaling behavior in the relevant length and 

On ' energy scales. Despite the dilute nature of the gases, interactions profoundly modify the static as 

well as the dynamic properties of the system; the predictions of mean-field theory are in excellent 
agreement with available experimental results. Effects of superfluidity including the existence of 
quantized vortices and the reduction of the moment of inertia are discussed, as well as the conse- 
quences of coherence such as the Josephson effect and interference phenomena. The review also 
assesses the accuracy and limitations of the mean-field approach. 
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CO \ Bose-Einstein condensation (BEG) (Bose, 1924; Einstein, 1924) was observed in 1995 in a remarkable series of 
■ experiments on vapors of rubidium (Anderson et ai, 1995) and sodium (Davis et at, 1995) in which the atoms were 
confined in magnetic traps and cooled down to extremely low temperatures, of the order of fractions of microkelvins. 
The first evidence for condensation emerged from time of flight measurements. The atoms were left to expand by 
switching off the confining trap and then imaged with optical methods. A sharp peak in the velocity distribution was 
then observed below a certain critical temperature, providing a clear signature for BEG. In Fig. |l|, we show one of the 
"j^ first pictures of the atomic clouds of rubidium. In the same year, first signatures of the occurrence of BEG in vapors 
of lithium were also reported (Bradley et ai, 1995). 
I ■ Though the experiments of 1995 on the alkalis should be considered a milestone in the history of BEG, the experi- 
' mental and theoretical research on this unique phenomenon predicted by quantum statistical mechanics is much older 
C ,: and has involved different areas of physics (for an interdisciplinary review of BEG see Griffin, Snoke and Stringari, 
1995). In particular, from the very beginning, superfluidity in helium was considered by London (1938) as a possible 
manifestation of BEG. Evidences for BEG in helium have later emerged from the analysis of the momentum distri- 
bution of the atoms measured in neutron scattering experiments (Sokol, 1995). In recent years, BEG has been also 
investigated in the gas of paraexcitons in semiconductors (see Wolfe, Lin and Snoke, 1995, and references therein), 
5^] ' but an unambiguous signature for BEG in this system has proven difficult to find. 

I Efforts to Bose condense atomic gases began with hydrogen more than 15 years ago. In a series of experiments 
hydrogen atoms were first cooled in a dilution refrigerator, then trapped by a magnetic field and further cooled by 
evaporation. This approach has come very close to observing BEG, but is still limited by recombination of individual 
atoms to form molecules (Silvera and Walraven, 1980 and 1986; Greytak and Kleppner, 1984; Greytak, 1995; Silvera, 
1995). At the time of this review, first observations of BEG in spin polarized hydrogen have been reported (Fried et 
ai, 1998). In the '80s laser-based techniques, such as laser cooling and magneto-optical trapping, were developed to 
cool and trap neutral atoms [for recent reviews, see Ghu (1998), Gohen-Tannoudji (1998) and Phillips (1998)]. Alkali 
atoms are well suited to laser-based methods because their optical transitions can be excited by available lasers and 
because they have a favourable internal energy-level structure for cooling to very low temperatures. Once they are 
trapped, their temperature can be lowered further by evaporative cooling [this technique has been recently reviewed 
by Ketterle and van Druten (1996a) and by Walraven (1996)]. By combining laser and evaporative cooling for alkali 
atoms, experimentalists eventually succeeded in reaching the temperatures and densities required to observe BEG. 
It is worth noticing that, in these conditions, the equilibrium configuration of the system would be the solid phase. 
Thus, in order to observe BEG, one has to preserve the system in a metastable gas phase for a sufficiently long time. 
This is possible because three-body collisions are rare events in dilute and cold gases, whose lifetime is hence long 
enough to carry out experiments. So far BEG has been realized in ^^Rb (Anderson et ai, 1995; Han et ai, 1998; 
Kasevich, 1997; Ernst et ai, 1998a; Esslinger et at, 1998; Dalibard et ai, 1998), in ^^Na. (Davis et at, 1995; Hau, 
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1997 and 1998; Lutwak et al, 1998) and in '^Li (Bradley et al, 1995 and 1997). The number of experiments on BEC 
in vapors of rubidium and sodium is now growing fast. In the meanwhile, intense experimental research is currently 
carried out also on vapors of caesium, potassium and metastable helium. 

One of the most relevant features of these trapped Bose gases is that they are inhomogeneous and finite-sized 
systems, the number of atoms ranging typically from a few thousands to several millions. In most cases, the confining 
traps are well approximated by harmonic potentials. The trapping frequency, cjho, provides also a characteristic length 
scale for the system, aho = [?i/('TKjJho)]^''^, of the order of a few microns in the available samples. Density variations 
occur on this scale. This is a major difference with respect to other systems, like for instance superfluid helium, 
where the effects of inhomogeneity take place on a microscopic scale fixed by the interatomic distance. In the case 
of *^Rb and ^'^Na, the size of the system is enlarged as an effect of repulsive two-body forces and the trapped gases 
can become almost macroscopic objects, directly measurable with optical methods. As an example, we show in Fig. |^ 
a sequence of "in situ" images of an oscillating condensate of sodium atoms taken at the Massachusetts Institute of 
Technology (MIT), where the mean axial extent is of the order of 0.3 mm. 

The fact that these gases are highly inhomogeneous has several important consequences. First BEC shows up not 
only in momentum space, as happens in superfluid helium, but also in co-ordinate space. This double possibility of 
investigating the effects of condensation is very interesting from both the theoretical and experimental viewpoints and 
provides novel methods of investigation for relevant quantities, like the temperature dependence of the condensate, 
energy and density distributions, interference phenomena, frequencies of collective excitations, and so on. 

Another important consequence of the inhomogeneity of these systems is the role played by two-body interactions. 
This aspect will be extensively discussed in the present review. The main point is that, despite the very dilute 
nature of these gases (typically the average distance between atoms is more than ten times the range of interatomic 
forces), the combination of BEC and harmonic trapping greatly enhances the effects of the atom-atom interactions 
on important measurable quantities. For instance, the central density of the interacting gas at very low temperature 
can be easily one or two orders of magnitude smaller than the density predicted for an ideal gas in the same trap, as 
shown in Fig. |3[ Despite the inhomogeneity of these systems, which makes the solution of the many-body problem 
nontrivial, the dilute nature of the gas allows one to describe the effects of the interaction in a rather fundamental way. 
In practice a single physical parameter, the s-wave scattering length, is sufficient to obtain an accurate description. 

The recent experimental achievements of BEC in alkali vapors have renewed a great interest in the theoretical 
studies of Bose gases. A rather massive amount of work has been done in the last couple of years, both to interpret 
the initial observations and to predict new phenomena. In the presence of harmonic confinement, the many-body 
theory of interacting Bose gases gives rise to several unexpected features. This opens new theoretical perspectives 
in this interdisciplinary field, where useful concepts coming from different areas of physics (atomic physics, quantum 
optics, statistical mechanics and condensed matter physics) are now merging together. 

The natural starting point for studying the behavior of these systems is the theory of weakly interacting bosons 
which, for inhomogeneous systems, takes the form of the Gross-Pitaevskii theory. This is a mean-field approach for 
the order parameter associated with the condensate. It provides closed and relatively simple equations for describing 
the relevant phenomena associated with BEC. In particular, it reproduces typical properties exhibited by superfiuid 
systems, like the propagation of collective excitations and the interference effects originating from the phase of the 
order parameter. The theory is well suited to describing most of the effects of two-body interactions in these dilute 
gases at zero temperature and can be naturally generalized to explore also thermal effects. 

An extensive discussion of the application of mean-field theory to these systems is the main basis of the present 
review article. We also give, whenever possible, simple arguments based on scales of length, energy and density, in 
order to point out the relevant parameters for the description of the various phenomena. 

There are several topics which are only marginally discussed in our paper. These include, among others, coUisional 
and thermalization processes, phase diffusion phenomena, light scattering from the condensate and analogies with 
systems of coherent photons. In this sense our work is complementary to other recent review articles (Burnett, 1996; 
Parkins and Walls, 1998). Furthermore in our paper we do not discuss the physics of ultracold collisions and the 
determination of the scattering length which have been recently the object of important experimental and theoretical 
studies in the alkalis (Heinzen, 1997; Weiner et al., 1998). 

The plan of the paper is the following: 

In Sec. H we summarize the basic features of the noninteracting Bose gas in harmonic traps and we introduce the 
first relevant length and energy scales, like the oscillator length and the critical temperature. We also comment on 
finite size effects, on the role of dimensionality and on the possible relevance of anharmonic traps. 

In Sec. Ill we discuss the effects of the interaction on the ground state. We develop the formalism of mean-field 
theory, based on the Gross-Pitaevskii equation. We consider the case of gases interacting with both repulsive and 
attractive forces. We then discuss in detail the large N limit for systems interacting with repulsive forces, leading to 
the so called Thomas-Fermi approximation, where the ground state properties can be calculated in analytic form. In 
the last part, we discuss the validity of the mean-field approach and give explicit results for the first corrections, beyond 
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mean- field, to the ground state properties, including the quantum depletion of the condensate, i.e., the decrease in 
the condensate fraction produced by the interaction. 

In Sec. IV we investigate the dynamic behavior of the condensate using the time dependent Gross-Pitaevskii 
equation. The equations of motion for the density and the velocity field of the condensate in the large N limit, where 
the Thomas-Fermi approximation is valid, are shown to have the form of the hydrodynamic equations of supcrfluids. 
We also discuss the dynamic behavior in the nonlinear regime (large amplitude oscillations and free expansion), the 
collective modes in the case of attractive forces and the transition from collective to single-particle states in the 
spectrum of excitations. 

In Sec. we discuss thermal effects. We show how one can define the thermodynamic limit in these inhomogeneous 
systems and how interactions modify the behavior compared to the noninteracting case. We extensively discuss the 
occurrence of scaling properties in the thermodynamic limit. We review several results for the shift of the critical 
temperature and for the temperature dependence of thermodynamic functions, like the condensate fraction, the 
chemical potential and the release energy. We also discuss the behavior of the excitations at finite temperature. 

In Sec. VI we illustrate some features of these trapped Bose gases in connection with superfluidity and phase 
coherence. We discuss in particular the structure of quantized vortices and the behavior of the moment of inertia, as 
well as interference phenomena and quantum effects beyond mean-field theory, like the collapse-revival of collective 
oscillatio ns. 

In Sec. VII we draw our conclusions and we discuss some further future perspectives in the field. 

The overlap between current theoretical and experimental investigations of BEG in trapped alkalis is already wide 
and rich. Various theoretical predictions, concerning the ground state, dynamics and thermodynamics are found 
to agree very well with observations; others are stimulating new experiments. The comparison between theory and 
experiments then represents an exciting feature of these novel systems, which will be frequently emphasized in the 
present review. 



II. THE IDEAL BOSE GAS IN A HARMONIC TRAP 



A. The condensate of noninteracting bosons 

An important feature characterizing the available magnetic traps for alkali atoms is that the confining potential 
can be safely approximated with the quadratic form 



y..,(r) = y(-^-^+-,V+c.,V). (1) 



Thus the investigation of these systems starts as a textbook application of nonrclativistic quantum mechanics for 
identical point-like particles in a harmonic potential. 

The first step consists in neglecting the atom-atom interaction. In this case, almost all predictions are analytical 
and relatively simple. The many-body Hamiltonian is the sum of single-particle Hamiltonians whose eigenvalues have 
the form 

c + hujx + [riy + 2^ Tt-^v + {riz + fiuJz , (2) 

where {nx,ny,nz} are non-negative integers. The ground state 0(ri, ..,vn) of N noninteracting bosons confined by 
the potential (|l|) is obtained by putting all the particles in the lowest single-particle state [ux = Uy = Uz =0), namely 
4>{vi, ..,r7v) = Hi V'o(ri), where ipo{v) is given by 

<^o(r) = (^-^J exp 1^- — (a)a;a; ^ ^yV +i^zZ)^ , (3) 

and we have introduced the geometric average of the oscillator frequencies: 

t^ho = (wi^Wj^w^)^/^ . (4) 

The density distribution then becomes n(r) — A^|(po(r)P and its value grows with N . The size of the cloud is instead 
independent of A'^ and is fixed by the harmonic oscillator length 

Oho = (5) 
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which corresponds to the average width of the Gaussian . This is the first important length scale of the system. 
In the available experiments, it is typically of the order of aho ~ 1 M^- finite temperature only part of the 
atoms occupy the lowest state, the others being thermally distributed in the excited states at higher energy. The 
radius of the thermal cloud is larger than aho- A rough estimate can be obtained by assuming ksT ^ hLUi^o and 
approximating the density of the thermal cloud with a classical Boltzmann distribution nci{r) oc expl—Vcxtif) / ksT]. 
If Ve^tir) = {l/2)mLu^^r^ , the width of the Gaussian is Rt — ahoiksT /tiwhoy/^ , and hence larger than aho- The use 
of a Bose distribution function does not change significantly this estimate. 

The above discussion reveals that Bose-Einstein condensation in harmonic traps shows up with the appearance 
of a sharp peak in the central region of the density distribution. An example is shown in Fig. ^ where we plot the 
prediction for the condensate and thermal densities of 5000 noninteracting particles in a spherical trap at a temperature 
T = 0.9T°, where is the temperature at which condensation occurs (see discussion in the next section). The curves 
correspond to the column density, namely the particle density integrated along one direction, n(z) = J dx n(x,0, z); 
this is a typical measured quantity, the x direction being the direction of the light beam used to image the atomic 
cloud. By plotting directly the density n(r), the ratio of the condensed and noncondensed densities at the center 
would be even larger. 

By taking the Fourier transform of the ground state wave function, one can also calculate the momentum distribution 
of the atoms in the condensate. For the ideal gas, it is given by a Gaussian centered at zero momentum and having 
a width proportional to a^^. The distribution of the thermal cloud is, also in momentum space, broader. Using 
a classical distribution function one finds that the width is proportional to (fcsT)^/^. Actually, the momentum 
distributions of the condensed and noncondensed particles of an ideal gas in harmonic traps have exactly the same 
form as the density distributions no and nx shown in Fig. |^. 

The appearence of the condensate as a narrow peak in both co-ordinate and momentum space is a peculiar feature 
of trapped Bose gases having important consequences in both the experimental and theoretical analysis. This is 
different from the case of a uniform gas where the particles condense into a state of zero momentum, but BEG cannot 
be revealed in co-ordinate space, since the condensed and noncondensed particles fill the same volume. 

Indeed, the condensate has been detected experimentally as the occurrence of a sharp peak over a broader dis- 
tribution, in both the velocity and spatial distributions. In the first case, one lets the condensate expand freely, by 
switching-off the trap, and measures the density of the expanded cloud with light absorption (Anderson et at, 1995). 
If the particles do not interact, the expansion is ballistic and the imaged spatial distribution of the expanding cloud 
can be directly related to the initial momentum distribution. In the second case, one measures directly the density of 
the atoms in the trap by means of dispersive light scattering (Andrews et ai, 1996). In both cases, the appearence of a 
sharp peak is the main signature of Bose-Einstein condensation. An important theoretical task consists of predicting 
how the shape of these peaks is modified by the inclusion of two-body interactions. As anticipated in Fig. ||, the 



interactions can change the picture drastically. This effect will be deeply discussed in Sec. |I]. 

The shape of the confining field fixes also the symmetry of the problem. One can use spherical or axially symmetric 
traps, for instance. The first experiments on rubidium and sodium were carried out with axial symmetry. In this case 
one can define an axial co-ordinate z and a radial co-ordinate r± = (x^ -|-t/^)^/^ and the corresponding frequencies, ujz 
and uj± = uJx = The ratio between the axial and radial frequencies, A = LOzjio^^^ fixes the asymmetry of the trap- 
For A < 1 the trap is cigar-shaped while for A > 1 is disk-shaped. In terms of A the ground state (^ for noninteracting 
bosons can be rewritten as 



Ai/4 

^O(^) ^ 3/2 



(6) 



Here = (?i/ma;j^)^/^ is the harmonic oscillator length in the x-y plane and, since = A ^/^ajho, one has also 
a_L = A^/'^Oho- 

The choice of an axially symmetric trap has proven useful for providing further evidence of Bose-Einstein conden- 
sation from the analysis of the momentum distribution. To understand this point, let us take the Fourier transform of 
the wave function (o): i^o(p) oc exp[— a^(|5^ + A^^p^)/2?i^]. From this one can calculate the average axial and radial 
widths. Their ratio, 

^(p2)/(p2^)=VA, (7) 

is fixed by the asymmetry parameter of the trap. Thus, the shape of the expanded cloud in the x-z plane is an 
ellipse, the ratio between the two axis {aspect ratio) being equal to a/A. If the particles, instead of being in the lowest 
state (condensate), were thermally distributed among many eigenstates at higher energy, their distribution function 
would be isotropic in momentum space, according to the equipartition principle, and the aspect ratio would be equal 
to 1. Indeed, the occurrence of anisotropy in the condensate peak has been interpreted from the very beginning as 
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an important signature of BEC (Anderson et al., 1995; Davis et aL, 1995; Mewes et al, 1996a). In the case of the 
experiment at the Joint Institute for Laboratory Astrophysics (JILA) in Boulder, the trap is disk-shaped with A = 
The first measured value of the aspect ratio was about 50% larger than the prediction, \/A, of the noninteracting 
model (Anderson et at, 1995). Of course, a quantitative comparison can be obtained only including the atom-atom 
interaction, which affects the dynamics of the expansion (Holland and Cooper, 1996; Dalfovo and Stringari, 1996; 
Holland et ai, 1997; Dalfovo et al., 1997c). However, the noninteracting model already points out this interesting 
effect due to anisotropy. 



B. Trapped bosons at finite temperature: tliermodynamic limit 

At temperature T, the total number of particles is given, in the grand-canonical ensemble, by the sum 

N= {exp[/?(e«.n„«. -A^)]-l}~\ (8) 

while the total energy is given by 

^= H ^"xn„n. {exp[/3(£„,„„„, -/i)] - l}"\ (9) 

where is the chemical potential and j3 — {kBT)~^. Below a given temperature the population of the lowest 
state becomes macroscopic and this corresponds to the onset of Bose-Einstein condensation. The calculation of the 
critical temperature, the fraction of particles in the lowest state (condensate fraction) and the other thermodynamic 
quantities, starts from Eqs. (||) and (P) with the appropriate spectrum en^nyu^ (de Groot, Hooman and Ten Seldam, 
1950; Bagnato, Pritchard and Kleppner, 1987). Indeed the statistical mechanics of these trapped gases is less trivial 
than expected at first sight. Several interesting problems arise from the fact that these systems have a finite size and 
are inhomogeneous. For example, the usual definition of thermodynamic limit (increasing N and volume with the 
average density kept constant) is not appropriate for trapped gases. Moreover the traps can be made very anisotropic, 
reaching the limit of quasi-2D and quasi-lD systems, so that interesting effects of reduced dimensionality can be also 
investigated. 

As in the case of a uniform Bose gas, it is convenient to separate out the lowest eigenvalue eooo from the sum (||) 
and call TVq the number of particles in this state. This number can be macroscopic, i.e., of the order of A'^, when the 
chemical potential becomes equal to the energy of the lowest state, 

3 

M ^ /^c = 1 (10) 

where ui = (ujx + tOy + ujz)/^ is the arithmetic average of the trapping frequencies. Inserting this value in the rest of 
the sum, one can write 

^ ^0 ^ eyi-p[(3h{ujxnx+ujyny+u!znz)]-l' ^^^^ 

In order to evaluate this sum explicitly, one usually assumes that the level spacing becomes smaller and smaller when 
N ^ oo, so that the sum can be replaced by an integral: 

N-No= / r^^T " " 1 ■ (12) 

Jq exp[fm{u)xnx + LUyUy + aJzUz)] - 1 

This assumption corresponds to a semiclassical description of the excited states. Its validity implies that the relevant 
excitation energies, contributing to the sum (pT|), are much larger than the level spacing fixed by the oscillator 
frequencies. The accuracy of the semiclassical approximation (p^) is expected to be good if the number of trapped 



atoms is large and fc^T 3> ?iWho- It can be tested a posteriori by comparing the integral (12) with the numerical 
summation (|ll|). 

The integral (^2|) can be easily calculated by changing variables {fihuixUx = fix, etc.). One finds 

N-N, = ai)(^)\ (13) 
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where (^(n) is the Riemann C- function and Who is the geometric average (^). From this result one can also obtain the 
transition temperature for Bose-Einstein condensation. In fact, by imposing that A'o ^ at the transition, one gets 

keT^ = huji^o [y^J = 0-94 ?jwho A^'/' • (14) 

For temperatures higher than the chemical potential is less than fic and becomes A^-dependent, while the population 
of the lowest state is of the order of 1 instead of N. The proper thermodynamic limit for these systems is obtained 
by letting N —> oo and tJho ~^ 0, while keeping the product Nujf^^ constant. With this definition the transition 



temperature (14) is well defined in the thermodynamic limit. Inserting the above expression for into Eq. (fi^) one 



gets the T-dependence of the condensate fraction for T <T^: 

The same result can be also obtained by rewriting (|l^) as an integral over the energy, in the form 

Jo exp(/3e) - 1 

where p{e) is the density of states. The latter can be calculated b y u sing the spectrum (0) and turns out to be 



quadratic in e: p{e) = (l/2)(?ia;ho) e • Inserting this value into < \16\) , one finds again result (|lj). The integral 
E = Jq°° (ie/ci(e)e/[exp(/3e) — 1] gives instead the total energy of the system (^ for which one finds the result 



E _ 3C(4) / T ^ ^ 



NkeTO C(3) \T0 



(17) 



Starting from the energy one can calculate specific heat, entropy and the other thermodynamic quantities. 

These results can be compared with the well known theory of uniform Bose gases (see, for example, Huang, 1987). 
In this case, the eigenstates of the Hamiltonian are plane waves of energy e — /{2m), with the density of states given 
by p(e) = {2n)-^V{2m/h^f/^^, where V is the volume. The sum (|) gives Nq/N = 1 - {T/T^f/"^ and ksT^ = 
(27rn7m)[n/C(3/2)]2/3, with n = N/V, while the energy is given by E/{NkBT°) = 3C(5/2)/[2C(3/2)](T/T0)5/2. 

Another quantity of interest, which can be easily calculated using the semiclassical approximation, is the density 
of thermal particles nrijc). The sum of nT{v) and the condensate density, no(r) = Ao|<Po(r)P, gives the total density 
n(r) = no(r) + ^^(r)- At T < and in the thermodynamic limit, the thermal density is given by the integral over 
momentum space nT{v) — J dp(27r?i)^'^[exp(/3e(p, r)) — 1]^^, where £(p,r) ~ p'^/2m + Vc^t{^) is the semiclassical 
energy in phase space. The result is 

nT(r) = g^^^ (^g-/3Vo..(r)) ^ (18) 

where At — [2iih^ / {mkBT)Y^'^ is the thermal wavelength. The function .g3/2(r) belongs to the class of functions 
ga{z) = Y^^=i [see, for example, Huang (1987)]. By integrating nyW over space one gets again the number 

of thermally depleted atoms N — Nq — N{T/T^)^ , consistently with Eq. ([l5|). In a similar way one can obtain the 
distribution of thermal particles in momentum space: nT(p) = (ATWWho)" y3/2(6xp(— /3p^/2rn)). 

The above analysis points out the existence of two relevant scales of energy for the ideal gas: the transition 
temperature, ksT^, and the average level spacing, Tiwho- From expression (p^), one clearly sees that fcsT^ can be 
much larger than ?iwho- In the available traps, with A^ ranging from a few thousand to several millions, the transition 
temperature is 20 to 200 times larger than TiWho- This also means that the semiclassical approximation is expected 
to work well in these systems on a wide and useful range of temperatures. The frequency a;ho/(27r) is fixed by the 
trapping potential and ranges typically from tens to hundreds of Hertz. This gives ?iWho of the order of a few nK. 
In one of the first experiments at JILA (Ensher et ai, 1996) for example, the average level spacing was about 9 nK, 
corresponding to a critical temperature [see Eq. (p^j of about 300 nK with 40000 atoms in the trap. We also note 
that, for the ideal gas, the chemical potential is of the same order of ?iWho, as shown by Eq. (^0|). However, as we 
will see later on, its value depends significantly on the atom-atom interaction and shall consequently provide a third 
important scale of energy. 

The noninteracting harmonic oscillator model has guided experimentalists to the proper value of the critical tem- 
perature. In fact, the measured transition temperature was found to be very close to the ideal gas value (M|; the 
occupation of the condensate becoming macroscopically large below the critical temperature as predicted by (O) . As 
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an example, in Fig. ^ we show the first experimental results obtained at JILA (Ensher et at, 1996). The occurrence 
of a sudden transition at T/T^ '--^ 1 is evident. Similar results have been obtained also at MIT (Mewes et ai, 1996a). 
Apart from problems related to temperature calibration, a more quantitative comparison between theory and exper- 
iments requires the inclusion of two main effects: the fact that these gases have a finite number of particles and that 
they are interacting. The role of interactions will be analysed extensively in the next sections. Here we briefly discuss 
the relevance of finite size corrections. 



C. Finite size effects 



The number of atoms that can be put into the traps is not truly macroscopic. So far experiments have been carried 
out with a maximum of about lO'^ atoms. As a consequence, the thermodynamic limit is never reached exactly. A 
first effect is the lack of discontinuities in the thermodynamic functions. Hence Bose-Einstein condensation in these 
trapped gases is not, strictly speaking, a phase transition. In practice, however, the macroscopic occupation of the 
lowest state occurs rather abruptly as temperature is lowered and can be observed, as clearly shown in Fig. ^. The 
transition is actually rounded with respect to the predictions of the N —>■ oo limit, but this effect, though interesting, 
is small enough to make the words transition and critical temperature meaningful even for finite-sized systems. It 
is also worth noticing that, instead of being a limitation, the fact that N is finite makes the system potentially 
richer, because new interesting regimes can be explored even in cases where there is no real phase transition in the 



thermodynamic limit. An example is BEC in ID, as we will see in Sec. II D. 

In order to work out the thermodynamics of a noninteracting Bose gas, all one needs is the spectrum of single 
particle levels entering the Bose distribution function. Working in the grand-canonical ensemble for instance, the 
average number of atoms is given by the sum (||) and it is not necessary to take the N ^ oo limit. In fact, the explicit 
summation can be carried out numerically (Ketterle and van Druten, 1996b) for a fixed number of particles and a 
given temperature, the chemical potential being a function of N and T. The condensate fraction No{T)/N, obtained 
in this way, turns out to be smaller than the thermodynamic limit prediction ( |l5| ) and, as expected, the transition is 
rounded off. An example of an exact calculation of the condensate fraction for 1000 noninteracting particles is shown 
in Fig. ^ (circles). With their numerical calculation, Ketterle and van Druten (1996b) found that finite size effects 
are significant only for rather small values of N, less than about 10^. They calculated also the occupation of the first 
excited levels, finding that the fraction of atoms in these states vanishes for iV — > oo and is very small already for N 
of the order of 100. 

The first finite size correction to the law (|l^) for the condensate fraction can be evaluated analytically by studying 
the large N limit of the sum (^) (Grossmann and Holthaus, 1995; Ketterle and van Druten, 1996b; Kirsten and Toms 
1996; Haugerud, Haugset and Ravndal, 1997). The result for NaiT)/N is given by 



N \T^J 2c^ho[C(3)]2/3 \T0 

To the lowest order, finite size effects decrease as iV~^/^ and depend on the ratio of the arithmetic (lj) and geometric 
{ujho) averages of the oscillator frequencies. For axially symmetric traps this ratio depends on the deformation 
parameter A = uJz/uj± as w/who = (A -I- 2)/(3A^/'^). For N = 1000 prediction ( p^ ) is already indistinguishable from 
the exact result obtained by summing explicitly over the excited states of the harmonic oscillator Hamiltonian, apart 
from a narrow region near where higher order corrections should be included to get the exact result. This is well 
illustrated in Fig.^, where we plot the prediction ([l^ ) (solid line) together with the exact calculation obtained directly 
from (||) (circles). Both predictions are also compared with the thermodynamic limit, Nq/N = 1 — [T/T^Y' . 

Finite size effects reduce the condensate fraction and thus result in a lowering of the transition temperature as 
compared to the N ^ oo limit. By setting the left hand side of Eq. (|l^) equal to zero one can estimate the shift of 
the critical temperature to order iV"^/'^ (Grossmann and Holthaus, 1995; Ketterle and van Druten, 1996b; Kirsten 
and Toms, 1996): 

= ?mL/3 ^"^^^ - - 0-73 — N-^/^ . (20) 

Another problem, which deserves to be mentioned in connection with the finite size of the system, is the equivalence 
between different statistical ensembles and the problem of fluctuations. In the thermodynamic limit the grand canon- 
ical, canonical and microcanonical ensembles are expected to provide the same results. However, their equivalence is 
no longer ensured when N is finite. Rigorous results concerning the ideal Bose gas in a box and, in particular, the 
behavior of fluctuations, can be found in Ziff et al. (1977), and Angelescu et al. (1996). In the case of a trapped 
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gas, Gajda and Rz§,zewski (1997) have shown that the differences between the predictions of the micro- and grand 
canonical ensembles for the temperature dependence of the condensate fraction are small already at ~ 1000. The 
fluctuations of the number of atoms in the condensate are instead much more sensitive to the choice of the ensemble 
(Navez et ai, 1997; Wilkens and Weiss, 1997; see also Holthaus, Kalinowski and Kirsten, 1998, and references therein). 
Inclusion of two-body interactions can, however, change the scenario significantly (Giorgini, Pitaevskii and Stringari, 
1998). 



D. Role of dimensionality 

So far we have discussed the properties of the ideal Bose gas in three-dimensional space. Though the trapping 
frequencies in each direction can be quite different, nevertheless the relevant results for the temperature dependence 
of the condensate have been obtained assuming that fc^T is much larger than all the oscillator energies tiLOx, hiOy, Tiujz- 
In order to observe effects of reduced dimensionality, one should remove such a condition in one or two directions. 

The statistical behavior of 2D and ID Bose gases exhibits very peculiar features. Let us first recall that in a 
uniform gas Bose-Einstein condensation cannot occur in 2D and ID at finite temperature because thermal fluctuations 
destabilize the condensate. This can be seen by noting that, for an ideal gas in the presence of BEG, the chemical 
potential vanishes and the momentum distribution, n{p) cx [exp(/3p^/2m) — 1]""'^, exhibits an infrared l/p^ divergence. 
In the thermodynamic limit, this yields a divergent contribution to the integral /dp n{p) in 2D and ID, thereby 
violating the normalization condition. The absence of BEC in ID and 2D can be also proven for interacting uniform 
systems, as shown by Hohenberg (1967). 

In the presence of harmonic trapping, the effects of thermal fluctuations are strongly quenched due to the different 
behavior exhibited by the density of states p{e). In fact, while in the uniform gas p{e) behaves as e('^^2)/2^ where 
d is the dimensionality of space, in the presence of an harmonic potential one has instead the law p{e) ~ e'^^^ and, 
consequently, the integral (16) converges also in 2D. The corresponding value of the critical temperature is given by 

N 



kBT2D = huJ2D ' (21) 

where uj2d = {^x'^yY^'^ (see, for example, MuUin, 1997, and references therein). One notes first that in 2D the 
thermodynamic limit corresponds to taking N ^ oo and uj2d — > with the product Nuj2jj kept constant. In order 
to achieve 2D Bose-Einstein condensation in real 3D traps, one should choose the frequency ujz in the third direction 
large enough to satisfy the condition hoj2D ^ kBT2D < ^^z', this implies rather severe conditions on the deformation 
of the trap. The main features of BEG in 2D gases confined in harmonic traps and, in particular, the applicability of 
the Hohenberg theorem and of its extensions to nonuniform gases, have been discussed in details by MuUin (1997). 

In ID the situation is also very interesting. In this case, Bose-Einstein condensation cannot occur even in the 
presence of harmonic confinement because of the logarithmic divergence in the integral (|l^). This means that the 
critical temperature for ID Bose-Einstein condensation tends to zero in the thermodynamic limit if one keeps the 
product NuoiD fixed. In fact, in ID the critical temperature for the ideal Bose gas can be estimated to be (Ketterle 
and van Drutcn, 1996b) 

N 

ksT,n = nu:,n^^^ (22) 

with LOiB = LJz- Despite the fact that one cannot have BEG in the thermodynamic limit, nevertheless for finite values 
of N the system can exhibit a large occupation of the lowest single-particle state in a useful interval of temperatures. 
Furthermore, if the value of N and the parameters of the trap are chosen in a proper way, one observes a new interesting 
phenomenon associated with the macroscopic occupation of the lowest energy state, taking place in two distinct steps 
(van Druten and Ketterle, 1997). This happens when the relevant parameters of the trap satisfy simultaneously the 
conditions Tid < T^u and huj± < kBT^o, where T^jj coincides with the usual critical temperature given in Eq. ( p^ ) 
and Luj_ is the frequency of the trap in the x-y plane. In the interval Tid < T < T^d, only the radial degrees of 
freedom are frozen, while no condensation occurs in the axial degrees of freedom. At lower temperatures, below Tio, 
also the axial variables start being frozen and the overall ground state is occupied in a macroscopic way. An example 
of this two-step BEG is shown in Fig. |7[ It is also interesting to notice that the conditions for the occurrence of 
two-step condensation in harmonic potentials are peculiar of the ID geometry. In fact, it is easy to check that the 
corresponding conditions T2D < T^d and huz < ksT^o, which would yield two-step BEC in 2D, cannot be easily 
satisfied because of the absence of the IniV factor. 
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It is finally worth pointing out that the above discussion concerns the behavior of the ideal Bose gas. Effects of 
two-body interactions are expected to modify in a deep way the nature of the phase transition in reduced dimension- 
ality. In particular, interacting Bose systems exhibit the well known Berezinsky-Kosterlitz-Thouless transition in 2D 
(Berezinsky, 1971; Kosterlitz and Thouless, 1973). The case of trapped gases in 2D has been recently discussed by 
MuUin (1998) and is expected to become an important issue in future investigations. 



E. Non harmonic traps and adiabatic transformations 



A crucial step to reach the low temperatures needed for BEC in the experiments realized so far is evaporative cooling. 
This technique is intrinsically irreversible since it is based on the loss of hot particles from the trap. New interesting 
perspectives would open if one could adiabatically cool the system in a reversible way (Ketterle and Pritchard, 1992; 
Pinkse et ai, 1997). Reversible cooling of the gas is achieved by adiabatically changing the shape of the trap at a 
rate slow compared to the internal equilibration rate. 

An important class of trapping potentials for studying the effects of adiabatic changes is provided by power-law 
potentials of the form 

l/ext(r) = Ar" , (23) 
where, for simplicity, we assume spherical symmetry. The critical temperature for Bose-Einstein condensation in the 



trap (23) has been calculated by Bagnato, Pritchard and Kleppner (1987) and is given by 



(2m)3/2 r(l + ,5)C(3/2 + (5) 



(24) 



Here we have introduced the parameter S = 3/a, while r(a;) is the usual gamma function. By setting 5 = 3/2 and 
A = one recovers the result for the transition temperature in an isotropic harmonic trap. The result for a 

rigid box is instead obtained by letting 5 — > 0. 

It is straightforward to work out the thermodynamics of a noninteracting gas in the confining potential ( [2^ ) 
(Bagnato, Pritchard and Kleppner, 1987; Pinkse et al. 1997). For example, for the condensate fraction one finds: 
Nq/N = 1~{T/T^)'^^^^^ . More relevant to the discussion of reversible processes is the entropy which remains constant 
during the adiabatic change. Above Tc the system can be approximated by a classical Maxwell-Boltzmann gas and 
the entropy per particle takes the simple form 

^ +5-\nC{3/2 + S)) + (l + s)\n(^) . (25) 



NkB \2 ' 7 V2 / \T^, 

From this equation one sees that the entropy depends on the parameter A of the external potential (|2^) only through 
the ratio T/T^. Thus, for a fixed power-law dependence of the trapping potential {5 fixed), an adiabatic change of 
A, like for example an adiabatic expansion of the harmonic trap, does not bring us closer to the transition, since the 
ratio T/r" remains constant. A reduction of the ratio T /T^ is instead obtained by increasing adiabatically (5, that 
is, changing the power-law dependence of the trapping potential (Pinkse et al. 1997). For example, in going from a 

harmonic (i5i = 3/2) to a linear trap {82 = 3), one gets the relation t2 — Q.7t\^^ between the initial and final reduced 
temperature t = T/T^. In this case a system at twice the critical temperature (ti = 2) can be cooled down to nearly 
the critical point {t2 — 1.1). Using this technique it should be possible, by a proper change of S, to cool adiabatically 
the system from the high temperature phase without condensate down to temperatures below Tc with a large fraction 
of atoms in the condensate state. 

The possibility of reaching BEC using adiabatic transformations has been recently successfully explored in an 
experiment carried out at MIT (Stamper-Kurn et ai, 1998b). 



III. EFFECTS OF INTERACTIONS: GROUND STATE 



A. Order parameter and mean-field theory 



The many body Haniiltonian describing TV interacting bosons confined by an external potential Vext is given, in 
second quantization, by: 
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dr ¥(r) + Foxt(r) *(r) + - / drdr' ¥(rWUr')V(r - r')*(r')*(r) (26) 

2m J 2 J 

where ^'(r) and ^'''(r) are the boson field operators that annihilate and create a particle at the position r, respectively, 
and V{r — r') is the two-body interatomic potential. 

The ground state of the system as well as its thermodynamic properties can be directly calculated starting from 
the Hamiltonian (26). For instance, Krauth (1996) has used a Path Integral Monte Carlo method to calculate the 
thermodynamic behavior of 10** atoms interacting with a repulsive "hard-sphere" potential. In principle, this procedure 
gives exact results within statistical errors. However, the calculation can be heavy or even impracticable for systems 
with much larger values of N. Mean-field approaches are commonly developed for interacting systems in order to 
overcome the problem of solving exactly the full many-body Schrodinger equation. Apart from the convenience of 
avoiding heavy numerical work, mean-field theories allow one to understand the behavior of a system in terms of a 
set of parameters having a clear physical meaning. This is particularly true in the case of trapped bosons. Actually 
most of the results reviewed in this paper show that the mean-field approach is very effective in providing quantitative 
predictions for the static, dynamic and thermodynamic properties of these trapped gases. 

The basic idea for a mean-field description of a dilute Bose gas was formulated by Bogoliubov (1947). The key point 
consists in separating out the condensate contribution to the bosonic field operator. In general, the field operator 
can be written as ^'(r) = 5'a(r)aa, where ^'a(r) are single-particle wave functions and Aq are the corresponding 
annihilation operators. The bosonic creation and annihilation operators aj, and Oq, are defined in Fock space through 
the relations 

a'l\no,ni,...,na,...)= + 1 I "o, ni, . . . , -|- 1, . . .) (27) 
Qa I no,ni, . . . ,nQ, . . .) = y/n^ | no, ni, . . . , - 1, . . .) (28) 

where Ua are the eigenvalues of the operator fia = o-aO'a giving the number of atoms in the single-particle a-state. 
They obey the usual commutation rules: 

[aa,a^p] = Sa.p , [aa,a^]=0, [ajj,a^]=0. (29) 

Bose-Einstein condensation occurs when the number of atoms uq of a particular single-particle state becomes very 
large: uq = Nq ^ 1 and the ratio Nq/N remains finite in the thermodynamic limit N — > oo. In this limit the states 
with Nq and A^o ± 1 — A'o correspond to the same physical configuration and, consequently, the operators aq and 
aj can be treated like c-numbers: ao — aj^ — \/Nq. For a uniform gas in a volume V, where BEC occurs in the 



single-particle state 5*0 = 1/v V having zero momentum, this means that the field operator ^'(r) can be decomposed 
in the form \['(r) = \/Nq/V + ^'(r). By treating the operator ^' as a small perturbation, Bogoliubov developed the 
"first-order" theory for the excitations of interacting Bose gases. 

The generalization of the Bogoliubov prescription to the case of nonuniform and time dependent configurations is 
given by 

*(r,i) = $(r,i) -|-*'(r,i) , (30) 

where we have used the Heisenberg representation for the field operators. Here <i>(r,f) is a complex function defined 
as the expectation value of the field operator: $(r, i) = (>I'(r,t)). Its modulus fixes the condensate density through 
no(r, i) = |$(r,f)p. The function ^{r^t) possesses also a well-defined phase and, similarly to the case of uniform 
gases, this corresponds to assuming the occurrence of a broken gauge symmetry in the many-body system. 

The function <I>(r, t) is a classical field having the meaning of an order parameter and is often called "wave function of 
the condensate". It characterizes the off-diagonal long-range behavior of the one-particle density matrix pi(r',r,t) = 
{'^'^{v' ,t)'^{T,t)). In fact the decomposition ( |30| ) implies the following asymptotic behavior (Ginzburg and Landau, 
1950; Penrose, 1951; Penrose and Onsager, 1956): 

lim pi(r',r,t) -$*(r',i)$(r,t) . (31) 

r' — r| — >oc 

Notice that, strictly speaking, in a finite-sized system neither the concept of broken gauge symmetry, nor the one of 
off-diagonal long-range order can be applied. The condensate wave function $ has nevertheless still a clear meaning: it 
can be in fact determined through the diagonalization of the one-body density matrix, J dr'pi{r' , r)^i{r') — Ni^i{r), 
and corresponds to the eigenfunction, <I>i, with the largest eigenvalue, A^^. This procedure has been used, for example, 
to explore Bose-Einstein condensation in finite drops of liquid helium by Lewart et al. (1988). The connection 
between the condensate wave function, defined through the diagonalization of the density matrix and the concept 
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of order parameter commonly used in the theory of superfluidity, is an interesting and nontrivial problem in itself. 
Another important question concerns the possible fragmentation of the condensate, taking place when two or more 
eigenstates of the density matrix /5i(r',r) are macroscopically occupied. One can show (Nozieres and Saint James, 
1982; Nozieres, 1995) that, due to exchange effects, in uniform gases interacting with repulsive forces the fragmentation 
costs a macroscopic energy. The behavior can be however different in the presence of attractive forces and almost 
degenerate single-particle states (Nozieres and Saint James, 1982; Kagan, Shlyapnikov and Walraven, 1996; Wilkin, 
Gunn and Smith, 1998). 

The decomposition ( ^0| ) becomes particularly useful if ^' is small, i.e., when the depiction of the condensate is 
small. Then, an equation for the order parameter can be derived by expanding the theory to the lowest orders in ^f' 
as in the case of uniform gases. The main difference is that here one gets also a nontrivial "zeroth-order" theory for 

In order to derive the equation for the condensate wave function <l>(r, t), one has to write the time evolution of the 
field operator 5'(r, t) using the Heisenberg equation with the many-body Hamiltonian (pq): 



2m 



+ V;^t(r)+ / dr' ii\r',t)V{r' -r)i'{r',t) 



*(r,t) 



(32) 



Then one has to replace the operator with the classical field In the integral containing the atom- atom interaction 
V{r' — r), this replacement is, in general, a poor approximation when short distances (r' — r) are involved. In a dilute 
and cold gas, one can nevertheless obtain a proper expression for the interaction term by observing that, in this case, 
only binary collisions at low energy are relevant and these collisions are characterized by a single parameter, the 
s-wave scattering length, independently of the details of the two-body potential. This allows one to replace V{r' — r) 
in (32) with an effective interaction 



F(r' - r) = gS{r' - r) 



where the coupling constant g is related to the scattering length a through 
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(33) 



(34) 



The use of the effective potential (|3^) in 
closed equation for the order parameter: 



p2) is compatible with the replacement of ^ with $ and yields the following 



zr.|$(r,t) = 



2m 



+ Kxt(r)+.g|$(r,i)|^ <f(r,t) 



(35) 



This equation, known as Gross-Pitaevskii (CP) equation, was derived independently by Gross (1961 and 1963) and 
Pitaevskii (1961). Its validity is based on the condition that the s-wave scattering length be much smaller than 
the average distance between atoms and that the number of atoms in the condensate be much larger than 1. The 
GP equation can be used, at low temperature, to explore the macroscopic behavior of the system, characterized by 
variations of the order parameter over distances larger than the mean distance between atoms. 
The Gross-Pitaevskii equation (Bq) can also be obtained using a variational procedure: 



where the energy functional E is given by 



2m' 



iv$|- 



6E 



Kxt(r)|$| 



(36) 



(37) 



The first term in the integral (|3^) is the kinetic energy of the condensate E'kin, the second is the harmonic oscillator 
energy i?ho, while the last one is the mean- field interaction energy i^int. Notice that the mean- field term, -Bint, 
corresponds to the first correction in the virial expansion for the energy of the gas. In the case of non-negative and 
finite-range interatomic potentials, rigorous bounds for this term have been obtained by Dyson (1967) and Lieb and 
Yngvason (1998). 
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The dimensionless parameter controlling the validity of the dilute-gas approximation, required for the derivation of 
Eq. (^5|), is the number of particles in a "scattering volume" |a|'^. This can be written as n|a|^, where n is the average 
density of the gas. Recent determinations of the scattering length for the atomic species used in the experiments on 
BEC give: a = 2.75 nm for ^SNa (Tiesinga et at, 1996), a = 5.77 nm for ^^Rb (Boesten et ai, 1997) and a = -1.45 
nm for ^Li (Abraham et ai, 1995). Typical values of density range instead from 10^'^ to 10^^ cm~'^, so that n|ap is 
always less than 10^"^. 

When fl|ap ^ 1 the system is said to be dilute or weakly interacting. However, one should better clarify the 
meaning of the words "weakly interacting" , since the smallness of the parameter njap does not imply necessarily that 
the interaction effects are small. These effects, in fact, have to be compared with the kinetic energy of the atoms in the 
trap. A first estimate can be obtained by calculating the interaction energy, E'int, on the ground state of the harmonic 
oscillator. This energy is given by gNfi, where the average density is of the order of N/a\^, so that E^int oc N'^\a\/a\^. 
On the other hand, the kinetic energy is of the order of NTlloi^q and thus E'kin Na^^ . One finally finds 

^oc^. (38) 

E^n Oho 

This is the parameter expressing the importance of the atom-atom interaction compared to the kinetic energy. It can 
be easily larger than 1 even if n\a\^ <C 1, so that also very dilute gases can exhibit an important nonideal behavior, 
as we will discuss in the following sections. In the first experiments with rubidium atoms at JILA (Anderson et ai, 
1995) the ratio |a|/aho was about 7 x 10"'^, with TV of the order of a few thousands. Thus iVa/oho is larger than 
1. In the experiments with ^Li at Rice University (Bradley et ai, 1997; Sackett et ai, 1997) the same parameter 
is smaller than 1, since the number of particles is of the order of 1000 and |a|/aho ~ 0.5 x 10^'^. Finally, in the 
experiments with sodium at MIT (Davis et ai, 1995) the number of atoms in the condensate is very large (10^ - lO'^) 
and iV|a|/aho ~ 10^ - lO". 

Due to the assumption = 0, the above formalism is strictly valid only in the limit of zero temperature, when 
all the particles are in the condensate. The dynamic behavior and the generalization to finite temperatures will be 
discussed in Sees. [V and 0, respectively. Here we present the results for the stationary solution of the Gross-Pitaevskii 



(GP) equation at zero temperature. 



B. Ground state 



For a system of noninteracting bosons in a harmonic trap, the condensate has the form of a Gaussian of average 
width aho [see Eq. (^], and the central density is proportional to N. If the atoms are interacting, the shape of the 
condensate can change significantly with respect to the Gaussian. The scattering length entering the Gross-Pitaevskii 
equation can be positive or negative, its sign and magnitude depending crucially on the details of the atom-atom 
potential. Positive and negative values of a correspond to an effective repulsion and attraction between the atoms, 
respectively. The change can be dramatic when the interaction energy is much greater than the kinetic energy, that 
is, when N\a\/aho ^ 1- The central density is lowered (raised) by a repulsive (attractive) interaction and the radius 
of the atomic cloud consequently increases (decreases). This effect of the interaction has important consequences, not 
only for the structure of the ground state, but also for the dynamics and thermodynamics of the system, as we will 
see later on. 

The ground state can be easily obtained within the formalism of mean-field theory. For this, one can write the 
condensate wave function as $(r, t) = (j){r) exp^—ifit/h), where /i is the chemical potential and (f> is real and normalized 
to the total number of particles, j dv cjp = = N . Then the Gross-Pitaevskii equation (|3^) becomes 

("^ + Kxt(r) + g02(r)^ 0(r) ^^^(r) . (39) 

This has the form of a "nonlinear Schrodinger equation" , the nonlinearity coming from the mean-field term, propor- 
tional to the particle density n(r) = (/)^(r). In the absence of interactions (g = 0), this equation reduces to the usual 
Schrodinger equation for the single-particle Hamiltonian — ?i^/(2m)V^ -I- T4xt(r) and, for harmonic confinement, the 
ground state solution coincides, apart from a normalization factor, with the Gaussian function (^: 4>{y) = ^/Nipo(r). 
We note, in passing, that a similar nonlinear equation for the order parameter has been also considered in connection 
with the theory of superfluid helium near the A-point (Ginzburg and Pitaevskii, 1958); in that case, however, the 
ingredients of the equation have a different physical meaning. 

The numerical solution of the GP equation ( ^9| ) is relatively easy to obtain (Edwards and Burnett, 1995; Ruprecht 
et ai, 1995; Edwards et ai, 1996b; Dalfovo and Stringari, 1996; Holland and Cooper, 1996). Typical wave functions 
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(j)^ calculated from Eq. ( |39| ) with different values of the parameter A^|a|/aho, are shown in Figs. || and || for attractive 
and repulsive interaction, respectively. The effects of the interaction are revealed by the deviations from the Gaussian 
profile (H) predicted by the noninteracting model. Excellent agreement has been found by comparing the solution of 
the GP equation with the experimental density profiles obtained at low temperature (Hau et at, 1998), as shown in 
Fig. ^. The condensate wave function obtained with the stationary GP equation has been also compared with the 
results of an ab initio Monte Carlo simulation starting from Hamiltonian (|26|), finding a very good agreement (Krauth, 
1996). 

The role of the parameter A^|a|/aho, already discussed in the previous section, can be easily pointed out, in the 
Gross-Pitaevskii equation, by using rescaled dimensionless variables. Let us consider a spherical trap with frequency 
ujho and use Oho, OjTo and hujho as units of length, density and energy, respectively. By putting a tilde over the rescaled 
quantities, Eq. (p9[) becomes 



f2 



8^(iVa/aho)0'(f) 0(f) =2A0(f) . 



(40) 



In these new units the order parameter satisfies the normalization condition J dr\(f)\'^ = 1. It is now evident that the 
importance of the atom- atom interaction is completely fixed by the parameter Na/a\^o. 

It is worth noticing that the solution of the stationary GP equation (^9|) minimizes the energy functional (^^ for 
a fixed number of particles. Since the ground state has no currents, the energy is a functional of the density only, 
which can be written in the form 



E\n] 



dr 



2m 



|VV^|2 + nK=xt(r) 



2 



E^ 



ho 



E, 



(41) 



The first term corresponds to the quantum kinetic energy coming from the uncertainty principle; it is usually named 
"quantum pressure" and vanishes for uniform systems. In general, for a nonstationary order parameter, the kinetic 
energy in (|3^) includes also the contribution of currents in the form of an additional term containing the gradient of 
the phase of $. 



By direct integration of the GP equation (39) one finds the useful expression 

^i^iEMn + El,o + 2Eint)/N 



(42) 



for the chemical potential in terms of the different contributions to the energy functional (|4l| ). Further important 
relationships can be also found by means of the virial theorem. In fact, since the energy ( pT^ is stationary for any 
variation of (f) around the exact solution of the GP equation, one can choose scaling transformations of the form 
4>{x, y, z) —^ {1 + t^)^/^0((l + iy)x, y, z), and insert them in (|37|). By imposing the energy variation to vanish at first 
order in v, one finally gets 



{Eyin)x — {Eho)x + 2^int = ^ 



(43) 



where (iJkin)^ = (EjPL)/2w and {Ey,o)x = {m/2)ujl{J2, 



Analogous expressions are found by choosing similar 



scaling transformations for the y and z co-ordinates. By summing over the three directions one finally finds the virial 
relation: 



2-Ekin — 2£'ho + 3£'ii 







(44) 



The above results are exact within Gross-Pitaevskii theory and can be used, for instance, to check the numerical 
solutions of Eq. ( ^ ) . 

In a series of experiments the gas has been imaged after a sudden switching-off of the trap and the kinetic energy 
of the atoms has been measured by integrating over the observed velocity distribution. This energy, which is also 
called release energy, coincides with the sum of the kinetic and interaction energies of the atoms at the beginning of 
the expansion: 



Erc\ — E] 



kin 



E. 



(45) 



During the first phase of the expansion both the quantum kinetic energy (quantum pressure) and the interaction 
energy are rapidly converted into kinetic energy of motion. Then the atoms expand at constant velocity. Since 
energy is conserved during the expansion, its initial value (45), calculated with the stationary GP equation, can 
be directly compared with experiments. This comparison provides clean evidences for the crucial role played by 
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two-body interactions. In fact, the noninteracting model predicts a release energy per particle given by £'rei/-/V = 
(1/2)(1 + \/2)fiL0\io, independent of N. Conversely, the observed release energy per particle depends rather strongly 
on N , in good agreement with the theoretical predictions for the interacting gas. In Figs. |l^ and |ll|, we show the 
experimental data obtained at JILA (Holland et aL,1997) and MIT (Mewes et ai, 1996a), respectively. 

Finally, we notice that the balance between the quantum pressure and the interaction energy of the condensate 
fixes a typical length scale, called the healing length, ^. This is the minimum distance over which the order parameter 
can heal. If the condensate density grows from to n within a distance ^, the two terms in Eq. ( ^9[ ) coming from the 
quantum pressure and the interaction energy are ~ ?i^/(2m^^) and ~ 4:Trh^an/m, respectively. By equating them, 
one finds the following expression for the healing length: 

C = (8^na)-i/2 . (46) 

This is a well known result for weakly interacting Bose gases. In the case of trapped bosons, one can use the central 
density, or the average density, to get an order of magnitude of the healing length. This quantity is relevant for 
superfluid effects. For instance, it provides the typical size of the core of quantized vortices (Gross, 1961; Pitaevskii, 
1961). Note that in condensed matter physics the same quantity is often named "coherence length", but the name 
"healing length" is preferable here in order to avoid confusion with different definitions of coherence length used in 
atomic physics and optics. 



C. Collapse for attractive forces 



If forces are attractive (a < 0), the gas tends to increase its density in the center of the trap in order to lower the 
interaction energy, as seen in Fig. |^. This tendency is contrasted by the zero point kinetic energy which can stabilise 
the system. However, if the central density grows too much, the kinetic energy is no longer able to avoid the collapse 
of the gas. For a given atomic species in a given trap, the collapse is expected to occur when the number of particles 
in the condensate exceeds a critical value A^cr, of the order of aho/|a|. It is worth stressing that in a uniform gas, 
where quantum pressure is absent, the condensate is always unstable. 

The critical number A^cr can be calculated at zero temperature by means of the Gross-Pitaevskii equation. The 
condensates shown in Fig. || are metastable, corresponding to local minima of the energy functional ( |37| ) for different 
A^. When A^ increases, the depth of the local minimum decreases. Above A^cr the minimum no longer exists and the 
Gross-Pitaesvkii equation has no solution. For a spherical trap this happens at (Ruprccht et at, 1995) 

— ^ = 0.575 . (47) 

For the axially symmetric trap with ^Li used in the experiments at Rice University (Bradley et al., 1995 and 1997; 
Sackett et al., 1997), the GP equation predicts A^cr — 1400 (Dalfovo and Stringari, 1996; Dodd et al., 1996); this value 
is consistent with recent experimental measurements (Bradley et al., 1997; Sackett et al., 1997). The same problem 
has been investigated theoretically by several authors (Kagan, Shlyapnikov and Walraven, 1996; Houbiers and Stoof, 
1996; Shuryak, 1996; Pitaevskii, 1996; Bergeman,1997). 

A direct insight into the behavior of the gas with attractive forces can be obtained by means of a variational 
approach based on Gaussian functions (Baym and Pethick, 1996). For a spherical trap one can minimize the energy 
( |37| ) using the ansatz 

N 



where w is a dimensionless variational parameter which fixes the width of the condensate. One gets 

#^ = + - (2.r)-V2^ ^-3 . (49) 

Nnujho 4 flho 

This energy is plotted in Fig. |l^ as a function of w, for several values of the parameter A"|a|/aho. One clearly 
sees that the local minimum disappears when this parameter exceeds a critical value. This can be calculated by 
requiring that the first and second derivative of E{w) vanish at the critical point (w = Wcr and A^ = A'cr). One finds 
Wci- = 5^^/^ « 0.669 and A'crlal/oho ~ 0.671. The last formula provides an estimate of the critical number of atoms, 
for given trap and atomic species, reasonably close to the value ( |47| ) obtained by solving exactly the GP equation. 
The Gaussian ansatz has been used by several authors in order to explore both static and dynamic properties of 
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the trapped gases. The stability of a gas with a < has been explored in details, for instance, by Stoof (1997), 
Perez-Garci'a et al. (1997), Shi and Zheng (1997a), Parola, Salasnich and Reatto (1998). The variational function 
proposed by Fetter (1997), which interpolates smoothly between the ideal gas and the Thomas-Fermi limit for positive 
a, also reduces to a Gaussian for a < 0. 

The behavior of the gas close to collapse could be significantly affected by mechanisms not included in the Gross- 
Pitaevskii theory. Among them, inelastic two- and three-body collisions can cause a loss of atoms from the condensate 
through, for instance, spin exchange or recombination (Hijmans et al, 1993; Edwards et al., 1996b; Moerdijk et at, 
1996; Fedichev et at, 1996). This is an important problem not only for attractive forces but also for repulsive forces 
when the density of the system becomes large. 

Recent discussions about the collapse, including quantum tunneling phenomena, can be found, for instance, in 
Sackett, Stoof and Hulet (1998), Kagan, Muryshev and Shlyapnikov (1998), Ueda and Leggett (1998), Ueda and 
Huang (1998). 



D. Large limit for repulsive forces 

In the case of atoms with repulsive interaction (a > 0), the limit Na/a\^o S> 1 is particularly interesting, since this 
condition is well satisfied by the parameters a and Oho used in most of current experiments. Moreover, in this 
limit the predictions of mean-field theory take a rather simple analytic form (Edwards and Burnett 1995; Baym and 
Pethick 1996). 

As regards the ground state, the effect of increasing the parameter Na/a^o is clearly seen in Fig. ||: the atoms 
are pushed outwards, the central density becomes rather flat and the radius grows. As a consequence, the quantum 
pressure term in the Gross-Pitaevskii equation ( |3^ ) , proportional to y'n(r), takes a significant contribution only 
near the boundary and becomes less and less important with respect to the interaction energy. If one neglects 
completely the quantum pressure in (^9|), one gets the density profile in the form 

n(r) = 02(r) - g-'[^i - 14xt(r)] (50) 

in the region where fi > T4xt(r), and n = outside. This is often referred to as Thomas- Fermi (TF) approximation. 
The normalization condition on n(r) provides the relation between chemical potential and number of particles: 



^^^^ 

2 V flho / 

Note that the chemical potential depends on the trapping frequencies, entering the potential Vcxt given in (|l|), only 
through the geometric average Who [see Eq. (^]. Moreover, since /i = dE/dN, the energy per particle turns out to 
be E/N — (5/7)/i. This energy is the sum of the interaction and oscillator energies, since the kinetic energy gives a 
vanishing contribution for large N . Finally, in the same limit, the release energy (|4^) coincides with the interaction 
energy: E,^i/N = (2/7)^. 

The chemical potential, as well as the interaction and oscillator energies obtained by solving numerically the GP 
equation (^9|) become closer and closer to the Thomas-Fermi values when N increases (see for instance, Dalfovo and 
Stringari, 1996). For sodium atoms in the MIT traps, where N is larger than 10^, the Thomas-Fermi approximation 
is practically indistinguishable from the solution of the GP equation. The release energy per particle measured by 
Mewes et al. (1996a) is indeed well fitted with a 7V^/^ law, as shown in Fig. |Tl|. The same agreement is expected to 
occur for rubidium atoms in the most recent JILA traps, having N larger than 10^ (Matthews et al., 1998). 



The density profile (50) has the form of an inverted parabola, which vanishes at the classical turnin g p oint R 

2 } 

ho-* 



defined by the condition fj, = VcxtCR)- For a spherical trap, this implies /i = mujf^^R'^ /2 and, using result (fell) for ^, 
one finds the following expression for the radius of the condensate 

1/5 



^ = flho 

V aho / 



(52) 



which grows with N . For an axially symmetric trap, the widths in the radial and axial directions are fixed by the 
conditions ^ = = tow^Z^/2. It is worth mentioning that, in the case of the cigar-shaped trap used at 

MIT, with a condensate of about 10'' sodium atoms, the axial width becomes macroscopically large {Z 0.3 mm), 
allowing for direct in situ measurements. 

The value of the density ( pO| ) in the center of the trap is nTF(O) = n/g- It is worth stressing that this density is 
much lower than the one predicted for noninteracting particles. In the latter case, using Eq. (^) one gets nho(O) = 
A^/(7r^/^a^Q). The ratio between the central densities in the two cases is then 
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?^ho(0) 



-3/5 



a\\oJ 



(53) 



and decreases with N . For the available traps with ^'^Na and ^^Rb, where Na/aho ranges from about 10 to 10^, the 
atom-atom repulsion reduces the density by one or two orders of magnitude, which is a quite remarkable effect for 
such a dilute systems. An example was already shown in Fig. ||; in that case, the number of particles is about 80000 
and Na/uho ^ 300. 

In Fig. ^|a we show the density profile for a gas in a spherical trap with Na/ay^o — 100. The comparison with the 
exact solution of the GP equation ( p9| ) shows that the TF approximation is very accurate except in the surface region 
close to R. In part b of the same figure, we plot the column density, n(z) — J dx n{x,0,z), which is the measured 
quantity when the atomic cloud is imaged by light absorption or dispersive light scattering. Using the TF density ( [50| ) 
with Vext = {l/2)mujl^r'^, one finds n{z) = (4/3)[2/(ma;^Q)]^/^5'~^[/^ - {l/2)mu!l^z'^]^^'^ . One notes that the accuracy 
of the Thomas-Fermi approximation is even better in the case of the column density, because the extra integration 
makes the cusp in the outer part of the condensate smoother. 

The only region where the Thomas- Fermi density ( |50| ) is inadequate is close to the classical turning point. This 
region plays a crucial role for the calculation of the kinetic energy of the condensate. The shape of the outer part of 
the condensate is fixed by the balance of the zero point kinetic energy and the external potential. In particular, this 
balance can be used to define an effective surface thickness, d. For a spherical trap, for instance, one can assume the 
two energies to have the form h^/{2md^) and muj^^Rd, respectively. One then gets (Baym and Pethick, 1996) 
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(54) 



this ratio is small when TF approximation is valid, i.e., when R ^ Cho- It is interesting to compare the surface 
thickness d with the healing length (|4^). In terms of the ratio Uho/R one can write ^/i? = (aho/R)^, showing that 
the healing length decreases with TV more rapidly than the surface thickness d. 

A good approximation for the density in the region close to the classical turning point, can be obtained by a suitable 
expansion of the GP equation (p9h. In fact, when |r — i?| ^ R, the trapping potential Vcxtir) can be replaced with 



a linear ramp, mLij^^R{r — R), and the GP equation takes a universal form (Dalfovo, Pitaevskii and Stringari, 1996; 
Lundh, Pethick and Smith, 1997), yielding the rounding of the surface profile. 

Using the above procedure it is possible to calculate the kinetic energy which, in the case of a spherical trap, is 
found to follow the asymptotic law 
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(55) 



where C ~ 1.3 is a numerical factor. Analogous expansions can be derived for the harmonic potential energy, E^o, 
and interaction energy, i?int, in the same large N limit (Fetter and Feder, 1997). A straightforward derivation is 
obtained by using nontrivial relationships among the various energy components -Ekim -E'ho and Ei^t of Eq. (41). A 
first relation is given by the virial theorem (^). A second one is obtained by using expression (^|) for the chemical 
potential and the thermodynamic definition /i — dE/dN. These two relationships, together with the asymptotic 
law (^) for the kinetic energy, allow one to obtain the expansions E^,^/N = (3/7)AiTF + n^limR^) ln[i?/(Caho)] and 
E\nt/N = (2/7)yLtTF — ?i^/(mi?2) ln[_R/(Caho)]- From them one gets the useful results 



/.t = ^TF 



R^ \Ca^o 



(56) 



and 



E ■ 
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R 



W«ho 



(57) 



for the chemical potential and the total energy, respectively. In these equations /ztf and R are the Thomas-Fermi 
values (51) and ( p2| ) of the chemical potential and the radius of the condensate. Equations (55)-(|^, which apply to 
spherical traps, clearly show that the relevant small parameter in the large N expansion is ai^o/R = {\hNa/a\^o)~^/^ . 

The Thomas-Fermi approximation ( |50| ) for the ground state density of trapped Bose gases is very useful not only 
for determining the static properties of the system, but also for dynamics and thermodynamics, as we will see in 
Sees. IV and ^ It is worth noticing that this approximation can be derived more directly using local density theory 
as we are going to discuss in the next section. 
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E. Beyond mean-field theory 



Before cfosing this discussion about the effect of interactions on the ground state properties, we wish to come back 
to the basic question of the vaUdity of the Gross-Pitaevskii theory. All the results so far presented are expected to 
be valid if the system is dilute, that is, if n|ap <C 1. In order to estimate the accuracy of this approach we will now 
calculate the first corrections to the mean-field approximation. Such corrections have been recently investigated in 
several papers as, for instance, by Timmermans, Tommasini and Huang (1997) and by Braaten and Nieto (1997). 
Here we limit the discussion to the case of repulsive interactions and large N, where analytic results can be found. In 
fact, in this limit the solution of the stationary GP equation (|3^) for the ground state density can be safely replaced 
with the Thomas-Fermi expression ( |50| ) and the energy of the system is given by E/N — (5/7)/xtf, where /^tf is the 
TF chemical potential (^. 

Let us first discuss the behavior of the ground state density. For large N one can use the local density approximation 
for the chemical potential: 

M = Wocai[n(r)] + Vc^t{r) ■ (58) 

The use of the local density approximation for fi is well justified in the thermodynamic limit N oo, lo where the 
profile of the density distribution is very smooth. Equation ( |5^ ) fixes the density profile n(r) of the ground state once 
the thermodynamic relation /xiocai('^) for the uniform fluid is known, the parameter fj, in the l.h.s. of Eq. (|5^ ) being 
fixed by the normalization of the density. For example, in a very dilute Bose gas at T = 0, one has fj-iocsnin) = gn 
and immediately finds the mean- field Thomas- Fermi result (|50|). The first correction to the Bogoliubov equation of 
state is given by the law (Lee and Yang, 1957; Lee, Huang and Yang, 1957) 



A*iocai(") = gn 



3V7r 



(59) 



which includes nontrivial effects associated with the rcnormalization of the scattering length. Using expression (|5^ 
for /iiocai, one can solve equation (Esh by iteration. The result is 



A-TTl / ^ I 

n(r) = g-^ - V,^,{v)\ - -—^ [^, - V,^,{r)f' , (60) 



with fi given by 



(l + y/naMO)) ■ (61) 



Then the energy can be also evaluated through the thermodynamic relation /_* = dE/dN, and one finds 

5 



E : 



yiVA^TF (^1 + ^V^^MO)) (62) 



where, in the second term, we have safely used the lowest order relation /ixF ~ .9"(0). In an equivalent way, results 
(|60|)-(p2|) can be derived using a variational procedure by writing the energy functional of the system in the local 
density approximation. 

Equations (|6l|)-(|62|) show that, as expected, the corrections to the mean-field results are fixed by the gas parameter 
a?n evaluated at the center of the trap. This quantity can be directly expressed in terms of the relevant parameters 
of the system: 

1 c2/5 / „ \ 12/5 

a3„(0) = i£_ (7Vi/6^) . (63) 



StT \ flho 

Inserting typical values for the available experiments, the corrections to the chemical potential and the energy turn 
out to be of the order of 1%. These corrections to the mean- field predictions should be compared with the ones due to 
finite size effects (quantum pressure) in the solution of the Gross-Pitaevskii equation [see Eqs. (|5^) and (|57|)], which 
have a different dependence on the parameters N and a/oho- One finds that finite size effects become smaller than 



the corrections given by Eqs. (61)-(|62[) when N is larger than about 10^. 

Another important quantity to discuss is the quantum depletion of the condensate. This gives the fraction of atoms 
which do not occupy the condensate at zero temperature, because of correlation effects. The quantum depletion 
is ignored in the derivation of the Gross-Pitaevskii equation. It is consequently useful to have a reliable estimate 
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of its value in order to check the validity of the theory. Also in this case we can use local density approximation 
(Timmermans, Tommasini and Huang, 1997) and write the density of atoms out of the condensate, riout(r), using 
Bogoliubov's theory for uniform gases at density n — n{r) (see for example Huang, 1987). One gets nout(r) = 
(8/3)[n(r)a'^/7r]^/^. Integration of Tiout yields the result: 



N 



y/a^n{Q) . (64) 



for the quantum depletion of the condensate. Similarly to the correction to the mean- field energy (^), this effect is 
very small (less than 1%) in the presently available experimental conditions. 

The above results justify a posteriori the use of the Bogoliubov prescription for the Bose field operators and the 
perturbative treatment of the noncondensed part at zero temperature. We recall that this situation is completely 
different from the one of superfluid ''He where quantum depletion amounts to about 90% (Griffin, 1993; Sokol, 1995). 

IV. EFFECTS OF INTERACTIONS: DYNAMICS 
A. Excitations of the condensate and time dependent GP equation 

The study of elementary excitations is a task of primary importance of quantum many-body theories. In the case 
of Bose fluids, in particular, it plays a crucial role in the understanding of the properties of superfluid liquid helium 
and was the subject of pioneering work by Landau, Bogoliubov and Feynman (for a recent discussion on the dynamic 
behavior of interacting Bose superfluids see, for instance, GrifRn, 1993). 

After the experimental realization of BEG in trapped Bose gases, there has been an intensive study of the excitations 
in these systems. Measurements of the frequency of the lowest modes have soon become available and the direct 
observation of the propagation of wave packets has been also obtained. In the meanwhile, on the theoretical side, a 
variety of papers has been written to explore several interesting features exhibited by the dynamic behavior of trapped 
Bose gases. 

Let us start our discussion recalling that for dilute Bose gases an appropriate description of the excitations can 
be obtained from the time dependent GP equation (pq) for the order parameter. This equation has been already 



used in Sec. IH for evaluating the stationary solution (/)(r) characterizing the ground state. In the low temperature 
limit, where the properties of the excitations do not depend on temperature, the excited states can be found from the 
"classical" frequencies uj of the linearized GP equation. Namely, one can look for solutions of the form 

$(r, t) ^ e-'f'*/'' [0(r) -I- u(r)e~^^* -I- v* (r)e*"*] (65) 

corresponding to small oscillations of the order parameter around the ground state value. By keeping terms linear in 
the complex functions u and v, Eq. ( p5[ ) becomes 

huju{r) = [Hq - /i + 2.g02(r)]M(r) + g(l)'^{r)vir) (66) 
-hujv{r) ^ [Ho - fi + 2g(f>'^{r)]v{r) + g(t)'^{r)u{r) . (67) 

where Hq — — (?i^/2m)V^ + Kxt(r). These coupled equations allow one to calculate the eigenfrequencies uj and hence 
the energies e — hio of the excitations. This formalism was introduced by Pitaevskii (1961), in order to investigate 
the excitations of vortex lines in a uniform Bose gas. 

This procedure is also equivalent to the diagonalization of the Hamiltonian in Bogoliubov approximation, in which 
one expresses the field operator ^' in terms of the quasiparticle operators aj and Oj through (Fetter, 1972 and 1996) 

*'(r) ^ ^K(r)a,(t) + z;*(r)a](t)] . (68) 
j 

By imposing the Bose commutation rules to the operators aj and aj, one finds that the quasiparticle amplitudes u 
and V must obey the normalization condition 

dr [u*{r)u,{r)-v*{r)v,{r)] = S.,, . (69) 

In a uniform gas, the amplitudes u and v are plane waves and the resulting dispersion law takes the most famous 
Bogoliubov form (Bogoliubov, 1947) 
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(M'= ('??^) (^ + 29n] (70) 



2m J \ 2m 

where q is the wavevector of the excitation and n — j^p is the density of the gas. For large momenta the spectrum 
coincides with the free-particle energy h^q^/2m. At low momenta Eq. (jT^) instead yields the phonon dispersion 
uj = cq, where 

(71) 

TO 

is the sound velocity. It is worth noticing that this velocity coincides with the hydrodynamic expression c = 
[{l/m)dP/dn]^^'^ for a gas with equation of state P = {l/2)gn? [see also the discussion after Eq. (fz^)]. 

In the case of harmonic trapping, an important role is played by the ratio Na/a\^o, and one expects different 
behaviors in the two opposite limits Na/aho ^ 1 and Na/ai^o 3> 1. In the first case, one recovers the excitation 
spectrum oj — UxLOx + nyUjy + UzUJz of the nonintcracting harmonic potential [see Eq. In the second case, one 

obtains a different dispersion law for the excitations of the system which are the analog of phonons [see Eq. ( |80| ) 
below] . 

The coupled equations (|66[)-(|67|) were first used to calculate numerically the excitations of trapped gases by Burnett 
and co-workers (Ruprecht et ai, 1996; Edwards et al, 1996a and 1996c). Similar calculations have been also performed 
by other authors, for both spherical and anisotropic configurations (Singh and Rokhsar, 1996; Esry, 1997; Hutchinson, 
Zaremba and Griffin, 1997; Hutchinson and Zaremba, 1997; You, Hoston and Lewenstein,1997; Dalfovo et ai, 1997a). 

For spherical traps, the solutions of Eqs. (^6|)-(|67|) are characterized by the quantum numbers n^, i and m, where 
rir is the number of radial nodes, £ is the angular momentum of the excitation and m its z component. For axially 



symmetric traps the third component to of angular momentum is still a good quantum number. In Fig. 14 we report 
the lowest solutions of even parity with m = and to = 2, obtained for a gas of rubidium atoms confined in an 
axially symmetric trap [ojx = ojy = u!±). The asymmetry parameter of the trap (A = uJz/ijJ± = VS) corresponds to 
the experimental conditions of Jin et al. (1996) and values of N up to 10* are considered. Actually the results are 
reported as a function of the dimensionless parameter Na/a± where a±_ = {muj±). The theoretical predictions are 
compared with the experimental results. In the experiments these oscillations are observed by shaking the condensate 
through the modulation of the trapping magnetic fields. The general agreement between theory and experiments is 
good and reveals the important role played by two-body interactions. In fact, in the absence of interactions, the 
eigenfrequencies would be the ones predicted by the ideal harmonic oscillator, which gives uj = 2uj± for both modes. 

Among the various excitations exhibited by these trapped gases, special attention should be devoted to the dipole 
mode. This oscillation corresponds to the motion of the center of mass of the system which, due to the harmonic 
confinement, oscillates with the frequency of the harmonic trap (this frequency can of course be different in the three 
directions). Two-body interactions cannot affect this mode because, in the presence of harmonic trapping, the motion 
of the center of mass is exactly decoupled from the internal degrees of freedom of the system. This is best understood 
by considering Eq. ( ^ ) and looking for solutions of the form 

e^^^^'^ nx,y,z + ait)) , (72) 

where, for simplicity, we have considered only oscillations along the z-axis. By a proper change of variables, z z + a, 
one finds that ( [72[ ) corresponds to an exact solution of the time dependent equation ( p5| ) oscillating with frequency 
ujz- This property holds not only in the context of the Gross-Pitaevskii equation, but is valid for any interacting 
system confined in a harmonic potential at zero as well as finite temperature, and is independent of statistics (Fermi 
or Bose). For example, such a decoupling is a well known property of shell model theory in nuclear physics (Elliott 
and Skyrme, 1955; Brink, 1957). It also exhibits interesting analogies with Kohn's theorem for electrons in a static 
magnetic field, stating that the cyclotron frequency is not affected by interactions (Kohn, 1961) [see also Dobson 
(1994) and references therein for discussions about the generalization of Kohn's theorem to the case of electrons 
confined in harmonic traps] . 

The fact that the dipole frequency is not affected by two-body interactions offers a direct test on the numerical 
accuracy of the various methods used to solve the equations of motion. On the other hand the experimental deter- 
mination of the dipole frequency turns out to be a very useful procedure to check the harmonicity of the trap and to 
determine accurately the value of the trapping frequencies. The properties of the dipole excitation in the framework 
of Bogoliubov theory have been discussed in detail by Fetter and Rokhsar (1998) [see also Kimura and Ueda (1998)]. 



Of course the coupled equations (|6^)-(67) provide a full series of solutions, with different values of the corresponding 
quantum numbers. So far experiments have provided direct information only on the low energy modes which can be 
directly excited by suitable modulation of the harmonic trap. These excitations will be further discussed in the next 
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section using the formalism of coUisionless hydrodynamic equations. States at higher energy and multipolarity are 
also important, since they characterize the thermodynamic behavior of the system, as we will see later on. 

Finally, we note that, starting from the solutions of Eqs. (|66|)-(|67|), one can also evaluate the density of particles out 
of the condensate at zero temperature (quantum depletion) by summing the square modulus of the "hole" amplitude 
V over all the excited states: nout(r) — J2j (Fetter, 1972). The results (Hutchinson, Zaremba and Griffin, 

1997; Dalfovo et al, 1997a) are in agreement with the local density estimate (|64[). 



B. Large Na/ai^o limit and coUisionless hydrodynamics 



When the number of atoms in the trap increases, the eigenfrequencies of the coupled equations (|6^)-(67) approach an 
asymptotic value. The new regime is achieved when the condition Na/a\^o S> 1 is satisfied. In this limit the excitations 
are properly described by the hydrodynamic theory of superfluids in the coUisionless regime at zero temperature. In a 
dilute gas this theory can be explicitly derived starting from the time dependent GP equation (|35|). To this purpose, 
it is convenient to write the complex order parameter $ in terms of a modulus and a phase, as follows: 



<i>(r,i) = v/;^e^^("-'*) . (73) 

The phase fixes the velocity field 

n(r,t)v(r,i) = -^($*V$ - $V$*) , (74) 



so that 



v(r,i) = -V5(r,t). (75) 



The GP equation (^5|) can hence be rewritten in the form of two coupled equations for the density and the velocity 
field: 

d 

—n + V • (vTi) = (76) 

and 

m-. + V + gn - ^V^V^ + — = . (77) 



Equation ( |76D is the equation of continuity, while ( |77| ) establishes the irrotational nature of the superfluid motion. It 
is worth noticing that, at this stage, Eqs. (|76| ) and ( [ttI ) do not involve any approximation with respect to the GP 
equation (^5|) and can be used in the linear as well as nonlinear regimes. 

If the repulsive interaction among atoms is strong enough, then the density profiles become smooth and one can 
safely neglect the kinetic pressure term, proportional to fi^ , in the equation for the velocity field, which then takes 
the form 

d ( TJlXp" \ 

TO^v + V ( Vext +.g?^+ — j = . (78) 

This result corresponds to the equation of potential flow for a fluid whose pressure and density are related by the 
equation of state P = (\l2)gr?. Equations (fz^ ) and ( ^8|) have the typical structure of the dynamic equations of 
superfluids at zero temperature (see, for example. Pines and Nozieres, 1966, Vol.11) and can be viewed as a particular 
case of the more general Landau's theory of superfluidity. According to this theory, which is valid if the relevant 
physical quantities change slowly on distances larger than the healing length, a complete description of the dynamics 
of the fluid is obtained by coupling the equation for the superfluid velocity field with a Boltzmann-type equation 
for the distribution function of elementary excitations [see, Lifshitz and Pitaevskii, 1981, §77]. At high temperature, 
when the mean free path of elementary excitations is short, one gets a system of two-fluid hydrodynamics equations. 
Conversely, at low temperature, where the role of thermally excited states is negligible, the same equations reduce to 
the hydrodynamic-type equations ( fz^ ) and (^8|) involving only the superfluid velocity. In this sense, equations can 
be referred to as the hydrodynamic equations of superfluids. They should not be confused with the hydrodynamic 
equations valid in the coUisional regime at high temperature. 
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The stationary solution of Eq. (|7g) coincides with the Thomas- Fermi density ( |50D while the time dependent equations 
( fldj ) and (|7^), after linearization, take the following simplified form: 

— 5n = V-\c\r)V5n] (79) 

where mc^{r) — dP/dn — n — Vi3xt(r), the quantity c having the meaning of a local sound velocity. 

The validity of the equation ( [79| ) is based on the assumption that the spatial variations of the density are smooth 
not only in the ground state, but also during the oscillation. In a uniform system {Vcxt = 0) this is equivalent to 
imposing that the collective frequencies be much smaller than the chemical potential. In this case, the solutions of 
( [tqI ) are sound waves propagating with the Bogoliubov velocity (|7^). Sound waves can propagate also in nonuniform 
media, provided we look for solutions varying rapidly with respect to the size of the system, so that one can assume 
a locally uniform sound velocity (Landau and Lifshitz, 1987, §67). This is possible if both the conditions qL ^ 1 and 
hq <?C mc are satisfied, where L is the size of the condensate and q is the wavevector of the sound wave. Furthermore, if 
the system is highly deformed and cigar-shaped, one can simultaneously satisfy the conditions qZ ^ 1 and qR± <^ 1, 
characterizing one-dimensional waves propagating in the z direction. Here Z and R± are the radii of the condensate 
in the axial and radial directions, respectively. In this case, one can show (Zaremba, 1998) that the sound velocity 
in the central region of the trap is given by a/ /i/2rn, instead of the usual Bogoliubov value •^//xT'^i where n = gn{0) 
and n(0) is the value of the central density. The occurrence of the extra factor 2 follows from the fact that, in the 
"one-dimensional" geometry, the sound velocity is fixed by the density averaged over the radial direction, which is of 
course smaller than its central value. 

In the experiments of Andrews, Kurn et al. (1997), one-dimensional sound waves are generated by focusing a laser 
pulse in the center of the trap. A wave packet forms in this way, propagating outwards. It is then imaged at different 
times so that the value of the sound velocity can be directly measured. In Fig. [l^ we show the observed values of c 
at different densities. The agreement with the theoretical predictions is reasonably good especially at high density. 
Possible sources of inaccuracy at low density are discussed by the same authors. The theoretical analysis of the 
propagation of wave packets and sound waves in the elongated geometry has been the object of several recent works 
(Zaremba, 1998; Kavoulakis and Pethick, 1998; Stringari, 1998). 

In these nonuniform condensates, as already said, oscillations having wavelength much smaller than the size of the 
system or, equivalently, frequency much larger than the trapping frequency Whoj propagate as usual sound waves. 



Conversely, solutions of (79) at lower frequency, of the order of Who, involve a motion of the whole system (Baym and 



Pethick, 1996). They coincide with the low energy solutions of Eqs. (|66|)-(|67|) discussed in the previous section. For a 

spherical trap these solutions are defined in the interval < r < R and have the form Sn{r) = p^*-^"''-' (r/R) Yim{S, if)) 

where p'f'"''^ are polynomials of degree 2n, containing only even powers. The dispersion law of the discretized normal 
modes is given by the formula (Stringari, 1996b) 

uj{nr, i) = cjiio(2n^ 4- 2nr£ + 3n^ + i)^^'^ . (80) 
This result can be compared with the prediction for noninteracting particles in harmonic potential: 

Uj{nr,e) =^LUho{2nr + l) (81) 

with 2rt,. + 1 = rix + ny + riz [see Eq. (H)] . Of particular interest is the case of the so called surface excitations (n^ = 0) 
for which Eq. ( pO| ) predicts the dispersion law uj — \/£ lui^q. The frequency of these modes is systematically smaller 
than the harmonic oscillator result iuiho- Notice that in the dipole case (n^ = 0, i = 1) the prediction ( |80| ) coincides 
with the oscillator frequency, in agreement with the general considerations discussed in the previous section. 

As concerns compressional modes (n^ 7^ 0), the lowest solution of ( |79| ) is the monopole oscillation, also called the 
breathing mode, characterized by the quantum numbers = 1 and i = 0. The formula (|80| ) gives the result \/5 Who, 
higher than the corresponding prediction of the noninteracting model, which gives 2Li;ho- 

For a fixed value of N the accuracy of prediction (|8^) is expected to become lower and lower as and £ increase. 
In fact, for large rir and i the oscillations of the density have shorter wavelength and neglecting the kinetic energy 
pressure in (^^ is no longer justified. In analogy with the case of uniform Bose gases, the condition for the applicability 



of the hydrodynamic theory of superfluids is expected to be huj < fi. However, as discussed in Sec. [VE, more severe 
restrictions are imposed when one considers surface excitations. 

The result ( |80| ) reveals that, in the Thomas- Fermi limit Na/a\^o ^ 1, the dispersion relation of the normal modes of 
the condensate has changed significantly from the noninteracting behavior, as a consequence of two-body interactions. 
However it might appear surprising that in this limit the dispersion does not depend any more on the value of the 
interaction parameter a. This differs from the uniform case where the dispersion law, in the corresponding phonon 
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regime, is given hy lu = cq and depends explicitly on the interaction through the velocity of sound. The behavior 
exibited in the harmonic trap is well understood if one notes that the values of q are fixed by the boundary and vary 
as 1/L where L is the size of the system. While in the box this size is fixed, in the case of harmonic confinement 
it increases with N due to the repulsive effect of two-body interactions: L ^ {Na/aho)'^^^{'muJho)~'^^^- On the other 
hand the value of the sound velocity, calculated at the center of the trap, is given by c = (A^a/aho)2^^(who/™)^^2 and 
also increases with N. One finally finds that in the product cq both the interaction parameter and the number of 
atoms in the trap cancel out, so that the collective frequency is proportional to the oscillator frequency cjho- 

The results for the spherical trap can be generalized to the case of anisotropic configurations. Let us consider the 
case of a harmonic oscillator trap with axial symmetry along the z axis. In this case the differential equation 
takes the form 

m^Sn = V .{[/.- I {ulrl + ulz^)] V6n] (82) 

where we have used mc^{Y) — /i — ycxt(r). [We notice, in passing, that the corresponding Eq. (21) in Stringari (1996b) 
was misprinted, since it contains the chemical potential counted twice.] 

Because of the axial symmetry of the trap the third component m of the angular momentum is a good quantum 
number. However, in contrast to the spherical case, the dispersion law depends on m. Explicit results are available in 
some particular cases. For example, quadrupole solutions of the form 5n ~ r'^Y2m{0, (f>) satisfy Eq. ( ^2[ ) for m = ±2 
and m = ±1. The resulting dispersion laws are: 

Lu'^{i = 2,m = ±2) = 2ujl (83) 

and 

uj'^{e^2,m^±l) =ljI+luI. (84) 

Conversely the £ = 2,m = mode is coupled to the monopolc £ = excitation and the dispersion law of the two 
decoupled modes is given by (Stringari, 1996b) 

uj^{m = 0) = 2ujI + ^ujI ip iy^9wf - 16^2^^ + IGcji . (85) 

When = one recovers the solutions for the quadrupole and monopole excitations in the spherical trap. The 
occurrence of analytic solutions for the excitation spectrum, like Eqs. ( po|) and (|8^)-(|85|), is the result of nontrivial 
underlying symmetries of the Hamiltonian that have been exploited by Fliesser et al. (1997). The result ( |85| ) can 
be generalized to a triaxially deformed trap of the form (^. In this case, the collective frequencies are given by the 
solution of the equation 

From Fig. p^one can see that the experiments at JILA do not fully fall in the asymptotic Na/a^o ^ 1 regime, where 
the frequencies are given by Eqs. (|8^^. Conversely the experimental results obtained on sodium vapors at MIT 
(Stamper-Kurn et ai, 1998c) represent a very clear example of excitations belonging to the Thomas-Fermi regime. 
In this experiment the magnetic trap is highly asymmetric, with A = Wz/wx = 17/230 (cigar-shaped geometry). 
Furthermore the number of atoms is very high, so that the condition Na/uho S> 1 is well satisfied and the energies 
of the collective oscillations along the axial direction are much smaller than the chemical potential, ^ « 200wz. This 
explains the excellent agreement between the observed frequency for the lowest axial m = mode of even parity 
{lu/uJz — 1.569(4)) and the theoretical prediction (lj/uJz = a/5/2 — 1.581) given by Eq. ( |85| ) with uj^ ^ uj±. In 
Fig. ^ we show the oscillations observed in the MIT experiment (see also Fig. ||). These measurements correspond 
to nondestructive in situ images of the oscillating condensate, while the ones at JILA (Jin et ai, 1996 and 1997), as 
well as the first experiments carried out at MIT (Mewes et al., 1996b) were taken after switching-off the trap and 
letting the gas expand. 

For highly deformed traps it is possible to obtain simple analytic results also for the excitations with higher quantum 
numbers. For example, in the case of cigar-shaped traps {coz ^ ojx) one finds the dispersion law (Fhesser et al., 1997; 
Stringari, 1998) 

u;^{k)^h{k + 3)Lul (87) 

where k is the relevant quantum number characterizing the spatial shape of the density oscillation 6n{z) — {z^ + 
^^fe-i _ y Equation (|87|) is valid if uj{k) <^ uj±. It includes, as special cases, the dipole {k = l,uj = uj^) and 
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"quadrupole" {k = 2,uj = ui^) modes already discussed. It also permits one to understand the transition 

between the discretized (small k) and "continuum" {k ^ 1) regimes, through the identification k = qz where q is 
the wave vector of "one-dimensional" phonons propagating with sound velocity ^ iijlvn. As already discussed, these 
phonons can be considered one-dimensional only if the conditions qZ 3> 1 and qR\_ ^ 1 are satisfied. The first 
condition implies large values of fc, the second one is equivalent to imposing lo <C 

Analogously for disk-shaped traps ^ w^) the dispersion law of the lowest modes takes the analytic form 
(Stringari, 1998) 

p/ N /4 , 4 \ , 

Lo (rir-, m) — I -n^ -I- —nj-m 2nr + ml , (88) 

where = 0, 1, .. is the number of radial nodes and m is the z-component of the angular momentum. 



C. Sum rules and collective excitations 



In the previous section we have discussed the excitations of the condensate when atoms interact with repulsive 
forces (a > 0). In the opposite case of attractive interactions (a < 0), one expects a different behavior. For example, 
interesting effects can originate from the fact that the system becomes more and more compressible when approaching 
the critical number, A^cr, for collapse. In terms of the excitation spectrum, this means a lowering of the frequency of 
the monopole oscillation. For repulsive forces, we have previously discussed the Thomas-Fermi Na/uho 3> 1 limit, in 
which the time dependent GP equation takes the form of the zero temperature hydrodynamic theory of superfluids. 
When the interaction is attractive, the large A^ limit is never reached, since the collapse occurs at A^a/oho of the 
order of 1, and one has to solve numerically the GP equation (see, for example, Dodd et ai, 1996) or use different 
theoretical schemes, as shown in the following. 

A useful physical insight on the behavior of collective oscillations for both positive and negative a can be obtained 
using the formalism of linear response and sum rules [see, for instance, Bohigas, Lane and Martorell (1979), Lipparini 
and Stringari (1989)]. This approach allows one to evaluate the energy weighted moments, nip = SF{E)EPdE, of 
the strength distribution function (dynamic form factor) associated with a given operator F: 

SF{E) = J2\{j\FmHiE-E,o) (89) 
j 

where the quantity Ejo = {Ej — Eq) is the excitation energy of the eigenstate |j) of the Hamiltonian. Consequently, 
the method provides information on the dynamic behavior of the system. Quantities like TOp+i/mp or (mp+2/mp)^/^ 
correspond to rigorous upper bounds for the energy of the lowest state excited by the operator F. They are close to 
the exact energy when this state is highly collective, that is, when the strength distribution is almost exhausted by a 
single mode. This is often true in the case of trapped gases, as we will see below. 

A major advantage of sum rules is that they can be often evaluated in a direct way, avoiding the full solution of the 
Schrodinger equation for the eigenstates of the Hamiltonian. For example, using the completeness of the eigenstates 
\j), the energy weighted moment, mi, can be easily transformed into the calculation of commutators involving the 
operator F and the Hamiltonian: 

m, = ^{0\[F\[H,F]]\0) . (90) 

Furthermore, if the operator _F depends only on spatial co-ordinates then only the kinetic energy gives a contribution 
to mi, whose calculation becomes straightforward. In a similar way, one can write the cubic energy weighted moment, 
ma, in the form — H], [H, [H, F]]]). Unlike mi and m^, the inverse energy weighted moment, m^i, cannot 

be expressed in terms of commutators; it can be however written in the useful form 

m-i = ^ X , (91) 

where x is the linear static response of the system. 

Let us first consider the case of compressional modes. The natural monopole operator is given by the choice 
F = and, from (|90|), one gets the result mi — 2N'h'^{r^)/m for the energy weighted sum rule. Furthermore, in 

the monopole case one can easily evaluate also the inverse energy weighted sum rule through Eq. (91). In fact, the 
static response xm (monopole compressibility) is fixed by the linear change d{r^) = exM of the mean square radius 
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induced by the external field — er^. Adding this field to the Hamiltonian is equivalent to renormalizing the trapping 
harmonic potential and hence, for isotropic confinement, one can write 

XM = ■ 92 



can express exactly the derivative d{r^) /duj^^ in terms of the square radius, (r^), and its derivative with respect to 



where Who is the frequency of the harmonic oscillator. Using the properties of the Gross-Pitaevskii equation (|3q), one 
N. One finds (ZambeUi, 1998) 

N 



XM 



moj. 



ho 



(93) 



where the term depending on the derivative arises from two-body interactions; in the case of an ideal gas, this term 
vanishes and the mean square radius, (r^) = (8/2)0^^, is independent of N . 

Using the moments mi and to_i one can define an average excitation energy, Tioj, through the ratio 



yielding the useful result (ZambeUi, 1998) 



{nuf - ^ , (94) 



^M=^^lo ,,._Vj^ ■ (95) 

V / 2 dNV 



In the noninteracting case, one recovers lom — '^'-^ho- When Na/aho is large and positive, the Thomas-Fermi ap- 
proximation ( pO| ) for the density provides the analytic behavior of the radius, (r^) cx TV^^^, and hence the result 
ujm = tJho already discussed in the previous section. For negative a and close to the critical size A^cr, the 
monopole frequency goes to zero because the compressibility of the system becomes larger and larger. Actually the 
A^-dependence of ujm ne ar A'cr can be determined analytically. For example, using the Gaussian variational procedure 



developed in Sec. IIIC , one finds the result ((r^) — (r^)cr) = (?'^)cr-\/8/5(l — N/Ncr)^/"^, where (r^)cr is the square 
radius of the condensate at the critical value A^cr- As a consequence of this peculiar A^ dependence, the monopole 
compressibility diverges near A'cr and the monopole frequency vanishes as (Singh and Rokhsar, 1996; ZambeUi, 1998; 
Ueda and Leggett, 1998) 

/ W \ 1/4 

c^M-c^ho(160)i/M 1-— j . (96) 

By using the numerical solution of the Gross-Pitaevskii equation to calculate the A-dependence of the square radius 
one finds a slightly smaller value for the numerical coefficient in (^6|), namely 3.43 instead of 3.56. It has been suggested 
that the behavior of the monopole frequency near A'cr might play an important role in the decay mechanism of the 
condensate for A^ very close to A^cr, due to quantum tunnehng (Ueda and Leggett, 1998). 



In Fig. |17| we show the frequency lom obtained from Eq. (|95| ) as a function of the parameter Na/aho (solid line). 
The square radius (r^) has been calculated by solving numerically the stationary GP equation (p9) in a spherical 
trap. As already said, the ratio ( p^ ) between moments of the strength distribution function Sf{uj) provides a rigorous 
upper bound to the lowest monopole frequency. The comparison with the numerical solutions of the time dependent 
GP equation (circles) shows that the sum rule estimate actually gives an excellent approximation to the collective 
frequency for both positive and negative values of a, practically indistinguishable from the exact result. This means 
that the strength distribution of the monopole operator F almost coincides with a 5-function located at the energy 
of the lowest compressional mode. For the same reason, also the ratio m^/mi turns out to be very close to mi/m_i 
(Zambelh, 1998). 

Unlike the monopole frequency, the quadrupole frequency increases with A^ when a < 0, due to the increase of the 
kinetic energy of the condensate. This behavior is well understood by calculating the quadrupole frequency through 
the ratio [hLoY — "13/mi, where mi and are the energy and cubic energy weighted moments for the natural 
quadrupole operator F = '"^^am- By explicitly working out the commutators of the two sum rules one finds the 

following result for the quadrupole frequency (Stringari, 1996b): 



w^^2< (1 + ^) . (97) 
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In the noninteracting gas one has E'kin = E^o and (|9^) gives the harmonic oscillator result ujq = 2(jJho- In the 
Thomas- Fermi limit Na/aho ^ 1, the kinetic energy term is negligible and one finds the value uiq = Who, while 
for negative a the kinetic energy term is larger than _Eho and one finds an enhancement of the quadrupole frequency. 
The numerical results are reported in Fig. |l^. Also in the quadrupole case the sum rule estimate (p7|) turns out to be 
very close to the exact numerical solution of the linearized time dependent GP equations (|6^)-(p7|), indicating that 
the lowest quadrupole mode almost exhausts the strength distribution of F, as already found for the monopole mode. 

We have here applied the sum rule approach to the case of spherical traps, but the same analysis can be easily 
generalized to more complex geometries, the main physical arguments remaining unchanged. Calculations of sum 
rules for axially symmetric traps have been carried out by Kimura and Ueda (1998), finding accurate predictions for 
the collective frequencies. 



D. Expansion and large amplitude oscillations 

So far we have discussed the behavior of normal modes of the condensate and sound propagation. It is also interesting 
to investigate nonlinear features associated, for example, with the dynamics of the expansion of the gas, following the 
switching-off of the trap, as well as with the frequency shifts of large amplitude oscillations. The dynamics of the 
expansion is an important issue because much information on these Bose condensed gases is obtained experimentally 
from images of the expanded atomic cloud. This includes in particular the temperature of the gas (which is extracted 
from the tail of the thermal component), the release energy and the aspect ratio of the velocity distribution. Nonlinear 
phenomena are also crucial in the analysis of the large amplitude oscillations which are produced and detected in 
current experiments. Phenomena like mode coupling, harmonic generation, frequency shifts and stochastic behavior 
may become interesting subjects of research in these systems. 

From the theoretical viewpoint one can again attack the problem starting from the time dependent Gross-Pitaevskii 
equation. Indeed, the GP equation ( ^5|) for the order parameter of the condensate can be applied to the nonlinear 
regi me an d it is important to check the validity of its predictions through a direct comparison with experiments. In 
Sec. IV A we have linearized this equation in order to obtain the coupled equations (|66|)-(|67|) for the excitations. The 



numerical solution in the nonlinear regime is also feasible (Holland and Cooper, 1996; Holland et ai, 1997; Ruprecht 
et ai, 1996; Smerzi and Fantoni, 1997; Morgan et ai, 1998; Brewczyk et ai, 1998). For instance, Holland et al. 
(1997) obtained results for the density and energy of the expanding gas of ^^Rb in good agreement with the first 



measurements at JILA. In the inset of Fig. 10 their results for the axial and radial widths are plotted as a function 
of expansion time. 

When the number of atoms in the trap is large, the time dependent GP equation ( ^ ) reduces to the equations 
of continuity ( [76| ) for the density and the Euler equation (|7^) for the velocity field. These equations can be used 
to investigate nonlinear phenomena in a simplified way. Let us write the external potential in the form Vcxtij) — 
(m/2) a^j^rf, with = x,y,z. In general, the trapping frequencies can depend on time, uji — LOilt)] their static 
values, ujQi = a;i(0), fix the initial equilibrium configuration of the system, corresponding to the Thomas- Fermi density 
(^). One can easily prove that the equations of motion admit a class of analytic solutions having the density in the 
form 

n(r, t) = ao(0 - a:,{t)x^ - ay{t)y'^ - a^{t)z'^ , (98) 
within the region where n(r, t) is positive and n(r, t) = Q elsewhere, and the velocity field as 

v(r,i) = ]^V[a,{t)x^ + ay{t)y^ + a,{t)z^] . (99) 

These results, combined with Eq. ([7^), allow one to obtain an explicit expression for the order parameter ([7^). In 
particular, its phase S takes the form 

S{v,t) = ^[a,{t)x^ + ay{i)y^ + a,{t)z^] . (100) 



Notice that, while the velocity field (|99|) is governed by the classical equations (|7^) and (|78|), the phase of the order 
parameter depends explicitly on the Planck constant Ti. 



The results (98) and (^9|) include the ground state solution (^0|) in the Thomas- Fermi limit. This is recovered by 
putting = and = muj^ i /{ 2q), with i = x,y, z, while oq — fJ^/g- In general, one can insert expressions ( p8[ ) and 
( p9| ) into the equations ([76])- ([Tq), getting six coupled differential equations for the time dependent coefficients ai(t) 
and ai{t), while the relation aq — {\hN/^iT)'^/^{axayazY/^ is fixed, at any time, by the normalization of the density 
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to the total number of particles. Instead of writing these equations, we note that the assumptions ( |98| ) and for 
the density and velocity distributions correspond to assuming a scaling transformation of the order parameter. This 
means that, at each instant, the parabolic shape of the density is preserved, while the classical radii Ri, where the 
density (|9q) vanishes, scale in time as 



R,it) = i?,(0)5,(<) = hit) . (101) 



equations (76)- ([781) then give at — hi/bi and 



The relation among the coefficients of Eq. ( P8|) and the variables bi is found to be — muiQ^/ [2gbxbybzb'l) and the 

-.^6.-n4^ = o. (102) 



bibxbyb^. 



These are three coupled differential equations for the scaling paramete rs bj ( t) , which in turn give the time evolution of 



the classical radii, Ri{t), of the order parameter. The second term in (102) comes from the confining potential, while 
the third one originates from the atom-atom interaction. Equations (102) have been derived and used by different 
authors (Kagan, Surkov and Shlyapnikov, 1996, 1997a and 1997b; Castin and Dum, 1996; Dalfovo et ai, 1997b and 
1997c). Their major advantage is that they are ordinary differential equations, very easy to solve, giving results close 
to the solutions of the time dependent GP equation in most situations. Kagan, Surkov and Shlyapnikov (1996) have 
shown that in 2D the scaling transformation of the order parameter, starting from the stationary configuration at the 
initial time, a ctua lly corresponds to an exact solution of the GP equation. 



Equations (102) can be used to simulate the expansion starting from a gas in equilibrium in the trap, by dropping 
at a certain time, t = 0, the term linear in bi associated with the confining potential. For an axially symmetric trap 
one can define b± = b^ — by and introduce a dimensionless time r = uj^t, with ll>j_ = luqx — ^Oy = ^~^^0z- Then 



Eqs. (102) take the form 



d2 1 , d2 A2 

By solving these equations, one can look, for instance, at the time evolution of the aspect ratio R±/Z ~ Xb±/bz. 
When T is large, both b± and b^ increase linearly with r and the parameters a± and a^, characterizing the velocity 
field (|99|), behave as 1/t, consistently with the classical equation of motion for free particles, v = r/t. In Fig. |l^ 
we show the results of this calculation in two cases where accurate experimental data are available: atoms of *^Rb 
released from a trap with uj± = 27r x 247 Hz and ujz = 2tt x 24 Hz at Konstanz (Ernst, 1998b), and sodium atoms 
released from a trap with ujj_ — 2ti y. 248 Hz and = 27r x 16.23 Hz at MIT (Stamper-Kurn and Ketterle, 1998). 
Both traps are cigar-shaped and the number of atoms is large enough for applying the Thomas- Fermi approximation. 
The agreement between theory (solid lines) and experimen ts (p oints) is remarkable. It is also worth mentioning that, 
as shown by Castin and Dum (1996), the two equations ( |103| ) can be solved analitycally for A ^ 1, leading to the 
useful expressions 



b^{T) = VT+^ (104) 
fez(r) = 1 + A2[r arctanr - In Vl + t2] . (105) 

The corresponding aspect ratio is plotted in Fi g, [l^ as a dashed line. As one can see, the analytic small-A limit 



pract i cally coincides with the exact solution of (103) for the two traps here considered. Moreover, from expressions 
(104)-( |l05| ) one also gets the asymptotic value limT^oo(^_L/^) = 2/(7rA). In the case of Fig. |l8|, this asymptotic limit 
is approximately 6.5 and 9.7 for the Konstanz and MIT data, respectively, but is far from being attained even after 
tens of ms; in fact, it takes a relatively long time to reach the regime of constant speed for the motion along the 
direction of weaker initial confinement, due to the slow acceleration induced by the mean-field potential. 

The general agreement between theory and experiments in Fig. |l^ is even better appreciated if one considers the 
predictions for the expansion of noninteracting particles. The aspect ratio obtained from the dispersion of a free atomic 
wave packet is represented by the two dot-dashed lines. The asymptotic limit for r ^ oo is \~^^^. The comparison 
with the behavior of the interacting gas shows, once again, the important role of the atom-atom interaction. 

The same formalism allows one to calculate the time evolution of the various contributions to the release energy of 
the condensate. In terms of the scaling parameters bi the release energy takes the form 
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1 viL' 



(106) 

This quantity is conserved during the expansion. At i = 0, when bi — 1 and bi — 0, the release energy is equal to 
(2/7) fi. During the expansion the mean-field energy (first term i n th e bracket) is converted into kinetic energy (second 
term). After a certain time, which can be estimated using Eq. ( |l06| ), the mean-field energy becomes negligible, since 
the system is mor e an d more dilute, and the expansion proceeds at constant speed in each direction. 

The same Eqs. ( |102| ) allow one also to study the effects of a sinusoidal driving force which simulates the modulation 
of the con fining potential used in experiments to generate collective modes in the trap. In the small amplitude limit, 



Eqs. (102) yields the frequencies of the normal modes in the regime of collisionless hydrodynamics already discussed 



in Sec. IV B . In particular for axially symmetric traps, expression (|98| ) includes the lowest m — and m — 2 modes, 
which have been measured experimentally, while other modes could be investigated by adding terms in xy, xz and 
yz. When the amplitude of the oscillations grows, the frequency of the modes can shift and the modes themselves can 
couple. For the experiments carried out at JILA and MIT (Jin et ai, 1996 and 1997; Stamper-Kurn et ai, 1998c) both 
the frequency shift and the mode coupling are small. An example is given in Fig. |l^, where we show the frequency 
of the lowest m = mode, observed in the cigar-shaped MIT trap (Stamper-Kurn et a t, 1 998c; Stamper-Kurn and 
Ketterle, 1998), as a function of its amplitude. The solid line is the prediction of Eqs. ( |l02| ). The overall agreement 
is good for both the zero amplitude frequency and its shift. It is worth recalling that the theory has no fitting 
parameters. 

An effect which deserves to be mentioned is the large enhancement of nonlinear effects for special values of the 
asymmetry parameter A. First, one notes that the frequencies of the collective modes depend on the shape of the 
trapping potential and hence on A. For certain values of this parameter, it may happen that different modes have 
the same frequency; this has been shown to occur, in the linear regime, for |m| > 2 by Ohberg et al. (1997) and 
a systematic investigation of the level crossing has been done by Hutchinson and Zaremba (1997). In the nonlinear 
regime, one finds strong mode-coupling via harmonic generation when the frequency of a mode becomes equal to 
the one of the second harmonics of other modes (Dalfovo et ai, 1997b; Graham et ai, 1998). The conditions for 



this degeneracy can be found numerically from Eqs. (102). In the limit of small amplitude one can also expand the 



solutions finding analytical results (Dalfovo et ai, 1997b). In particular, one gets a quadratic shift in the form 

ij (A) ^ ij{0)[l + S{\)A^] (107) 

where A is the relative amplitude of the oscillation and (5(A) is an analytic coefficient depending on the anisotropy of 
the trapping potential and on the mode considered. For instance, in the case of the m = 2 mode, one finds: 



The divergence at A = 1/16/7 is due to the degeneracy between the frequencies of the high-lying m = mode and 
the second harmonic of the m = 2 mode. In this case, it is difficult to drive the system in a pure mode and, even 
for relatively small amplitudes, the motion is rather complex and the resulting trajectories can exhibit a chaotic-like 
behavior. The coefficient (5(A) can be calculated also for other modes. For the low-lying m = mode, for instance, 
similar divergences are found when A — (Vl25 ± \/29)/\/72 (i.e., A ~ 0.683 and A w 1.952). They occur because the 
frequency of the high-lying mode becomes equal to the second harmonics of the low-lying mode. It would be very 
interesting to study experimentally the system in these conditions. 

As a final remark we note that, by means of a variational appro ach based on Gaussian wave functions, Perez-Garcia 
et al. (1996, 1997) have derived equations of motion of the form ( |102| ), but with an additional term included, propor- 
tional to 1/fof, accounting for the quantum pressure in the Gross-Pitaevskii equation. Even though the equilibrium 
configuration in the Thomas- Fermi regime Na/aho S> 1 is not exactly recovered, because of the Gaussian ansatz, 
these equations represent a good approximation to the GP equation for finite N, interpolating between the noninter- 
acting and strongly interacting systems. The frequencies of the lowest m = and m = 2 modes, calculated with the 
parameter of the JILA trap, differ from the exact solutions of (|6^)-(^7|) by less than 1% over all the relevant range 
of N . Since this method includes quantum pressure effects, it can be used to explore the nonlinear dynamics of the 
gas also in the case of attractive forces (Perez-Garci'a et al., 1997). 



E. Density of states: collective vs. single-particle excitations 

In the previous sections we have discussed several features of collective excitations pointing out the crucial role played 
by two-body interactions. We may ask whether these collective modes are relevant for the statistical properties of 
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these many-body systems. One knows, for instance, that the thermodynamic behavior of superfluid ^He is dominated 
by the thermal excitation of phonons and rotons up to the critical temperature. For the trapped gas the situation 
is very different. First the system is very dilute and one expects that the effects of collectivity shall be less relevant 
except at very low temperature. Second the harmonic confinement leaves space for excitations of single-particle nature 
which actually dominate the thermodynamic behavior even at low temperature. 

The simplest way to understand the role of these single-particle excitations is to look at the spectrum obtained by 
solving numerically the two Bogoliubov-like equations (|6^)-(|67|). In Fig. ^ we show the eigenstates evaluated for a 
condensate of 10* atoms of *^Rb in a spherical trap (Dalfovo et ai, 1997a). Each state, having energy e and angular 
momentum £, is represented by a thick solid bar. For a given angular momentum, the number of radial nodes, i.e., 
the quantum number n^, increases with energy. 

By looking at the eigenstates at high-energy and multipolarity in the spectrum of Fig. |2^, one notes that the 
splitting between odd and even states is approximately ?iWho and the spectrum resembles the one of a 3D harmonic 
oscillator. Actually, the states with the same value of {2nr + 1) would be degenerate in the harmonic oscillator case, 
while here they have slightly different energies, the states with lowest angular momentum being shifted upwards as 
a results of the mean-field produced by the condensate in the central region of the trap. Indeed the high-energy 
part of the spectrum is expected to be well reproduced by a single-particle description in mean-field approximation. 
The single-particle picture is obtained by neglecting the coupling between the positive (u) and negative [v) frequency 
components of the order parameter ( |65| ) in the Bogoliubov-type equations (|66|)-(|67|), which is responsible for the 
collectivity of the solutions. This corresponds to setting w = in Eq. (|6^), which then reduces to the eigenvalue 
problem (i?sp ^ l^)'^' = Tt-lou, for the single-particle (sp) Hamiltonian 

= -{f,^/2m)V^ + yext(r) + 2gn{v) . (109) 

In this case, the eigenfunctions u{r) satisfy the normalization condition J dr u*(r)uj{r) ~ Sij. This approximation is 
directly related to Hartree-Fock theory, as we will discuss in Sec. |VB| . 

Once the condensate density and the chemical potent ial a re calculated from the stationary GP equation (|39|), the 



single-particle excitation spectrum of the Hamiltonian (10£) can be easily calculated. The eigenstates are shown as 



dashed orizontal bars in Fig. gO|. One sees that the general structure of the spectrum is very similar to the one 
obtained with the Bogoliubov-type equations (|6^)-(|67|) apart from the states with low energy and multipolarity. The 
lowest levels, with energy well below ^ and small angular momentum, are in fact the collective modes discussed in 
the previous sections (for instance, the lowest states with Ur — are the monopole, dipole, and quadrupole modes 
for which the theory, in the limit Na/a\^o ^ 1? predicts e = V5, 1, and v^, respectively, in units of Who)- The single- 
particle spectrum, which does not account for collective motion of the condensate, fails to describe these states. It is 
worth noticing, however, that even below there are many states, with relatively high ^, which are well approximated 



by the single-particle Hamiltonian (109). Actually, the numerical analysis reveals that, for these states, the condition 



\v\ <C |u| is well satisfied (You, Boston and Lewenstein, 1997; Dalfovo et al, 1997a). These excitations are mainly 
located near the surface of the condensate, where Hsp has a minimum. The existence of such a minimum is evident 
in the large TV limit, where the Thomas-Fermi approximation for the condensate density is accurate. In this case, one 
has 

i/,p - ^ = -(aV2m)V2 + Kntofjr^ - R^\ , (110) 

where R = [2ii/{mui^^)]^/^ is the classical radius of the condensate and we have taken, for simplicity, a spherical trap. 
Of course, for finite values of N the minimum of the single-particle potential is rounded. 

The fact that the Bogoliubov-type spectrum exhibits states of single-particle nature localized near the surface 
represents an important difference with respect to the uniform Bose gas, where no single-particle states are present at 
energy lower than the chemical potential. The transition between the collective and single-particle character can be 
understood in terms of length scales. In fact, an excitation inside the condensa te can no longer be phonon-like when 



its wavelength is of the order of, or shorter than, the healing length ^ [see Sec. [II B and Eq. ([16[)]. This happens for 



states with a large number of radial nodes and energy larger than fj,. Conversely, for st ates localized mainly at the 



surface, the appropriate length scale is the surface thickness d, introduced in Sec |[IID| [see Eq. (|5^)]. In this case, 
excitations cannot be collective, and hence cannot be described by the equations of coUisionless hydrodynamics, if 
their wavelength is smaller than d. This happens when their angular momentum is larger than i ^ R/d ^ N^^^^. 
This critical value of ^ corresponds to an energy of the order of ii{aho/R)^^^, so that the transition from the collective 
to the single-particle behavior occurs at energies smaller than fi in states with high multipolarity. These states can 
be viewed as atoms rotating in the outer part of the condensate (Lundh, Pethick and Smith, 1997; Dalfovo et at, 
1997a). 
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In order to discuss the relevance of single-particle excitations in the statistical behavior of these trapped Bose gases 
it is useful to evaluate the density of states. For a finite system one can easily count the number of available states, 
with energy e' and angular momentum £, each one multiplied by its degeneracy (2£ + 1), up to a given energy e: 

7V(e)=^(2^+l). (Ill) 



The density of states is the derivative of (|lll| ). In Fig. ^we show the quantity N{e) obtained by summing the levels 
of the two spectra of Fig. ^ The agreement between the results of the Bogoliubov-type equations (solid circles) and 
of the single-particle theory (open circles) is remarkable even at low energy, indicating that the effects of collectivity 



are not relevant in the sum (111). Indeed, the number of states which are badly reproduced by the single-pa rticle 



Hamiltonian is small and their degeneracy factor {2i + 1) is also small, so that their contribution to the sum (111) 
is negligible. The effects of two-body forces on the density of states are nevertheless sizable, as emerges from the 
comparison with the prediction of the noninteracting model (dashed line). We also report the results obtained using 
the dispersion relation ( |80| ) for the excitations in the Na/aho S> 1 limit. This gives a poor approximation for N(s), 
revealing that the hydrodynamic picture becomes completely inadequate for excitation energies of the order of /i. 



The number of states N{e) associated with the single-particle Hamiltonian (109) can be also calculated using the 



semiclassical approximation. In this case one counts the available states through a simple integration over phase space 

Nie) = 1^ de' I -^Sie' - e^^p, r)) , (112) 

where e^P(p, r) — /2m + VJ;xt(r}_f 2g7i(r) — /i is the semiclassical energy corresponding to the Hamiltonian H^-p — fJ,- 
In Fig. |2l| the prediction of Eq. (|ll2| ) is shown as a solid line. The semiclassical approximation is expected to be valid 
only for e 3> huho- However, the low energy states which are not reproduced by this approximation give a negligible 
contribution to N{e) and the semiclassical prediction is practically indistinguishable from the Bogoliubov spectrum 
in the whole range of energies. 

As we will see in Sec. [V], the relevance of single-particle excitations in determining the density of states makes 
Hartree-Fock theory and the semiclassical approximation very effective tools for the investigation of the thermody- 
namic properties of these trapped gases as well as their dynamic behavior at finite temperature. 



V. EFFECTS OF INTERACTIONS: THERMODYNAMICS 



A. Relevant energy scales 



The occurrence of Bose-Einstein condensation is revealed by an abrupt change in the thermodynamic properties of 
the system below the critical temperature. In the presence of harmonic trapping a sharp peak appears in both the 
density and velocity distributions superimposed on the broader distribution of the thermal component. By further 
lowering the temperature the height of the condensate peak increases, while the tails of the thermal component are 
reduced, until they completely disappear at very low temperatures, as shown in Fig. At the transition, the 
temperature dependence of the energy shows a sudden change in slope which reflects the occurrence of a maximum 
in the specific heat. 

In Sect. H we have discussed the thermodynamic behavior of the noninteracting gas. In this model BEC takes place 
below the critical temperature fc^T^ = HuJhoiN / C{3)y . The fraction of atoms in the condensate and their energy 
obey the simple laws Nq/N = 1 — (T/T^Y ^-^id E cx T"', respectively [see Eqs. ( p^ ) and (p^)]. A major question is to 
understand whether the predictions of the ideal gas are adequate and under which conditions the effects of interactions 
become sizable. This is the main purpose of the present section. 

The effects of two-body interactions in a dilute Bose gas are expected to be significant only in the presence of the 
condensate, since only in this case can the density become relatively high due to the occurrence of the peak in the 
center of the trap. A first important consequence of repulsive forces is the broadening of the condensate peak. This 
effect, already discussed in Sect. HI at zero temperature, provides a dramatic change in the density distribution also 
at finite T and its experimental observation is an important evidence of the role played by two-body forces. The 
opposite happens in the presence of attractive forces which produce a further narrowing of the peak and a consequent 
increase of the peak density. In the following we will mainly discuss the case of systems composed by a large number 
of particles interacting with repulsive forces. 

Let us discuss the effects of a repulsive interaction by estimating the relevant energies of the system. At zero 
temperature the interaction energy per particle can be simply estimated using the Thomas-Fermi approximation 
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Eint/N = {2/7)n where fi = (l/2)?i(jJho(15A^a/aho)^^^ is the value of the chemical potential [see Eq. (51)]. It is useful 
to compare Ei^t/N, or equivalently ^, with the thermal energy ksT. If ksT is smaller than fi, then one expects 
to observe important effects in the thermodynamic behavior due to interactions. If instead /cbT is larger than fx, 
interactions will provide only perturbative corrections. Thus for repulsive forces the chemical potential provides an 
important scale of energy lying between the oscillator energy and the critical temperature: Tiwho < fJ, < ksT^- A 
useful parameter is the ratio 



7/ 



-ii-=af^V6^]''' (113) 



between the chemical potential calculated at T = in Thomas-Fermi approximation and the critical temperature for 
noninteracting particles in the same trap. Here a = 15^/^(^(3))^/'^/2 ~ 1.57 is a numerical coefficient. If one uses the 
typical values for the parameters of current experiments, one finds that 77 ranges from 0.35 to 0.40. Thus, one expects 
that interaction effects will be visible also at values of T of the order of T°. 

It is worth discussing the dependence of the parameter 77 on the relevant parameters of the system. First one 
should point out that this dependence is different from that of the interaction parameter TVa/aho already introduced 



in Sect. [II to account for the effects of two-body interactions in the GP equation for the condensate. The parameter 
Na/ a\io determines the value of the chemical potential in units of the oscillator energy, while rj fixes it in units of the 
critical temperature. This brings a different dependence of rj on N which turns out to be very smooth (77 ~ N^^^^). 
Thus, in order to change the value of this parameter, and consequently the effects of interactions on the thermodynamic 
behavior, it is much more effective to modify the ratio a/aho between the scattering and oscillator lengths rather than 
the value of N. 

Another important feature of the parameter rj is that it can be expressed explicitly in terms of the traditional "gas 
parameter" a^n, through the relation 77 = 2.24[a'^7iT=o(0)]^^^ [see Eq. (p3|)]. Notice that in this formula 7tT=o(0) is 
the density at the center of the trap evaluated at zero temperature. Due to the 1/6-th power entering this relation, 
the value of rj can be easily of t he or der of 1 even if the gas parameter is very small. For example, taking a'^77 = 10~^ 



one finds i] = 0.33. Equation (113) can be also written in terms of the ratio between the transition temperature 
ksT^ and the energy /ma^; one has, in fact, rj — 1.59(fcBT°)^/^(7i^/7na^)~^/^. This expression reveals that in the 
thermodynamic limit, where N 00 and Who with Nujf^^ kept fixed, the parameter 77 has a well defined value. 

In the absence of the condensate (T > Tc) interaction effects are less important because the system is very dilute. In 
this case one can estimate the interaction energy using the expression Ei^t/N ~ gN/R^ where Rt ~ {2kBT /muj'^^Y/'^ 
is the classical radius of the thermal cloud. For temperatures of the order of Tc one finds, 

i^int ^1/6^ _ ^5/2 (^^4) 



This ratio depends on the interaction parameter 77 through a high er power law as compared to the analogous ratio 



for the energy of the condensate, which is linear in 77 [see Eq. ( 113 )], and the effect of E^^t is hence much smaller for 
noncondensed atoms. 



The above discussion emphasizes the importance of the dimensionless parameter (113) which can be used to discuss 



the effe cts o f interactions on the thermodynamic behavior of the system both at low and high temperatures. Actually, 



in Sec. VD we will show that in the thermodynamic limit the system exhibits a scaling behavior on this parameter. 



B. Critical temperature 



The first quantity we discuss is the critical temperature. As anticipated in the previous section, at the onset of 
BEG the system is very dilute and one does not expect atom interactions to give large corrections to thermodynamics. 
Nevertheless the role of interactions on critical phenomena is an important question from a conceptual viewpoint. 
It is interesting to understand, in particular, the differences between the behavior of uniform and nonuniform Bose 
gases. As concerns the comparison with experiments, one should however note that finite size corrections to Tc [see 
Eq. (pO|)] cannot be in general ignored, being in many cases of the same order as the ones due to interactions. 

In the noninteracting model the system can be cooled, remaining in the normal phase, down to the temperature 
T^ which satisfies the condition n{Q)\% = C(3/2) ~ 2.61. Here At = [2T:h^ / {mkBT)Y/'^ is the thermal wavelength 
and n(0) is the density at the center of the trap which, at the critical temperature, is given by the thermal density 
(p^. The presence of repulsive interactions has the effect of expanding the atomic cloud, with a consequent decrease 
of density. The opposite happens for attractive forces, which tend to compress the system. Lowering (increasing) 
the peak density has then the consequence of lowering (increasing) the value of the critical temperature. This effect 
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is absent in the case of a uniform gas where the density is kept fixed in the thermodynamic Hmit. It consequently 
represents a typical feature of trapped Bose gases that is worth discussing in some detail. 

The shift in the critical temperature due to the mechanism described above can be easily estimated by treating 
the interaction in mean-field approximation. The simplest scheme is Hartree-Fock (HF) theory, which consists in 
assuming the atoms to behave as "noninteracting" bosons moving in a self-consistent mean-field; 



H 



HF 



2m 



+ Kfr(r) 

+ V"cxt(r) + 2g?i(r) , 



(115) 



where the last term, 2gn{v) is a mean-field generated by the interactions with the other atoms. This method has 
been first applied to the study of trapped Bose gases by Goldman, Silvera and Leggett (1981) and Huse and Siggia 
(1982) and it has since been adopted in many papers (Bagnato, Pritchard and Kleppner, 1987; Oliva, 1989; Giorgini, 
Pitaevskii and Stringari, 1996, 1997a and 1997b ; Chou, Yang and Yu, 1996; Minguzzi, Conti and Tosi, 1997; Shi and 
Zheng, 1997b). 

In Eq. (115) the quantity n(r) is the total density of the system, the sum of the density of both the condensate and 
thermal components. The single-particle energies and the density n(r) are obtained by solving a Schrodinger equation 
with a density-dependent effective potential. In the presence of Bose-Einstein condensation, the equations for the 
single-particle excitations are coupled to the equation for the order parameter and the whole set of equations must be 
solved using a self-consistent procedure. At zero temperature the Hamiltonian Hhf coincides with the single-particle 
Hamiltonian ( |109| ), which describes the excitations of time dependent Gros s-Pita cvskii theory after neglecting the 
quasi-particle amplitude Vj in the equations of motion (see discussion in Sec. IVE ). 

In the semiclassic al ap proximation (see Sec. II B) one can easily calculate the thermal averages over the eigenstates 
of the Hamiltonian ( 115 ) . The thermal density of the system is given by the ideal gas formula 



nT(r) 



-[Vottir)-IJ.]/kBT 



(116) 



where we have replaced V^xt 
normalization condition 



with [VcS — f-i]- Bose-Einstein condensation starts at the temperature for which the 



N 



J dr nT{r,Tc,fic 



(117) 



can be satisfied with the value of the chemical potential /ic corresponding to the minimal eigenvalue of the Hamiltonian 
(115). For large systems the leading contribution arises from interaction effects 



= 2gn(0) 



(118) 



where, working to the lowest order in g, one can calculate the central density n(0) using the noninteracting model. 
Equation ( |ll8| ) ignores finite size effects, given by ( [lO| ) for the ideal gas. 

By expanding the right hand side of (117) around Hc — and Tc ~ T° one obtains the following result for the shift 
STc = Tc — T^ of the critical temperature (Giorgini, Pitaevskii and Stringari, 1996) 



5T^ 
TP 



Oho 



(119) 



Equation (119) shows that, to lowest order in the coupling constant, the shif t of Tc is linear in the scattering length 
and is negative for rep ulsiv e interactions (a > 0). In this case, the ratio (|ll9| ) can be expressed in terms of the 
parameter 77 defined in (112), and one has STc/T9 — —0.43 ry^/^. For a typical configuration with rj — 0.4, the shift 
is ~ 4%; this can be compared with the shift ( po| ) arising from the finite size correction. Unlike the shift (119) due 
to interactions, the finite size effect (|2^) depends on the anisotropy of the trap and decreases with N. Taking, for 
example, N = 10^ and A = -^8 one finds that finite size effects provide a negativ e cor rection of ~ 2%. For larger 
values of N these corr ections become negligible and one can safely use prediction ( |119| ) . For attractive interactions 
(a < 0), equation (11£) predicts instead a positive shift. However, in this case, finite size effects are always important 
because the value of N cannot be large. 

First measurements of the critical temperature (Ensher et al., 1996), as shown in Fig. |^, indicate the occurrence 
of a negative shift with respect to r° by about 6%, in agreement with the theoretical predictions. However the 
experimental uncertainties are at present too large to draw definitive conclusions from this comparison. 
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Let us conclude this section by recalling that in the mean-field approach discussed above the relation between Tc and 
the critic al de nsity in the center of the trap remains the same as for the noninteracting model: n(0)Ai^^ = C(3/2) = 2.61 
[see Eq. ( |116| )] and it is interesting to look for effects which violate this relation. These can be either finite size or 
many-body effects beyond mean-field theory. These latter effects have been recently calculated in the uniform gas 
through a Path Integral Monte Carlo simulation of the homogeneous hard-sphere Bose gas (Griiter, Ceperley and 
Laloe, 1997). This work has shown that the critical temperature as a function of the gas parameter na^ first 
increases from the noninteracting value T° = (27r?i^/TOfcB)[rt/C(3/2)]^/'^, reaches a maximum for na^ ~ 0.01 where 
Tc/T^ ~ 1.06 and finally decreases for larger values of na^. For the densities relevant for the experiments in traps the 
effects on Tc calculated by these authors are much smaller than the mean-field correction (|119|). 



C. Below Tc 



Below the critical temperature Tc, Bose-Einstein condensation results in a sharp enhancement of the density in 
the centr al region of the trap. This makes interaction effects much more important than above Tc, as discussed in 
Sec. VA. In this section we will be dealing only with systems interacting with repulsive forces and we will consider 
the limit of large N where finite size effects can be ignored. The main purpose is to develop a perturbative scheme 
which permits one to obtain simple analytic formulas for the temperature dependence of the condensate fraction and 
of the energy of the system, providing a useful guide to understand the role of two-body interactions. We will use the 
finite temperature Hartree-Fock scheme already presented in the previous section. An important result emerging from 
this analysis is that, to lowest order in the coupling constant, the corrections to the thermodynamic quantities due 
to interaction effects are linear in the parameter 77 defined in (113). A more complete analysis of the thermodynamic 



behavior, based on self-consistent numerical calculations will be presented in Sec. VP 



A first important problem concerns the temperature dependence of the order parameter and of the chemical po- 
tential. As long as No{T)a/ato ^ 1 and one ignores the interaction with the thermal component, the Thomas- Fermi 
approximation (|5^) to the GP equation provides a good description for the condensate also at T > 0. Equation 
( pT| ) then permits one to estimate the temperature dependence of the chemical potential whose value is fixed by the 
number of atoms in the condensate. One can write 



tiiNo,T) ^l{N,T^O) 



knTO 



knTO 
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N 



2/5 



3^2/5 



(120) 



In order to express the condensate fraction in terms of the reduced temperature t ~ T/T^ we have used the non- 
interacting prediction A^o = A^(l — t^ )- Inclusion of corrections to this law would yield higher order effects in the 



interaction parameter. Equation (120) provi des a useful estimate of /i, which is expected to be accurate in the range 
fi < T < T^. For smaller temperatures, Eq. (120) misses the thermal contributions arising from collective excitations. 
These effects represent however very small corrections and will be ignored in the present discussion. 

Concerning the uncondensed atoms, at high temperature they can be treated as free particles governed by the 
effective mean-field potential T4ff(r) given by (115). The form of this potential can be simplified by ignoring the 
contribution to the density n(r) due to the dilute thermal component and by evaluating the condensate density in the 
Thomas-Fcrmi approximation. This yields the simple result Kff (r) ~ fJ. — \Vcxt{^) — /^| [see also Eq. (115)]. In practice 
most of the thermal atoms occupy regions of space lying outside the condensate where V^xt > fJ- and Vctt — Vcxt ■ As a 
consequence, to a first approximation the effective potential felt by thermal atoms is the same as without interaction. 
However, this does not mean that interaction effects are negligible. In fact, these atoms have a chemical potential 
(120) quite different from the noninteracting value and the corresponding contribution to the thermodynamic averages 



is modified. 

Let us first discuss the problem of thermal depletion. Using the semiclassical picture one can write 



Nt 



drdp 

(27r?i)3 



{exp[(pV2™ + Kff (r) - fi)/kBT] - 1} ' 



(121) 



Explicit integration of (121), using the Thomas- Fermi approximation for the effective mean-field potential, T4ff (r)— /i 
|14xt(r) — fJ-\, leads to the result 



No ^ ^ 
N 



C(3) 



(122) 



valid to the lowest order in the interaction parameter 77. 
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Equation (122) shows that the effects of the interaction depend linearly on rj and are consequently expected to be 
much larger than the ones in the shift of the critical temperature (119) which behave like 77^/^. For example, taking 



rj = 0.4 and t — 0.6 one finds that interactions reduce the value of Nq by 20% as compared to the prediction of 
the noninteracting model. It is worth noting that the quenching of the condensate represents a peculiar behavior of 
trapped Bose gases and takes place because the thermal component of the gas is, in large part, spatially separated 
from the condensate. In a uniform system one has an opposite mechanism. In fact, in this case, the condensate and 
the ther mal c omponents completely overlap and the effective potential is enhanced due to the interaction term 2gn 
[see Eq. ( |ll5| )]. This effect is only partially compensated by the presence of the chemical potential and the final result 
is a suppression of the thermal component. 

In a similar way one can calculate the energy of the system. The main effects of temperature and interactions are 
twofold: on the one hand the number of atoms in the condensate is reduced at finite temperature and the density 
in the central region of the trap decreases. As a consequence the atom cloud becomes larger but more dilute and 
the interaction energy is reduced as compared to the zero temperature case. On the other hand the particles out of 



the condensate are thermally distributed with a modified Bose factor as in (121). By explicitly calculating the two 



contributions, one finds that the total energy of the system exhibits the following temperature dependence: 

E -^^(%i + \(^i^t^)V5^5 + iQt^y (123) 



NkeTO C(3) 7 

Noti ce that the contribution of the interaction, which is again linear in rj, can be obtained directly starting from r esul t 



(120) for the chemical potential, through the use of general relations of thermodynamics. Analogously to Eqs. (|12' 



and (122), expression (|123|) is valid in the temperature regime fi < T < and to the lowest order in the parameter 
V- 

Another useful quantity is the release energy, which corresponds to the energy of the system after switching off the 
trap. Using the same approximations as discussed above, one can easily calculate also this quantity, for which one 
finds the result 

Erel _ E - Eho _ 3C(4) 4 1 , _ 3,2/5 fn , ^'^4 



NkeTO NkeT^ 2C(3) 7 



^4 1^(1 _ i^)^/5 2 + — . (124) 



The release energy can be extracted from time of flight measurements and, consequently, equation (124) provides a 
useful formula to check the effects of two-body interactions directly from experiments. 

The formulas presented in this section account for first order effects in the coupling constant 77. Their validity is 
ensured only for relatively high temperatures and weakly interacting gases. In order to appreciate the accuracy of 



these predictions, in Fig. we compare the energy predicted in (123) with the one obtained by means of a self- 



consistent calculation based on the Popov approximation (see Sec. VE). The agreement is excellent over a wide range 
of temperatures except, of course, very close to where higher orders in 77 give the leading contribution of two-body 
interactions to thermodynamics. 

Expansions similar to the ones discussed in this section can be carried out also in the opposite limit of low tem- 
perature t < rj, which is the analog of the phonon regime of uniform superfluids. Though this regime is not easily 
reachable in current experiments, its theoretical investigation is rather interesting. The low temperature properties 
of trapped Bose gases are deeply influenced by the thermal excitation of the single-particle states localized near the 



surface of the condensate, already discussed in Section IV E (see also, Giorgini, Pitaevskii and Stringari, 1997b) 



D. Thermodynamic limit and scaling 



The thermodynamic limit for the noninteracting gas confined in harmonic traps has been discussed in Sec. II B. 
This limit is reached by letting the total number of particles N increase to infinity and the oscillator frequency Who 
decrease to zero, with the product lj^oN^^'^ kept fixed. Here Who is the geometrical average of the three frequencies. 
This procedure provides a natural extension of the usual thermodynamic limit used in uniform systems where one 
takes N — + 00, V —> 00, and keeps the density n = N/V fixed. In harmonic traps the quantity Lo-hoN^^'^ represents, 
together with T, the relevant thermodynamic parameter of the system and replaces the role played by the density 
in uniform systems. In particular it fixes the value of the critical temperature fc^T^ = ?iti;ho-/V^^^/(C(3))^^'^. In the 
thermodynamic limit all the thermodynamic properties of the noninteracting model can be expressed in terms of 
the critical temperature T° and the reduced temperature t = T/T^. Of course dimensionless quantities, like the 
condensate fraction or the entropy per particle, will depend only on t. 

The thermodynamic limit discussed above applies also in the presence of repulsive interactions. As discussed in 



Sec. V A the parameter rj depends on N and Who only through the transition temperature r° and is consequently 
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well defined in the thermodynamic limit. It is also worth noticing that the dimensionless parameter Na/ai^o, which 
characterizes the effects of two-body interactions in the Gross-Pitaevskii equation for the ground state behaves as 
Na/a\io A^^/^ry^/^. Thus, in the thermodynamic limit, the condition Noa/a^o ^ 1 which ensures the validity of the 
Thomas- Fermi approximation for the condensate, is always guaranteed below Tc- 

The above discussion suggests that in the thermodynamic limit the relevant functions of the system will depend 
on T^, t and 77 (Giorgini, Pitaevskii and Stringari, 1997a). In the previous section we have already anticipated such 



a be h avior by calculating some relevant thermodynamic functions to the lowest order in 77, as done in Eqs. (120) 



(122)-(124). This points out a scaling behavior exhibited by these systems. Different configurations, corresponding to 
different values of A^, m, trapping frequencies and scattering length, will be characterized by the same thermodynamic 
behavior provided they correspond to the same value of ry. 

The scaling behavior can be proved in a general way by noting that in the limit cjjjo — > the size of the system 
increases and the density n(r) changes very slowly. As a consequence in the thermodynamic limit the density is fixed 
by the condition ( [5^ ) of local equilibrium, //(T) — fJ-ioca.i{n,T) + VIsxtCr), where ^i\ac&\{n,T) is the chemical potential 
of an interacting uniform system at density n = n(r). By inverting the above condition one can write the density 
of the gas in the form n(r) = — Vc^t{'r),T), where n{^,T) is the density of the uniform gas as a function of 
chemical potential and temperature. Notice that the inversion of the function /iiocai(f^) requires that fiiocai be a 
monotonous function of the density. This condition is satified by interacting systems where the stability condition 
implies d^/dn > 0. The total number of atoms N is obtained by integrating the density over space co-ordinates. 
Introducing the new variable ^ — Ksxt(r), one can write 
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The parameters of the trap enter the above equation through the combination Who-^^^^ T^- On the other hand the 
integral on the right hand side requires the knowledge of the density of the uniform system as a function of chemical 



potential and T . Equation (125) shows that the knowledge of the thermodynamic behavior of the interacting uniform 
system would permit one, in principle, to determine the thermodynamics of the trapped gas (Damle et al., 1996). 
For a dil ute Bose gas, where the interaction is accounted for by a single parameter (the scattering length a), the 



integral (125) depends on the quantities /i, T and jmcP' ^ the latter bein g the only energy one can construct with the 



mass m and the scattering length a. As a consequence, inversion of ( 125 ) yields the following general dependence for 
the chemical potential /i ~ fi{T,T^, h^/nia^). Due to dimensionality arguments the above expression can be always 
written in the form 

lJi = kBT^J{t,7^) . (126) 

Here / is a dimensionless function depending on the reduced temperatu re t = T/T^ and on the scaling parameter 77, 
fixed by the ratio between ksT^ and jma^ [see discussion after Eq. (113)]. A similar scaling behavior applies to 
the other thermodynamic functions. 

The above discussion applies to the thermodynamic limit where N 00. An important question is to understand 
whether in the available experimental conditions, where N ranges between 10'* and 10^, this limit is reached in practice 
or finite size effects are still significant. In Fig. ^ we show, as an example, the behavior of the condensate fraction. 
This quantity depends, in the thermodynamic limit, only on the variables t and rj. In the figure we plot the numerical 



results obtained from a self-consistent mean- field calculation based on the Popov approximation (see Sec. VE) for 
two very different configurations, both corresponding to the same value of 77 = 0.4. Open squares correspond to 
= 5 X 10"' rubidium atoms in a trap with a/oho = 5.4 x 10^^ and A = VS, while solid squares correspond to 
= 5 X 10^ sodium atoms in a trap with a/aho — 1-7 x 10^'^ and A = 0.05. One sees that both set of data coincide 
with the asymptotic scaling function (solid line), calculated with the same value of rj, by taking the thermodynamic 
limit in the equations of the Popov approximation (Giorgini, Pitaevskii and Stringari, 1997b). The figure points 
out, in an explicit way, how very different configurations can give rise to the same thermodynamic behavior, if the 
corresponding scaling parameter 77 is the same. It is also interesting to notice that the scaling behavior is reached 
faster in the presence of two-body interactions than for noninteracting particles. In the latter case, in fact, finite size 
effects, which are responsible for the deviations from the scaling law (1 — t^) (dashed line), are more visible (open and 
solid circles). 

The scaling behavior is very well reproduced also by the other thermodynamic quantities. For this reason, in the 
next section, we will discuss the behavior of interacting Bose gases confined in harmonic traps calculating directly the 
various physical quantities in the thermodynamic limit. 
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E. Results for the thermodynamic functions 



In Sec. V C we have used Hartrcc-Fock theory to estmiate the temperature dependence of the chemical potential, 
condensate fraction, release energy to lowest order in the scaling parameter 77. The full equations of Hartree-Fock 
theory can be solved numerically in a self-consistent way (Minguzzi, Conti and Tosi, 1997; Shi and Zheng 1997b; 
Giorgini, Pitaevskii and Stringari, 1997b) going beyond the perturbative scheme. Hartree-Fock theory is expected to 
be quite accurate at high temperature except, of course, very close to Tc, where mean-field theories are inadequate 
(Shi, 1997; Shi and Griffin, 1998). The accuracy of Hartre e-Foc k at high T is justified by the crucial role played by 



single-particle excitations, as we have already seen in Sec. 
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for the density of states. This theory is instead less 



accurate at low T, since it ignores the effects of collectivity which characterize the low energy part of the excitation 
spectru m. Su ch collective effects are instead properly included in time dependent Gross-Pitaevskii theory, as discussed 
in Sec. IV A. A mean-field scheme, describing correctly both the high and low temperature regimes, is provided by 



the so called Popov approximation (Popov, 1965 and 1987; Griffin, 1996; Shi, 1997) whose application to interacting 
bosons in harmonic traps has been considered by several authors in the last few years. This mean-field scheme is 
based, on the one hand, on a finite T extension of the Gross-Pitaevskii equation, in which the interaction between 
condensed and noncondensed atoms is explicitly accounted for, and, on the other hand, on Bogoliubov-type equations 
for the excitations of the system. The corresponding equations have the form 
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where rto(r) = |0(r)P is the condensate density, while the ther mal d e nsity n-r is calculated through the relation 
ut = + bl^)[6xp(/3£j) — 1]~^ with Uj, Vj and Sj solutions of (^28|)-(|l29|). These quantities are now temperature 

dependent. The sum n{r) — no(r) -I- nT(r) is the total density. The functions Ui, Vi entering Eqs. (128)-(129) are 
normalized according to condition (|69|). Notice that in this approximation the thermal component is treated as a 
thermal bath generating an additional static external field in the equation for the condensate. One also ignores here 
the T ~ quantum depletion nout(r) = which has been shown to be very small in these trapped gases 



(see Sec. [HE). Finally, at low temperature, whe n ut i s negligible compared with uq, equation (127) coincides with 
the stationary GP equation (|39|) while Eqs. (|l28|)-(|l29|) reduce to (|6|)-(|3). 

From the solution of Eqs. (127)-( p^9| ) one obtains density profiles in good agreement with experimental data. As 
an example, we take two of the density profiles already shown in Fig. |2^ and we plot them again in Fig. ^ together 
with the theoretical prediction from Eqs. (127)-( p^ ), using the number of particles, N, and the temperature, T, as 
fitting parameters. The same equations have been used also for fitting the experimental data by Hau et al. (1998). 

Using the distribution function of the excited states, fj = [exp(/3ej) — 1]~^, and the combinatorial expression for 
the entropy, S — ks'^jiP^jfj — ln[l — exp(— /3ej)]}, one can work out all the thermodynamic quantities (Giorgini, 
Pitaevskii and Stringari, 1997b). The comparison between the predictions of Hartree-Fock and Popov theories has 
revealed that there are no significant differences between the two approaches for most thermodynamic quantities. 
Only increasing the value of the interaction parameter 77 one can observe some differences. This is a further evidence 
of the negligible role played in thermodynamics by the collective modes of the condensate, which are ignored in the 
Hartree-Fock scheme. This behavior, in accordance with the analysis of the density of states made in Sec 
represents a significant difference with respect to the case of uniform superfiuids. 

In the following we present results for various thermodynamic quantities (chemical potential, condensate fraction 
and release energy) obtained using the Popov approximation in the thermodynamic limit. As discussed in the 
preceeding section, this limit is well achieved in the configurations realized in present experiments. The results are 
consequently presented as a function of the reduced temperature t for different values of the scaling parameter rj. We 
have considered the value rj = 0.4, which corresponds to the typical configurations realized in actual experiments, and 
the value r/ = 0.6, which would correspond to a more correlated gas. 

First, in Fig. Pq, we show the chemical potential in units of fcsT^ as a function of the reduced temperature t. 
Notice that for t — > the plotted quantity coincides, by definition, with the parameter 77 [see definition (112)]. In the 
classical limit, T ^ T^, the chemical potential approaches the ideal gas value ^/UbT^ = t\n{C,{3)/t^). 
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The results for the condensate fraction Nq/N are given in Fig. The open circles are the experimental points 
(Ensher et at, 1996). In the experiment the number of atoms N varies with t and the corresponding value of i] ranges 
from 0.39 to 0.45. The data are compared with the predictions of the mean-field theory for rj = 0.4 (solid line) and 
the noninteracting gas model (dotted line). The experimental points are shifted from the noninteracting case to lower 
temperature, but not as much as predicted by mean-field theory. However, the experimental uncertainties are still 
too large to draw any definitive conclusion. 

In Fig. we show the results for the release energy. The dots are the experimental data (Ensher et al, 1996) 
which, below Tc, lie well above the noninteracting curve, showing again a clear evidence for the effects of two-body 
interactions. We notice that, by differentiating the total energy with respect to T, one could calculate the specific heat. 
In the present experiments however, this quantity is not directly available because only the release energy is measured. 
What one can see from a fit to the experimental data of the release energy is the occurrence of a characteristic bump 
in the derivative near the transition temperature (Ensher et al., 1996); this behavior is in good agreement with the 
prediction obtained by taking the derivative of the theoretical curves in Fig. 

Finally, the above results of the mean-field theory at finite temperature can be also compared with the ones of 
Quantum Monte Carlo calculations (Krauth, 1996; Holzmann, Krauth and Naraschewski, 1998). It is worth noticing 
that the possibility of making a close comparison between exact Monte Carlo simulations, experimental data and 
mean-field calculations is a rather rare event in the context of interacting many-body systems and represents a further 
nice feature of BEC in traps. In Fig. the condensate fraction obtained with path integral Monte Carlo calculations 
by Krauth (1996) is represented by solid circles with error bars. The simulation has been done with 10000 atoms 
interacting through an hard-core potential, and corresponds to -q = 0.35. The results are very close to the mean-field 
prediction. The value of rj used for the solid curve is actually 0.4, since this value is closer to the experimental 
situation, but the same calculation for rj = 0.35 gives an even better agreement, crossing precisely the three Monte 
Carlo data at high temperature. A detailed comparison between Monte Carlo results and mean-field theory has been 
recently performed by Holzmann, Krauth and Naraschewski (1998), including the analysis of the density profiles of 
the gas at different temperatures. 



F. Collective modes at finite temperature 

In Sec. we have studied the collective excitations of a trapped Bose gas at zero temperature. In this case, all the 
atoms are in the condensate and there are no collisions. In the collisionless regime the force acting on a given particle 
comes from the mean-field created by the other particles; this field generates a collective oscillation of the system, 
which is sometimes called Bogoliubov sound and is the analog of zero sound for normal Fermi liquids. 

At finite temperature the situation is more complicated. On the one hand both the condensate and the thermal 
cloud can oscillate. On the other hand collisions between excitations can play an important role and one must 
distinguish between a coUisional and a collisionless regime. 

So far the temperature dependence of these oscillations has been analysed experimentally both at JILA and MIT. In 
the first case, Jin et al. (1997) investigated the m = and m = 2 modes for a system of ^ 10"* Rb atoms in a trap with 

asymmetry parameter A = -\/8. At low temperatures (T ~ 0.4r°) they found that the oscillations of the condensate 



have frequencies in good agreement with the predictions of the T = Gross-Pitaevskii equation (see Sec. [IV| ), while 
at higher temperatures the frequency exhibits an unexpected temperature dependence with very different behavior 
for the two modes. In the second case, Stamper-Kurn et al. (1998c) studied the low-lying m = mode for a much 
larger system (TV ^ 10^ Na atoms) in a cigar-shaped trap with A <C 1. Similarly to the JILA group, they observed 
a shift of the collective frequency with respect to the T = value. Both groups measured also the damping of the 
collective modes of the condensate, finding a rather strong T-dependence. They also observed the oscillations of the 
thermal cloud below and above the critical temperature. 

Mean-field approaches have been used to predict the properties of the collective excitations in the collisionless regime. 
This regime is achieved at low temperatures and for low densities of the thermal cloud. In this case, the oscillations 
of the condensate behave similarly to the T = case and can still be called Bogoliubov's sound modes. Several 
authors have used finite T extensions of the time-dependent Gross-Pitaevskii equation in the Popov approximation 



[see Eqs. (|127|)-(|29|)] (Hutchinson, Zaremba and Griffin, 1997; Dodd et al, 1998a and 1998b). This approach, in 
which the thermal component is treated as a static thermal bath, does not account for any damping mechanism. In 
order to include damping, a dynamic description of the oscillations of both the condensate and the thermal cloud is 
needed (Minguzzi and Tosi, 1997; Giorgini, 1998). The dynamic coupling between the motion of the two components 
might also be important in the determination of the temperature dependence of the frequency shift. There are still 
open questions on this problem and several aspects of the theory are expected to be clarified in the next future. 
Damping processes in the collisionless regime have been investigated using perturbative approaches (Liu and Shieve, 
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1997; Pitaevskii and Stringari, 1997; Liu, 1997; Fedichev, Shlyapnikov and Walraven, 1998). An important mechanism 
of coUisionless damping is provided by the Landau damping. In Fermi Hquids, Landau damping originates from the 
couphng between single-particle excitations and zero sound at T = (Lifshitz and Pitaevskii, 1981, §30; Pines 
and Nozieres, 1966, Vol.1). In a Bose gas, an analogous damping occurs because the thermal bath of elementary 
excitations can absorb quanta of the collective oscillation. Landau damping increases with temperature, because 
of the larger number of elementary excitations available at thermal equilibrium. For temperatures larger than the 
chemical potential it increases linearly with T and for a uniform system the ratio between the imaginary, 7 — — Im(ct;), 
and real, Re(aj), parts of the collective frequency takes the analytic form (Szepfalusy and Kondor, 1974; Shi, 1997; 
Pitaevskii and Stringari, 1997) 

7 _^^ksTa ^^3^^ 
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This equation is valid for excitations with energ y KR b{uj) <^ 1^1 <ti ksT, which in a uniform system are phonons 
propagating with the sound velocity c. Equation ( |130| ) can be used for a rough estimate of the damping of collective 
excitations in traps by taking for c the value of the sound velocity in the center of the trap. As an example, in Fig. ^ 
we show the dampin g ra te, 7, measured for the m = and m — 2 modes by Jin et al. (1997), compared with the 
theoretical estimate ( |130| ). Taking into account that this estimate, which is expected to apply at high temperatures, 
is very rough and that the effects of geometry and multipolarity are completely neglected, the agreement between 
theory and experiment can be considered reasonable. Recently Fedichev, Shlyapnikov and Walraven (1998) have 
indeed argued that Landau damping is strongly influenced by the geometry of the traps, being particularly effective 
in anisotropic traps, due to the randomness of the excitation spectrum. 

At high temperature and/or high density, collisions are more important and can affect the nature of collective 
excitations. In Bose superfluids the coUisional regime is described by the equations of two-fluid hydrodynamics and 
is characterized by the occurrence of two distinct oscillations: first and second sound. In liquid ^He, first sound is 
a density wave with in phase oscillations of the superfluid and normal components, while second sound is an almost 
pure entropy wave with opposite phase motion of the two components. For a dilute Bose gas the situation is different 
because the interaction between the condensate and the thermal cloud is very weak. In particular, except at very 
low temperature, first sound mainly involves the thermal cloud and reduces to the usual hydrodynamic sound above 
Tc] conversely, second sound is essentially the oscillation of the condensate and disappears above Tc (Lee and Yang, 
1959; Griffin and Zaremba, 1997). Similar features have been pointed out also in the presence of harmonic trapping 
(Zaremba, Griffin and Nikuni, 1997; Shenoy and Ho, 1998). 

An important limiting case is represented by the collective motion of the gas above the critical temperature, 
where the gas exhibits an almost classical behavior. In the coUisional regime, one can then use the equations of 
hydrodynamics for classical gases in order to obtain explicit results for the collective frequencies. For example the 
coupled quadrupole and monopole modes with to = obey the following dispersion law (Griffin, Wu and Stringari, 
1997; Kagan, Surkov and Shlyapnikov, 1997a): 

1, 



oj^ ^- [5luI + AljI ± (25^1 + 16uj^ - 32^1^2)1/2] _ (131) 
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If the trap is spherical {uj± = lUz — Who), tbe two solutions have frequency ui — \/2 tJho (quadrupole) and uj — 2uJho 
(monopole). Notice that the frequency of the surface (quadrupole) mode is equal to the one of the oscillation of the 
condensate at zero temperature, given by Eq. (p3|). In the limit of highly deformed cigar-shaped traps (ujz ^ 
the lowest frequency of (|l3lD becomes tu = ^/Yljh 



This coUisional regime, yielding the dispersion law (131), is achieved if ojt <C 1, where t is a typical collision time. 



In the opposite limit (wr 3> 1) one finds the coUisionless regime, where the gas oscillates with frequencies fixed by 
the trapping potential, corresponding to the predictions of the noninteracting model. Oscillations of this type have 
been observed by Jin et al. (1997) for both the gas above Tc and the thermal component below Tc\ the frequency of 
the TO = and to = 2 modes was found to be roughly twice the trap frequency, in agreement with the prediction for 
the ideal gas. 

In the intermediate regime where lot ^ 1 one expects a smooth cross-over from coUisionless to coUisional hydrody- 
namic modes. A useful interpolation formula is provided by the law (Kavoulakis, Pethick and Smith, 1998) 

^1 MiD^:^ ^ (132) 
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predicted by general theory of relaxation phenomena (Landau and Lifshitz, 1987, §81). Here loc and w/^d are the 
frequencies of the mode in the coUisionless and coUisional hydrodynamic regimes respectively. 
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An estimate of the collisional time r ca n be obtained by considering again a classical picture of the system. One 
finds r ~ Imfp/vT, where vt = y^2kBT/m is the thermal velocity and ^mfp = {nrcr)'^ is the mean free path which is 
fixed by the s-wave cross section a — Stto^ and by the density. For frequencies of the order of the oscillator frequency 
the condition lut <C 1 is equivalent to requiring that the mean free path be much smaller than the thermal radius 
Rt = y/2kBT/mLj^^. Near Tc o ne finds that the collisional frequency l/r behaves like rj^kBTc/h where 77 is the 
scaling parameter defined in (|113| ). Despite the smallness of the factor 77^, the coUisionless frequency can be of the 
same order of the collective frequencies of the system because of the factor N^^^ contained in so that increasing 
N favours the achievement of the collisional hydrodynamic regime. Notice that the multipolarity of the excitation 
can play an important role in characterizing the relaxation of the collective oscillation. For example the dipole mode 
cannot have any relaxation mechanism in the presence of harmonic trapping. The same happens in the case of the 
monopole excitation if the harmonic trap is isotropic. In both cases the coUisionless (coc) and hydrodynamic {lohd) 
frequencies coincide. 

In Fig. ^ we plot the imaginary part of uj against the real part, as given in ( |132| ), for the case of the low- lying 
m — mode observed in the cigar-shaped trap at MIT (Stamper-Kurn et ai, 1998c). In this experiment, the thermal 
cloud is found to oscillate with a frequency of 1.78^^, which is larger than the hydrodynamic prediction ■\/l2/5 lOz, 
but lower than the noninteracting value 2uJz- A damping rate of about 20 s~^ was also observed, corresponding to 
— Im(tj) — 0.190;^. In Fig. 30 the experimental results are represented by the solid circle with error bars, which turns 
out to be reasonably close to the theoretical curve, the difference being of the order of the experimental uncertainty. 
The part of the curve near the maximum corresponds to values of collision time such that Re(cj)r ^ 1 and this 
suggests that, differently from the JILA experiment (Jin et ai, 1997), the motion of the thermal cloud in the MIT 
experiment is affected by collisions. 

Finally, we mention that Stamper-Kurn et al. (1998c) observed also the opposite phase dipolar oscillation of the 
thermal cloud and the condensate, occurring below Tc (Zaremba, Griffin and Nikuni, 1998). This mode exhibits 
strong damping. 



VI. SUPERFLUIDITY AND COHERENCE PHENOMENA 



Superfluidity is one of the most spectacular consequences of Bose-Einstein condensation. However, the explicit 
connection between superfluidity and BEC is not trivial and has been the object of a longstanding and deep investi- 
gation in the last decades, mainly for its importance in understanding the physics of liquid helium. In macroscopic 
bodies superfluidity shows up with many peculiar features: absence of viscosity, reduction of the moment of inertia, 
occurrence of persistent currents, new collective phenomena (second sound, third sound, etc.), quantized vortices, 
and others. Several properties are usually interpreted as coherence effects associated with the phase, S, of the or- 
der parameter whose gradient fixes the velocity of the superfluid through Vj — (fi/m)'^ S. A major question is to 
understand whether some of these effects can be observed also in trapped gases. Of course in a mesoscopic system 
one expects the manifestations of superfluidity to be different from the ones exhibited by macroscopic bodies. In par- 
ticular, traditional experiments based on the study of transport phenomena are not easily feasible in trapped gases. 
On the other hand, interference patterns, associated with phase coherence, have been already observed (Andrews, 
Townsend et ai, 1997) and successfully compared with theory. This opens a promising field of research based on the 
investigation of coherence phenomena, including the realization of the so-called "atom laser" (Ketterle, 1998). 



A. Rotational properties: vortices and moment of inertia 



Among the several properties exhibited by superfluids, the occurrence of quantized vortices and the strong reduction 
of the moment of inertia represent effects of primary importance. 

In a dilute Bose gas the structure of quantized vortices can be investigated starting from the Gross-Pitaevskii 
equation. Indeed one of the primary motivations of the GP theory was the study of vortex states in weakly interacting 
bosons (Gross, 1961 and 1963; Pitaevskii, 1961). These studies were further developed by Fetter (1972) including 
higher order effects in the interaction. 

A quantized vortex along the z-axis can be described by writing the order parameter in the form 

0(r) = 0i,(rj_, z) exp[iK(y9] (133) 

where ip is the angle around the z-axis, k is an integer, and (/)„(r^,z) = n{r±, z). This vortex state has tangential 
velocity 
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The number k is the quantum of circulation and the angular momentum along z is NkH. The equation for the 
modulus of the order parameter is obtained from the GP equation (|39|). The kinetic energy brings a new centrifugal 
term arising from the velocity flow which pushes the atoms away from the z-axis. The GP equation then takes the 
form 



2m 



4'v{r±,z) = fi(f>v{r±, z) 



(135) 



Due to the presence of the centrifugal term, the solution of this equation for k 7^ has to vanish on the z-axis. An 
example is shown in Fig. where the solid line represents the condensate wave function, cf)v{x^ 0, 0), for a gas of 10^ 
rubidium atoms in a spherical trap and with vorticity k = 1. In the inset, we give the contour plot for the density in 
the xz-plane, n{x,0,z) = |0„(a;, 0, z)p. 

For noninteracting systems the solution of Eq. ( 135 ) is analytic and, for k = 1, has the form 
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In this case, the vortex state corresponds to putting all the atoms in the m = 1 single-particle state. Its energy is 
then NTiLU± plus the ground state energy. In Fig. gl|, the corresponding wave function is shown as a dashed line. 
Similarly to what happens for the ground state without vortices, the presence of repulsive forces reduces dramatically 
the density with respect to the noninteracting gas, the condensate wave function becoming much broader. 

The structure of the core of the vortex is fixed by the balance between the kinetic energy and the two-body 
interaction term. For a uniform Bose gas the size of the core is of the order of the healing length ^ = (STrna)"-'^/^, 
already introduced in Sec. III.B, where n is the density of the system. For the trapped gas, an estimate of the core 
size can be obtained using for n the central value of the density in the absence of vortices. If the trap is spherical, 
as in the case of Fig. the ratio between ^ and the radius R of the condensate takes the form (Baym and Pethick, 
1996) 
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where we have used the Thomas-Fermi approximation for the central density and the radius R. For the condensate 
in the figure, the radius is about 4.1 in imits of Oho and the ratio ^/R is then ~ 0.06. The actual core size depends, 
obviously, also on the position on the z-axis and becomes larger when the vortex line reaches the outer part of the 
condensate, where the density decreases. This can be clearly seen in the inset of Fig. |l]. 

The energy of the vortex can be evaluated through the energy functional (p7|). The difference between the energy 
of the vortex state and the one of the ground state allows one to calculate the critical frequency needed to create a 
vortex. In fact, in a frame rotating with angular frequency 51, the energy of a system carrying angular momentum Lz 
is given by {E — flL^), where E and are defined in the laboratory frame. At low rotational frequencies this energy 
is minimal without the vortex. If f2 is large enough the creation of a vortex can become favourable due to the term 
—flLz- This happens at the critical frequency 
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(138) 



where E'k is the energy of the system in the presence of a vortex with angular momentum Nfin. In Fig. |3^ we plot 
the critical frequency for the creation of a vortex with k = 1 as a function of the number of atoms, for rubidium in a 
spherical trap. As shown in the figure, the predicted critical frequency is a fraction of the oscillator frequency and, in 
typical experimental conditions, corresponds to a few Hz. It decreases when N increases, because for large systems 
the energy cost associated with the occurrence of a vortex increases as InA^, while the gain in VlLz is always linear 
in N . This behavior is similar to the one exhibited by uniform systems where, approximating the vortex core as a 
cylindrical hole of radius ^, one can calculate explicitly the critical frequency; one finds VLc = {fi/mR^) ln(i?/^), where 
R is the radius of the region occupied by the vortex flow. Analogous expressions can be derived also in the presence 
of harmonic trapping for large N . Baym and Pethick (1996) and Sinha (1997) have shown that the critical frequency, 
in units of Who, goes as ~ (aho/-R)^ ln(-R/^). Using the asymptotic solution of the GP equation in the large N limit, 
Lundh, Pethick and Smith (1997) have found a useful analytic expression for the critical velocity in the case of axially 
symmetric traps: 
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where R± is the Thomas-Fermi radius of the cloud in the xy-plane, orthogonal to the vortex line, while the healing 
length is defined by ^ = (87rna)~^/^, with n equal to the central density of the gas without vortex. This formula gives 
a critical frequency which differs significantly from the numerical result shown in Fig. ^ only for N smaller than 
about 2000, while for larger N it becomes more and more accurate. 

The above discussion regards the structure of vortices for repulsive interactions. An intriguing problem is how the 
angular momentum is distributed in systems with attractive forces. This question was recently addressed by Wilkin, 
Gunn and Smith (1998), who showed that the lowest eigenstates of fixed angular momentum do not exhibit vortex 
configurations. 

A major question concerning vortices in trapped Bose gases is whether they can be observed in experiments. So 
far no evidence about their existence has been reported. In principle, it should not be difficult to produce them in 
a steadily rotating trap. The value of the critical frequency is in fact easy to achieve in the laboratory. However, 
when one stops the rotation it is not obvious whether the vortex remains stable. The problem of stability of vortices 
is rather complex even in uniform superfluids, like ^He, where it has been the object of much experimental and 
theoretical work (Donnelly, 1991). In a recent paper, Rokhsar (1997) has argued that a vortex placed at the center of 
a nonrotating harmonic trap is unstable. Other discussions about the stability of vortex configurations can be found 
in Fetter (1998), Benakh et al. (1997), Isoshima and Machida (1998). 

Creating a vortex is only part of the problem. A second important problem is its detection. The excitation energy 
associated with a vortex is too small to be observed with measurements of the release energy. In fact the increase in 
the energy per particle is hflf,, a quantity much smaller than the energy per particle in the ground state, given by 
{5/7) ji. However, imaging the core of the vortex during the expansion, after switching-off the trap, should be feasible, 
as recently suggested by Lundh, Pethick and Smith (1998). Promising perspectives are also given by the effects of 
vortices on the shift of the collective frequencies of the condensate. These can be measured with high precision and 
the observation of a breaking of degeneracy between states of opposite angular momentum would represent a rather 
unambiguous evidence of the presence and the quantization of the vortex. The shift of the collective frequencies 
has been already investigated by several authors. Sinha (1997) used a large N semiclassical expansion in the Gross- 
Pitaevskii equation, Dodd et al. (1997a) carried out a direct numerical solution of the same equation. Very recently 
Zambelli and Stringari (1998) have developed a sum rule approach yielding an explicit analytic expression for the 
splitting of the quadrupole modes, while Svidzinsky and Fetter (1998) have developed a perturbative solution of the 
collisionless hydrodynamic equations. One should also recall that the existence of vortices gives rise to a new series of 
collective excitations localized near the vortex core. Similar modes are found in uniform superfluids. In that case, the 
corresponding dispersion law of the lowest mode is a; = (?ifc^/2m) ln(fc2^), where is the wavevector associated with 
a periodic motion of the vortex line along the z-axis. The effects of thermally excited vortex waves in trapped gases 
have been explored by Barenghi (1996). Another explicit proof of the existence of vortices would be the observation 
of an asymmetry in the velocity of sound when one considers wave packets propagating in the same, or in the opposite 
direction with respect to the vortex fiow. Such a test on the quantization of the superfiuid fiow requires a ring-type 
geometry for the confinement of the atomic cloud. The quantization of the superfiuid flow in the ring could be also 
revealed using interference experiments (see next section) . By letting the condensate expand one should in fact observe 
interference fringes associated with the modulations of the phase produced by the quantization of the circulation in 
the ring. In a similar way, one could detect a quantized vortex by looking at the phase slip in the interference fringes 
produced by two expanding condensates (Bolda and Walls, 1998). 

We have already pointed out that if we induce at zero temperature rotations on an axially symmetric system with 
angular frequency smaller than f2c, then the system remains in its ground state. In fact only the normal (nonsuperfluid) 
part of the system can participate in the rotational motion and, consequently, axially symmetric Bose systems can 
possess a moment of inertia only at flnite temperature. A deviation of the moment of inertia from the rigid value 
represents an important manifestation of superfluidity. In liquid *He such deviation has been directly observed below 
the lambda temperature, where the system becomes superfiuid (Donnelly, 1991). Measuring the moment of inertia of 
a trapped gas is actually a challenging problem, because direct measurements of angular momentum are difficult to 
obtain. 

The moment of inertia O relative to the z axis can be defined as the linear response of the system to a rotational 
field Hcxt = ~^Lz, according to the definition 

(L,) = c^e , (140) 

where = Tlii^iPi ~ UiPi) is the z component of the angular momentum and the average is taken on the state 
perturbed by -ffoxt • For a classical system the moment of inertia takes the rigid value 

e,ig = mN{rl) . (141) 

Vice versa, the quantum mechanical determination of O is much less trivial. It involves a dynamic calculation and, 
according to perturbation theory, can be written as 
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2^\{j\L^ ( 



e-^l^^f^exp --^), (142) 



2J 

where |j) and |i) are eigenstates of the unperturbed Hamiltonian, Ej and Ei are the corresponding eigenvalues, and 
Z is the partition function. 

The moment of inertia can be easily calculated if one considers the simplest case of an ideal gas trapped by a 
harmonic potential. The result is (Stringari, 1996a) 

e = elm{rl)QNQ{T) + m{rl)TNT{T) , (143) 

where the indices ()o and mean average taken over the densities of the Bose condensed and noncondensed com- 
ponents of the system, respectively. The quantity 

(x^ + r )o 

is the deformation parameter of the condensate given by {ujy — ujx)/{^y + uJx)- This quantity vanishes for an axially 
symmetric trap. 



The physics contained in (143) is very clear. In fact, the first term in the moment of inertia arises from the 
atoms in the condensate, which contribute with their irrotational flow and can be hence interpreted as the superfluid 
component. The second term arises instead from the particles out of the condensate which rotate in a rigid way 
(normal component). These two distinct contributions are at the origin of an int erest ing T-dependence of Q. In fact. 



above Tc, where A'o = 0, the moment of inertia takes the classical rigid value (141), while at T = 0, where all the 
atoms are in the condensate, it is given by the irrotational value Qmot — eo®iig- the limit of small deformation 
(i.e., small eo)j the deviation of the moment of inertia from its rigid-body value is given by the useful expression: 



^'^ (145) 



Orig No{rl)o + NT{rl)T ' 

It is worth discussing the behavior of the moment of inertia in the thermodynamic limit iV — s- cx), with Nujf^^ kept 
constant. In this limit, the ratio (?'i)o/(?'i)r tends to zero. In fact the square radius of the condensate increases as 
l/uj±, while the one of the thermal cloud as fcsT/w^. In this limit, one then finds 0/8rig 1 everywhere except at 
T ~ where Nx — 0. This behavior is not surprising. In fact if the radius of the condensate is much smaller than 
the one of the thermal component then there is no distinction between Q and Oiig since in both cases the leading 
contribution is given by the thermal component. For finite values of N the ratio Q/Oiig goes smoothly to zero as 
temperature decreases. An example is shown in Fig. ^ for a spherical trap and two different values of N (dashed and 
dot-dashed hnes). 

How do two-body interactions change the above picture? Result (145) is expected to be valid also in the presence 
of interactions, to the extent that the relevant excitations are well described by a single-particle picture and the 
condensate can be still identified with the superfluid component. For example in Hartree-Fock theory one obtains 
coupled equations for the condensate wave function and the single-particle excited states. The irrotational flow for the 
condensate follows from the definition of the superfluid velocity as the gradient of the phase of the order parameter. 
On the other hand, the flow of the atoms out of the condensate is rigid-like if one treats the thermally excited states 



in semiclassical approximation, as done in Sec. V [see Eq. (121)]. Only at very low temperature, where the effects of 



collectivity can be important and the superfluid component must be distinguished from the condensate, expression 



(145) for the moment of inertia is no longer correct. 

Interactions can affect the value of the moment of inertia by changing the temperature dependence of the condensate 
as well as the value of the square radii (r^). The change in the radii is particularly significant at large N. In fact, 
unlike for the noninteracting case, the ratio (''j^)o/(''1)t does not vanish in the thermodynamic limit and is fixed by 
the value of the scaling parameter rj. As a consequence interactions have the important effect of reducing the value 
of the moment of inertia with respect to the rigid value in the whole range of temperatures below Tc- This behavior 
is explicitly shown by the solid line in Fig. |3[ 

B. Interference and Josephson effect 

An important consequence of phase coherence in Bose-Einstein condensates is the occurrence of interference phe- 
nomena. A beautiful example is the experiment carried out at MIT (Andrews, Townsend et ai, 1997), where a laser 
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beam was used to cut a cigar-shaped atomic cloud into two spatially separated parts. After switching off the confining 
potential and the laser, the two independent atomic clouds expand and eventually overlap. Clean interference patterns 
have been observed in the overlapping region (see Fig. ^Jb) . 

From a qualitative viewpoint, one can imagine the initial condensates as two point-like pulsed sources placed at 
distance d on the z-axis. Let us consider the interference taking place in the region of space where the density of 
the gas is small enough and hence the condensate wave function is a linear superposition of two de Broglie waves. 



By using result (100) for the phase of each expanding condensate, one finds that the relative phase of the two waves 
behaves as [S{x, y,z + d/2) — S{x, y,z — d/2)] ~ {m/K)az{t)zd. If the time delay, t, between the switching off of the 
trap and observation is large, one has az{t) — > 1/t. One then finds straight interference fringes, orthogonal to the 
z-axis, with wavelength given by 

Using the spacing between the two initial condensates as an estimate of the distance d one gets fringe periods in 
reasonable agreement with the observed patterns (Andrews, Townsend et ai, 1997). Typical values are t ~ 40 ms, 
d ~ 40 /im, and A ~ 20/im. 

Gross-Pitaevskii theory is a natural framework for investigating interference phenomena in a quantitative way. In 
this theory, the phase coherence of the condensate is assumed from the very beginning. The interference patterns can 
be obtained by solving numerically the GP equation (|3^) for two condensates. This has been done for instance by 
Hoston and You (1996), Naraschewski et al. (1996), Rohrl et al. (1997), WaUis et al. (1997a and 1997b). Interference 
phenomena have been investigated also without using the concept of broken gauge symmetry by Javanainen and 
Yoo (1996). In Fig. ^ the experimental results are plotted together with the theoretical calculations by Rohrl et al. 
(1997). The good agreement between theory and experiment reveals that the concept of phase coherence, as assumed 
in GP theory, works very well. This was not obvious a priori, since the system is finite-sized and interacting and 
hence phase coherence is expected to be only approximate. 

Another interesting manifestation of phase coherence in trapped condensates is the possible occurrence of Josephson- 
type effects, in analogy with well known properties of Josephson junctions in superconductors and superfluids. The 
physical idea consists of considering a double-well trap, with a barrier between the two condensates. If the chemical 
potential in the two traps is different, a flux of atoms is generated. In Fig. ^ we show a simplified scheme. If one 
assumes the barrier between the two wells to be high enough, then Eq. ( |39| ) has two natural solutions, 4>i(r) and 02 (r), 
localized in each potential well, 1 and 2, and having chemical potentials /ii and fi2- A difference between the chemical 
potentials in the two traps can be achieved by filling them with a different number of atoms. The overlap between 
the condensates occurs only in the classically forbidden region, where the wave function is small and nonlinear effects 
due to interactions can be ignored. Thus in this region the linear combination 

0(r,i) = 0i(r)exp(-z^) + (/)2(r)exp(-z^) (147) 



is still a solution of the time dependent equation {p5\). If the two condensates are elongated in the z-direction, the 
current through the barrier can be written as 

I{z,t) = ^ dxdy(^{v,t) — 4>*{v,t)-4>*{v,t) — ^{v,t)] (148) 
2m J \ oz oz J 



Using the wave function ( 147 ) , the current can be easily calculated and takes the typical Josephson form 

I = Iosm[{fii- ^i2)t/h] (149) 

with Iq — (h/m) J dxdy{(f>i(f)'2 ~ 4'24'i)- The calculation of the critical current Iq is a difficult task, since it corresponds 
to a nonlinear 3D tunneling problem. If /ii and fj,2 differ from the average value /i by a small quantity Sfi and 
one treats the motion under the barrier in WKB semiclassical approximation, the estimated current Iq turns out 
to be proportional to exp(— S'o), where 5*0 = /^^ dz[2TO(ycxt(0, 0, z) — /i)/?i^]^/^ and the integral is taken between 



the points 1 and 2 located as in Fig. 35 (Dalfovo, Pitaevskii and Stringari, 1996). Zapata, Sols and Leggett (1998) 
have recently applied the same formalism to realistic 3D configurations, deriving the Josephson current in the form 
I ~ {ksT^ /K) ex-p{—So). Their results suggest that Josephson effects might be indeed observed in experiments. At 
finite temperatures one should also include possible contributions arising from the incoherent fiux of thermally excited 
atoms; this "normal" current is expected to be proportional to In order to observe the Josephson effect one must 
consequently work at low enough temperatures where the system is fully superfluid. It is also worth noticing that the 
geometry of the trapped gases allows one to realize qualitatively new Josephson-type effects, as suggested by Smerzi 
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et al. (1997), Raghavan et al. (1997) and Williams et al. (1998). The propagation of density solitons across regions 
of phase discontinuity in the collision of two condensates has been also considered as an analogue of a Josephson-like 
effect (Reinhardt and Clark, 1997). 

An open problem concerns the possible decoherence mechanisms, which could affect, or even destroy, the phase 
coherence in interference and Josephson-like experiments through phase diffusion processes. The fluctuations of the 
phase can have either a thermal or a quantum origin. Actually, even at T = the phase of the condensate must 
diffuse since having a fixed phase is inconsistent with atom number conservation. Many authors have investigated 
the problem of the quantum diffusion of the phase by describing the system as a coherent superposition of states 
with different N. This yields fluctuations in the chemical potential and hence in the phase of the order parameter. 
Discussions about this effect and on the general problem of phase coherence can be found in: Lewenstein and You 
(1996), Barnett, Burnett and Vaccarro (1996), Wright, Walls and Garrison (1996), Castin and Dalibard (1997), Walhs 
et al. (1997b), Imamoglu, Lewenstein and You (1997), and Dodd et al. (1997b), Javanainen and Wilkens (1998), 
Leggett and Sols (1998) [see also Parkins and Walls (1998) for more discussions and references]. 



C. Collapse and revival of collective oscillations 



In trapped gases one can predict another interesting "mesoscopic" phenomenon having no classical analog, namely 
the collapse and revival of collective excitations. This process should not be confused with the decay of coherence in the 
many-body wave function, which corresponds to the phase diffusion mentioned at the end of the previous section and 
which is also sometimes called "collapse" of the condensate. Conversely, the collapse-revival of collective excitations 
originates from a dephasing of an oscillation due to the quantum fluctuation of the number of quanta. Indeed, a 
classical oscillation of the condensate can be viewed as a coherent superposition of stationary states with different 
numbers of quanta of the oscillator. Fluctuations in the number of quanta cause a dephasing and a consequent 
decrease in the amplitude (collapse) of the oscillation. Since there is no dissipation of energy in this process, the 
oscillation can eventually reappear (revival) after a certain time interval. A schematic picture is shown in Fig. |3^. 
Similar processes of collapse-revival of coherent quantum states have been already observed in atomic Rydberg wave 
packets (Yeazell and Stroud, 1991; Meacher et al., 1991), molecular vibrations (Vrakking et al., 1996) as well as atoms 
and ions interacting with an electromagnetic field (Meekhof et al., 1996; Brune et al., 1996). 

Let E be the energy associated with a classical oscillation of the system induced, for example, by some external 
sinusoidal drive. By classical oscillation we mean that the number n of quanta of os cillation is very large. Let us 
further suppose the frequency lo to be weakly dependent on the amplitude as in Eq. (107). The energy of the oscillation 
is proportional to the square of the amplitude, E cx , and the coefficient of proportionality can be calculated, for 
instance, by solving the time dependent GP equation (|35|). Thus one can rewrite Eq. (107) as 



(150) 



with \k\E <C 1. Now, one can use the semiclassical approximation in order to express the energy E in terms of the 
number of quanta of oscillation, through Tito = (dEn/dn). One finds En/h = nojo + n'^{fiu)QK/2). The wave function 
describing the coherent state of the oscillator can be written in the form ip — Cnipn exp{—iEnt/h). The coefficients 



1 



exp 



(n — n)"^ 
2ri 



(151) 



characterize a coherent gaussian distribution, and n is the average value of quanta, which is supposed to be much larger 
than 1. Given a generic oscillator co-ordinate £,{t), its average over this Gaussian superposition of states can be easily 
estimated; it takes contributions from n n±l transitions and the result is (^(i)) oc J2n I P cos[(wo + fiujQKn)t]. 
For small enough values of t, one can replace the summation over n by an integral and one gets a Gaussian damping, 
or collapse, of the oscillation according to (£_) cx exp[— (^/tc)^], where 



= LUoiEhuJa/2y/^ 



(152) 



defines the collapse time. The periodicity of {^{t)) gives also the revival time, = 2-k / [fiLu'^K) . One finds Tc ~ 
\/{Tfn) Tr <^ Tr- Thcse expressions for the time scales were derived by one of us (Pitaevskii, 1997), and the theory 
of collapse-revival has been also developed by Kuklov et al. (1997) and Graham et al. (1998). The revival time is in 
agreement with the theory of Averbukh and Perelman (1989). 

An explicit estimate of the collapse time (152) can be obtained, for instance, using Gross-Pitacvskii theory within 
the coUisionless hydrodynamic scheme of Sec. [VD. One can solve the equations of motion for the lowest m = 
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and m = 2 mode s in a n axi ally symmetric trap and find the relation between the coefficients k and S{\) entering the 



expansions (150) and (107), respectively (Dalfovo et ai, 1997b). For the m = 2 mode in the first JILA trap, with 
about 5000 rubidium atoms, one finds a collapse time of the order of 5 s, if the relative amplitude is about 20%. This 
time is much larger than the lifetime reported in the experiment by Jin et al. (1996, 1997), which is of the order 
of 100 ms and hence clearly originates from other damping mechanisms. The collapse time would be even longer 
for larger N and this makes the collapse rather difficult to observe. It is not hopeless, however, to observe such an 
effect at very low temperatures where other damping processes, like Landau damping, become much less effective, 
and looking for special values of the anisotropy parameters, where nonlinear effects and frequency shifts are larger. 



as suggested in Sec. IV D 



VII. CONCLUSIONS AND OUTLOOK 

In this paper we have provided an introductory description of the properties of Bose condensed gases confined 
in harmonic traps. The main message emerging from our analysis is that, despite the dilute nature of these gases, 
two-body interactions have crucial consequences on most measurable quantities. This is the combined ef fect of Bose- 



Einstein condensation and of the nonuniform nature of the system. In particular the ground state (Sec. Ill) and the 



dynamic (Sec. IV) properties are affected by two-body forces in an essential way. Interactions can be included using 
fundamental many-body theories for the order parameter, depending on a single interaction parameter, the s-wave 
scattering length. Direct measurements of the density profiles, release energy and collective frequencies have already 
provided very accurate tests of the theoretical predictions. Concerning thermodynamics (Sec. |^ our analysis has 
pointed out the possibility of defining the thermodynamic limit for these nonuniform systems including the effects of 
two-body interactions. Such effects are less important than those of the ground state since at finite temperature the 
system is more dilute. Nevertheless significant corrections to the critical temperature and to the T-dependence of the 



release energy can be predicted and in some cases compared with experiments. In Sec. VI we have discussed possible 
superfluid and coherence phenomena exhibited by trapped Bose gases. This discussion could not be exhaustive because 
the evolution of current research in this field is very rapid. 

The mean-field picture of these interacting Bose gases turns out to be quite accurate in describing most of the 
available experimental results. Deviations from the mean- field predictions are expected to arise from "correlation" 
effects beyond Bogoliubov theory, when the gas parameter n|a|^ is not very small. They can also originate from 
"mesoscopic" effects associated with the fact that the concepts of order parameter and gauge symmetry breaking are 
only approximate in finite systems and, in particular, the fluctuations of the phase are not always negligible. There 
are no experimental evidences so far for these effects, but accurate theoretical predictions concerning both correlation 
and mesoscopic effects might stimulate new important experiments in the future. 

In our review we have been able to cover only part of the huge body of literature which arose after the experimental 
discovery of Bose-Einstein condensation in 1995. We would like to mention here some important issues that we have 
not discussed and that have been recently at the center of significant theoretical and/or experimental research. 

Kinetics of the condensate: An important question, not yet fully understood, is the kinetics of the condensate 
nucleation. The process of condensation of a uniform ideal gas was considered by Semicoz and Tkachev (1995) on 
the basis of the Boltzmann equation. They assumed that the distribution function depends only on the energy of the 
atoms and found that this function exhibits a divergence at zero energy after a finite time interval, corresponding to 
the onset of Bose-Einstein condensation. In previous investigations the mechanism of condensation was predicted to 
occur only asymptotically. The next stage of the process is the growth of the condensate. This was considered recently 
by Gardiner et al. (1997) and Jaksch et al. (1997) on the basis of quantum kinetic master equations. The kinetics of 
the Bose gas near critical conditions for condensation has also been studied by Monte Carlo simulations taking into 
account the Bose statistics under the random phase approximation (Wu, Arimondo and Foot, 1997). In a very recent 
experiment at MIT (Miesner et ai, 1998) the formation and growth of the condensate has been investigated by means 
of imaging techniques. This work has shown clear evidence for a behavior known as "bosonic stimulation" , which 
corresponds to an enhancement of the condensation rate induced by the condensate itself. Explicitly, if N atoms are 
in the condensate the condensation rate is proportional to {N + 1). This gain mechanism is familiar in the physics of 
optical lasers and, in the case of trapped atoms, can lead to matter-wave amplification. Another important question 
recently investigated experimentally is the decay of the trapped gas due to three-body recombination, yielding forma- 
tion of molecules and loss of atoms from the trap (Burt el ai, 1997; Stamper-Kurn et ai, 1998a). The corresponding 
rate turns out to be different for a condensate and a thermal cloud, confirming the theoretical predictions by Kagan, 
Svitsunov and Shlyapnikov (1985). 

Mixtures of condensates: Binary mixtures of condensates can be obtained experimentally by trapping at the same 
time two different atomic species (different isotopes or different alkalis), or two different spin states of the same atoms. 
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The confining potential of the two condensates may be centered at the same point, or not. Ground state and dynamic 
calculations of two interacting condensates have been carried out by Ho and Shenoy (1996b) using the Thomas- Fermi 
approximation and by Esry et al. (1997) by solving the GP equations. Similar calculations have also been done by 
Graham and Walls (1998) and by Pu and Bigelow (1998). Interesting behaviors can be predicted depending on the 
values of the intraspecies and interspecies scattering lengths. Hydrodynamic equations have been recently derived by 
Ho and Shenoy (1998). From the experimental viewpoint, mixtures of this kind have been created and observed at 
JILA (Myatt et al., 1997) with different spin states of *^Rb. In a very recent experiment (Matthews et al., 1998) all 
the atoms have been converted, via two-photon transitions, into a different hyperfine state. The system in the final 
configuration is no longer in equilibrium and will start oscillating. From the analysis of the subsequent oscillations 
(see Fig. ^7|) it has been possible to determine with high precision the ratio of the intraspecies scattering length 
relative to the final and initial states. Using the same apparatus with a mixture of the \F — l,mp ~ —1) and 
\F — 2,mp = 1) spin states, it was also possible to measure the relative phase of two condensates, thus realizing a 
"condensate interferometer" (Hall et al., 1998). Suggestions have been made to use the same mixture of states in 
order to observe nonlinear Josephson-type oscillations (Williams et al., 1998). 

Fermions: The study of degenerate Fermi gases in traps is expected to be an important issue of future research. 
Trapping of fermionic species has been reported for ^Li (McAlexander et al., 1995) and "'"K (Cataliotti et al., 1998). 
Sympathetic cooling of fermions by bosons might yield low temperature regimes overcoming the problem of the 
suppression of coUisional processes exhibited by polarized Fermi gases at low temperature. Degenerate Fermi gases 
behave quite differently from bosonic sytems. Effects of Fermi statistics can be observed in the behavior of the 
release energy below the Fermi temperature; for an ideal gas of N fully polarized atoms, the latter is given by 
kgTp = {6NX)^^^huj±, with the usual definition A — ll>z/^±- Figure ^ shows how the release energy of an ideal 
Fermi gas compares with the corresponding behavior of an ideal Bose gas confined in the same harmonic potential and 
with the same number of atoms. At very low temperatures interacting Fermi gases can undergo a phase transition to 
a superfluid phase. The resulting behavior in the presence of a harmonic trap has been the object of several studies 
[see, for instance, Baranov et al. (1996), Stoof et al. (1996), Houbiers et al. (1997) and references therein]. 

Optical confinement: The recent realization of Bose-Einstein condensation in optical traps (Stamper-Kurn et al, 
1998a) is also expected to open important perspectives. On the one hand one can obtain higher densities, useful, 
for example, to study three-body decay processes and more correlated configurations. On the other hand different 
geometrical configurations can be achieved, like for example quasi ID structures. Finally, by releasing the condition 
of spin polarization imposed by magnetic trapping, this new method of confinement will permit one to study in a 
systematic way the magnetic properties of these gases, including the spinor nature of the order parameter (Ho and 
Shenoy, 1996a; Ho, 1998; Ohmi and Machida, 1998). Spin domains in condensates of sodium, made by three hyperfine 
states of the F — 1 multiplet, have been recently observed by Stenger et al. (1998), who have demostrated the anti- 
ferromagnetic character of the spin-dependent interaction. A further advantage of the optical traps is that they allow 
one to observe Feshbach resonances for strong-field seeking states, as already done by Inouye et al. (1998) with 
sodium. Feshbach resonances are strong variations of the scattering length, induced by an external field, which occur 
when a quasibound molecular state has nearly zero energy and couples resonantly to the free states of the colliding 
atoms. The possibility of tuning the scattering length with external magnetic fields provides new perspectives in the 
manipulation of Bose condensates. 
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FIG. 1. Images of the velocity distribution of rubidium atoms in the experiment by Anderson et al. (1995), taken by 
means of the expansion method. The left frame corresponds to a gas at a temperature just above condensation; the center 
frame, just after the appearance of the condensate; the right frame, after further evaporation leaves a sample of nearly pure 
condensate. The field of view is 200/im x270/im, and corresponds to the distance the atoms have moved in about 1/20 s. The 
color corresponds to the number of atoms at each velocity, with red being the fewest and white being the most. Prom Cornell 
(1996). [note: this is B/W version of reduced quality for e-archive only; the original is a color jpeg file] 



FIG. 2. Collective excitations of a Bose-Einstein condensate. Shown are in-situ repeated phase-contrast images taken of 
a "pure" condensate. The excitations were produced by modulating the magnetic fields which confine the condensate, and 
then letting the condensate evolve freely. Both the center-of-mass and the shape oscillations are visible, and the ratio of their 
oscillation frequencies can be accurately measured. The field of view in the vertical direction is about 620/jm, corresponding 
to a condensate width of the order of 200-300/Lim. The time step is 5 ms per frame. Prom Stamper-Kurn and Ketterle (1998). 
[note: this is B/W version of reduced quality for e-archive only; the original is a color jpeg file] 
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FIG. 3. Density distribution of 80000 sodium atoms in the trap of Hau et al. (1998) as a function of the axial co-ordinate. 
The experimental points correspond to the measured optical density, which is proportional to the column density of the atom 
cloud along the path of the hght beam. The data well agree with the prediction of mean-field theory for interacting atoms (solid 
Conversely, a noninteracting gas in the same trap would have a much sharper Gaussian distribution 



[II 



line) discussed in Sec. 

(dashed line) . The same normalization is used for the three density profiles. The central peak of the Gaussian is found at about 
5500Aim-^ The fi gure points out the role of atom-atom interaction in reducing the central density and enlarging the size of 
the cloud. 



500 




FIG. 4. Column density for 5000 noninteracting bosons in a spherical trap at temperature T — 0.9Tc . The central peak is the 
condensate, superimposed on the broader thermal distribution. Distance and density are in units of a^o and a^^, respectively. 
The density is normalized to the number of atoms. The same curves can be identified with the momentum distribution of the 
condensed and noncondensed particles, provided the abscissa and the ordinate are replaced with pz, in units of a^^ , and the 
momentum distribution, in units of ajj^, respectively. 
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FIG. 5. Condensate fraction as a function of T/T^ . Circles axe the experimental results of Ensher et al. (1996), while the 
dashed line is the law 
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FIG. 6. Condensate fraction vs. temperature for an ideal gas in a trap. The circles correspond to the exact quantum 
calculation for = 1000 atoms in a trap with spherical symmetry and the solid line to the prediction (0). The dashed line 
refers to the thermodynamic limit (|l5|). 
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FIG. 7. Behavior of an ideal gas with A*' = 10^ particles in a highly anisotropic trap: ijj±_ ~ 5.6 x lO'^tJ^, corresponding to 
T^D = iTiD- Solid line: fraction of atoms in the ground state (ux — 0,ny — 0,nz =0), dashed line: fraction of atoms in the 
lowest radial state {rix = 0,ny — 0). 
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FIG. 8. Condensate wave function, at T = 0, obtained by solving numerically the stationary GP equation (js^) in a spherical 
trap and with attractive interaction among the atoms (a < 0). The three solid lines correspond to N\a\/aho = 0.1,0.3,0.5. 
The dashed line is the prediction for the ideal gas. Here the radius, r, is in units of the oscillator length aho and we plot 
{a'^^/Ny^^(j}{r), so that the curves are normalized to 1 [see also Eq. (ptc|)]. 
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FIG. 9. Same as in Fig. ^, but for repulsive interaction (a > 0) and Na/a^^o ~ 1, 10, 100. 
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FIG. 10. Comparison of the release energy as a function of interaction strength from the stationary GP equation (solid line) 
and the experimental measurements (solid circles). Inset shows the expansion of widths of the condensate in the horizontal 
(empty circles) and vertical (crosses) directions against the predictions of the time dependent GP equation (dashed and solid 
lines) for the data point at 10~''A'^i^^^^ — 0.53. Here v is the frequency of the trapping potential and the trapped gas is 
rubidium. From Holland et al. (1997). 
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FIG. 11. Release energy of the condensate as a function of the number of condensed atoms in the MIT trap with sodium 
atoms. For these condensates the initial kinetic energy is negligible and the release energy coincides with the mean-field energy. 
The symbol Uint is here used for the mean-field energy per particle. Triangles: clouds with no visible thermal component. 
Circles: clou ds with both thermal and condensed fractions visible. The solid line is a fit proportional to A'^'^^ (see discussion 
in Sec. |[IID[). From Mewes et al. (1996a). 
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FIG. 12. Energy per particle, in units of ftwho, for atoms in a spherical trap interacting with attractive forces, as a function of 
the effective width w in the Gaussian model of Eqs. (p8|)-(|4S|). Curves are plotted for several values of the parameter Ai'|a|/aho. 
The local minimum disappears at A'^ = A'cr- 
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FIG. 13. Density profile for atoms interacting with repulsive forces in a spherical trap, with TVa/sho = 100. Solid line: 
solution of the stationary GP equation (^). Dashed line: Thomas- Fermi approximation (^o|). In the upper part the atom 
density is plotted in arbitrary units, while the distance from the center of the trap is in units of aho. The classical turning point 
is at ~ 4.31aho. In the lower part the column density for the same system is reported. 
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FIG. 14. Frequency of the lowest collective modes of even parity, m = and m = 2, for rubidium atoms in the JILA trap 
(A = \/8)- The abscissa is the dimensionless parameter Na/a^, with = [?i/(mtJx)]^''^, while the frequency is given in units 
of u)_L. Points are taken from the experimental data of Jin et al. (1996). Solid lines are the predictions of the mean-field 
equations (p6[)-(|6^) [see, for instance, Edwards et al. (1996c); Esry, (1997); You, Hoston and Lewenstein (1997)]. Dashed lines 
are the asymptotic results for Na/a± <x (Stringari, 1996b), as discussed in Sec. IVB. 
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FIG. 15. Speed of sound, c, versus condensate peak density, n(0), for waves propagating along the axial direction in the 
cigar-shaped condensate at MIT. The experimental points are compared with the theoretical prediction c = [gn{0)/2m]^'''^ 
(solid line). From Andrews, Kurn et al. (1997). 
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FIG. 16. Oscillations of the axial width of the condensate in the cigar-shaped trap at MIT. The excited collective motion 
is the low-lying m — mode. Oscillations are shown at low (a) and high (b) temperature. Points show the axial width 
determined from fits to phase-contrast images, similar to the ones in Fig. y. Lines are fits to a damped sinusoidal oscillation. 
From Stamper-Kurn et al, (1998c). 
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FIG. 17. Frequencies of the monopole and quadrupole excitations of a condensate in a spherical trap as a function of the 
parameter Na/aho, for positive and negative values of a. The solid line for the monopole mode is obtained from the ratio 
(mi/m_i)^/^, as in Eq. (p5|). For the quadrupole mode it corresponds to the ratio (ms/mi)^^^ , as in Eq. (^^. Circles and 
squares are the eigenenergies of the linearized time dependent GP equation (|66|)-(|6^. 
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FIG. 18. Aspect ratio, R±/Z, of a freely expanding condensate as a function of time. The experimental points in part (a) 
correspond to *^Rb atoms initially confined in a trap with A = 0.099 (Ernst et al, 1998b). The points in part (b) are measure- 
ments on sod ium atoms, initially in a trap with A = 0.065 (Stamper-Kurn and Ketterle, 1998). The solid lines are obtained by 
solving Eqs. (lOJ), which are equivalent to the time dependent GP equat ion i n the Thomas- Fermi approximation. The dashed 
lines correspond to the A ^ 1 limit of the same equations, that is to Eqs. ( 104 )-( |l05|) , and are almost indistinguishable from the 
solid lines. The dot-dashed lines are the predictions for noninteracting atoms. Theoretical curves have no fitting parameters. 
In part (a), they have been corrected to include the effect of the observation angle, as explained by Ernst et al. (1998b). 
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FIG. 19. Frequency of the low-lying m — mode measured at MIT (Stamper-Kurn et ai, 19 98, S tamper-Kurn and Ketterle, 
1998) as a function of the amphtude of the osciUation. The solid line is the prediction of Eqs. (102). 
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FIG. 20. Excitation spectrum of 10000 atoms of *^Rb in a spherical trap with aho = 0.791 /im. The eigenenergies of 
the Hnearized time dependent GP equations (|66[)-(p7|) are represented by thick solid bars. Dashed bars correspond to the 
single-particle spectrum of Hamiltonian ( |l09| ). The thin orizontal line is the chemical potential, /i — 8.41 in units of hcoi^o, 
which is fixed by the solution of the stationary GP equation (|39[). 



60 



7500 



5000 



2500 




FIG. 21. Number of states N{e) vs. energy. Solid circles are obtained by counting the Bogoliubov-type states in the spectrum 
shown in Fig. ^ (thick solid bars). Open circles corresponds to counting the single-particle state s in the same figure (dashed 
bars). Both calculations are compared with the predictions of the semiclassical approximation (112) (solid line), as well as 



with the ones of the noninteracting harmonic oscillator (dashed line) and of the coUisionless hydrodynamic equations in the 
Thomas- Fermi regime (dot-dashed line). Chemical potential is /i = 8.41 in this scale. 
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FIG. 23. Energy per particle as a function of T/T^ for r; = 0.4. The solid line refers to the perturbati ye ex p ansi on (123); 
the dashed line is the result of the self-consistent calculation based on the Popov approximation [see Eqs. (127)-(12£)]. 
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FIG. 24. Theoretical predictions for the condensate fraction vs. T/T^ for interacting (squares) and noninteracting (circles) 
particles in two different traps. Concerning interacting particles, we show here the results obtained from a self-consistent 
mean-field calculation, within Popov approximation, for N = 5 x 10* rubidium atoms in a trap with a/aho ~ 5.4 x 10~^ and 
A — \/8 (open squares) and for N — 5 x 10^ sodium atoms in a trap with a/aho = 1-7 x 10"'^ and A = 0.05 (solid squares). 
The numerical results are compared with the prediction of the scaling limit for = 0.4 (solid line). Open and solid circles 
correspond to A'' = 5 x lO** and N = 5 x 10^ noninteracting particles, respectively, in the same two traps as the corresponding 
open and solid squares. The dashed line is the 1 — curve of the noninteracting model in the thermodynamic limit. 
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FIG. 25. Axial profiles of a cloud of sodium atoms. The thick solid lines are two of the profiles already shown in Fig. 03, 
obtained at MIT (Stamper-Kurn and Ketterle, 1998) from phase-contrast images at different temperat ures: T — 0.7^K 
(lowermost), T = 1.2^K (uppermost). Dot-dashed lines are theoretical predictions obtained from Eqs. (127)-(12!:), using 
A'^ and T as fitting parameters: the lower curve corresponds to N = 1.4 x lO'^ and T = 0.8 /iK, and the upper one to 
iV = 2.3 X 10^ and T = 1.1 fiK. In both cases, a difference in temperature of about 10 % between the experimental estimate 
and the result of the fit is consistent with the experimental uncertainty. 
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FIG. 26. Chemical potential as a function of T/T^ in the thermodynamic limit. Solid line: rj = 0.4, dashed line: r; — 0.6. 
The dotted line refers to the non-interacting model (77 = 0). 
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FIG. 27. Condensate fraction as a function of T/T^ in the thermodynamic limit. As in Fig. 26| the curves are the theoretical 
predictions for rj = 0.4 (solid line), rj — 0.6 (dashed line) and the noninteracting case r] — (dotted line). Open circles are the 
experimental data by Ensher et al. (1996), corresponding to r) in the range 0.39 — 0.45. Solid circles with error bars are the 
path integral Monte Carlo results by Krauth (1997), with r) — 0.35. 
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FIG. 28. Release energy as a function of T/T^ in the thermodynamic limit. The curves refer to the same values of r; as in 
Fig. Circles are the experimental data by Ensher et al. (1996). 
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FIG. 29. Temperature dependent damping rates, 7. rn easured for the m — (triangles) and m — 2 (circles) modes (Jin et 
ai, 1997). The solid line is the theoretical estimate ( [130| ) where we have used the Thom as- F ermi approximation for the local 
sound velocity, c(T) = {fj,{T)/m)^^^ , calculated in the center of the trap, and expression (12C) for the temperature dependence 
of the chemical potential. 



0.25 



0.2 





0.05 




1.7 1.8 1.9 
Re(a;)/ cj^ 



2.1 



FIG. 30. The imaginary part of lj against the real part, as given by the interpolation formula (132), in the case of the 
low-lying m = mode observed at MIT (Stamper- Kur n et at, 1998c). In the collisional hydrodynamic regime the frequency of 
the mode is given by ujhd = (12/5) 



1/2,, 



see Eq. (131)], while in the coUisionless regime it is given by the noninteracting value 



uc ~ 2ujz. In both cases the motion is undamped (Im(a;) = 0). Stamper-Kurn et al. (1998c) measured a frequency of about 
30 Hz with a damping rate of about 20 s~^; the corresponding values Re(Lj) ~ 1.78aJz and —lm{Lj) ~ O.WuJz are represented 
by the solid circle. The theoretical curve near this point corresponds to collision time such that Re(a;)r ~ 1. 
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FIG. 31. Condensate with a quantized vortex along the z-axis. The order parameter, 4>v{x, 0, 0), is plotted in the case ol 10* 
rubidium atoms confined in a spherical trap with Oho = 0.791 /im. Distances are in units of the oscillator length aho and the 
curves correspond to {a\^/NY^'^4'v{r), so that they are normalized to 1 wh en in normalized to A'^. The dot-dashed line is 

with «: = 0; the solid line is the profile of a vortex 
In the inset, the contour plot for the density 



the solution of the stationary GP equation (^9|) , or equivalently of Eq. ( 13 
with K = \, from ([l35|) ; the dashed line is the noninteracting wave function (136) 



in the a;z;-plane, n{x, 0, z) — \<f>^{x, 0, z)\ , is given. Luminosity is proportional to density, the white area being the most dense, 
[note: this is version of reduced quality for e-archive only] 
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FIG. 32. Critical angular velocity, in units of uj_l, for the formation of a k = 1 vortex in a spherical trap with atoms of 
''Rb and Oho ~ 0.791 /im. 
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FIG. 33. Moment of inertia, O, divided by its rigid value 6rig, as a function of T/T^. Solid line: interacting gas in the 
thermodynamic limit with q = 0.4. The dashed and dot-dashed lines are the predictions for 5 x lO'^ and 5 x 10* noninteracting 
particles, respectively, in a spherical trap. 
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FIG. 34. Density pattern for the interference of two expanding and overlapping condensates, (a) Theory by Rohrl et al. 
(1997), based on the solution of the time dependent GP equation, (b) Experimental data by Andrews, Townsend et al, (1997). 
(c) Theory including the effect of finite experimental resolution. 
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FIG. 36. Schematic picture for collapse-revival of collective excitations. The quantity ^ is a generic oscillator co-ordinate 
and the symbol (^) means an average over different replica of the system prepared in the same conditions. 



68 




time (units of ro ^) time (units of w ^) 

FIG. 37. Oscillation in the width of the cloud in both the axial and radial direction due to the instantaneous change in 
scattering length in the experiment by Matthews etal. (1998). Time is in units of uj^^ = '^±^ = 9.4 ms. The solid lines are 
the time dependent widths calculated using Eqs. ( [l02| ), with only the amplitude of the oscillation and the initial size as free 
parameters. 
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